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In the past decade, Resonant Inelastic X-ray Scattering (RIXS) has made remarkable 
progress as a spectroscopic technique. This is a direct result of the availability of high- 
brilliance synchrotron X-ray radiation sources and of advanced photon detection instru- 
mentation. The technique's unique capability to probe elementary excitations in com- 
plex materials by measuring their energy-, momentum-, and polarization-dependence 
has brought RIXS to the forefront of experimental photon science. We review both 
the experimental and theoretical RIXS investigations of the past decade, focusing on 
those determining the low-energy charge, spin, orbital and lattice excitations of solids. 
We present the fundamentals of RIXS as an experimental method and then review the 
theoretical state of affairs, its recent developments and discuss the different (approxi- 
mate) methods to compute the dynamical RIXS response. The last decade's body of 
experimental RIXS data and its interpretation is surveyed, with an emphasis on RIXS 
studies of correlated electron systems, especially transition metal compounds. Finally, 
we discuss the promise that RIXS holds for the near future, particularly in view of the 
advent of x-ray laser photon sources. 
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I. INTRODUCTION 




particle momentum/ft (A ') 



FIG. 1 (Kinetic) energy and momentum carried by the differ- 
ent elementary particles that are often used for inelastic scat- 
tering experiments. The scattering phase space -the range of 
energies and momenta that can be transferred in a scattering 
event- of x-rays is indicated in blue, electrons in brown and 
neutrons in red. 



Resonant Inelastic X-ray Scattering (RIXS) is a fast 
developing experimental technique in which one scatters 
x-ray photons inelastically off matter. It is a photon-in 
— photon-out spectroscopy for which one can, in prin- 
ciple, measure the energy, momentum, and polarization 
change of the scattered photon. The change in energy, 
momentum, and polarization of the photon are trans- 
ferred to intrinsic excitations of the material under study 
and thus RIXS provides information about those excita- 
tions. RIXS is a resonant technique in which the energy 
of the incident photon is chosen such that it coincides 
with, and hence resonates with, one of the atomic x-ray 
transitions of the system. The resonance can greatly en- 
hance the inelastic scattering cross section, sometimes 
by many orders of magnitude, and offers a unique way to 
probe charge, magnetic, and orbital degrees of freedom 
on selective atomic sites in a crystal. Early experimen- 



tal work, and some more recent reviews include (Ban 



nett and Freund 1975 Blume 1985 Eisenberger et al. 



1976a|b| IKotani and Shin[ |2001| |Schaike[ |2007[ |Sp~ 
19741). 



1. Features of RIXS as an experimental method 

Compared to other scattering techniques, RIXS has a 
number of unique features: it covers a huge scattering 
phase-space, is polarization dependent, element and 
orbital specific, bulk sensitive and requires only small 
sample volumes. We briefly illustrate these features 
below and discuss them more extensively in the sections 
to follow. 



1. RIXS exploits both the energy and momentum 
dependence of the photon scattering cross-section. 
Comparing the energies of a neutron, electron, and 
photon, each with a wavelength on the order of the 
relevant length scale in a solid, i.e. the interatomic 
lattice spacing, which is on the order of a few Angstroms, 
it is obvious that an x-ray photon has much more energy 
than an equivalent neutron or electron, see Fig. [T] 
The scattering phase space (the range of energies and 
momenta that can be transferred in a scattering event) 
available to x-rays is therefore correspondingly larger 
and is in fact without equal. For instance, unlike 
photon scattering experiments with visible or infrared 
light, RIXS can probe the full dispersion of low energy 
excitations in solids. 

2. RIXS can utilize the polarization of the photon: 
the nature of the excitations created in the material 
can be disentangled through polarization analysis of 
the incident and scattered photons, which allows one, 
through the use of various selection rules, to characterize 
the symmetry and nature of the excitations. To date, 
very few experimental facilities allow the polarization of 



the scattered photon to be measured ( jBraicovich et al. 



2007b Ishii et al. 2010), though the incident photon 



polarization is frequently varied. It is important to note 
that a polarization change of a photon is necessarily re- 
lated to an angular momentum change. Conservation of 
angular momentum means that any angular momentum 
lost by the scattered photons has been transferred to 
elementary excitations in the solid. 

3. RIXS is element and orbital specific: Chemical 
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FIG. 2 Energy of K, L\, L3, Mi and M5 x-ray absorption 
edges as a function of atomic number Z. X-ray energies below 
1 keV are referred to as soft, above as hard. 

sensitivity arises by tuning the incident photon energy to 
specific atomic transitions of the different types of atoms 
in a material. Such transitions are called absorption 
edges. RIXS can even differentiate between the same 
chemical element at sites with inequivalent chemical 
bondings, with different valencies or at inequivalent 
crystallographic positions if the absorption edges in 
these cases are distinguishable. In addition, the type of 
information that may be gleaned about the electronic 
excitations can be varied by tuning to different x-ray 
edges of the same chemical element (e.g., K edge for 
exciting Is core electrons, L edge for electrons in the 
n = 2 shell or M edge for n = 3 electrons), since the 
photon excites different core-electrons into different 
valence orbitals at each edge. The energies of these 
edges are shown in Fig. [2j 

4. RIXS is hulk sensitive: the penetration depth 
of resonant x-ray photons is material and scattering 
geometry-specific, but typically is on the order of a few 
jim. for photons of 10 keV in the hard x-ray regime, and 
on the order of 0.1 /im for photons of 1 kcV in the soft 
x-ray regime. 

5. RIXS needs only small sample volumes: the 
photon-matter interaction is relatively strong, compared 
to for instance the neutron-matter interaction strength. 
In addition, photon sources deliver many orders of 
magnitude more particles per second, in a much smaller 
spot, than do neutron sources. These facts make RIXS 
possible on very small volume samples, thin films, 
surfaces and nano-objects, in addition to bulk single 
crystal or powder samples. 

In principle RIXS can probe a very broad class of in- 
trinsic excitations of the system under study - as long 
as these excitations are overall charge neutral. This con- 
straint arises from the fact that in RIXS the scattered 



photons do not add or remove charge from the system un- 
der study. In principle then, RIXS has a finite cross sec- 
tion for probing the energy, momentum and polarization 
dependence of, for instance, the electron-hole continuum 
and excitons in band metals and semiconductors, charge 
transfer and cid-excitations in strongly correlated mate- 
rials, lattice excitations and so on. In addition magnetic 
excitations are also symmetry- allowed in RIXS, because 
the orbital angular momentum that the photons carry 
can in principle be transferred to the electron's spin an- 
gular moment. This versatility of RIXS is an advantage 
but at the same time a complicating factor, because dif- 
ferent types of excitations will generally be present in a 
single RIXS spectrum. An objective of this Review is to 
help disentangling these where needed. 

The generic advantages of the RIXS technique listed 
above perhaps raise the question as to why this spec- 
troscopic technique is not as widely used as, say, angle- 
resolved photoemission (ARPES) or neutron scattering. 
The main limitation is that the RIXS process is photon- 
hungry - that is, it requires a substantial incident pho- 
ton fiux to obtain enough scattered photons to collect 
spectra with a high enough resolution in energy and mo- 
mentum in a reasonable time. With a required resolving 
power (defined as the incident photon energy divided by 
the energy resolution) of four orders of magnitude, RIXS 
has been a real challenge. Up until a few years ago this 
has limited RIXS experiments to measuring energy losses 
on the order of half an electron volt or greater. Thus 
neutron scattering and ARPES offered a more direct ex- 
amination of the low energy excitations near the Fermi 
level. However, recent progress in RIXS instrumentation 
has been dramatic and this situation is now changing. 
One of the purposes of this review is to summarize this 
progress which is beginning to elevate RIXS into an im- 
portant condensed matter physics tool for probing ele- 
mentary excitations in solids. 



2. Progress of RIXS in the last decades 

As discussed above, the generic features of RIXS make 
it, in principle, an attractive technique to study the in- 
trinsic momentum dependent, low energy response of a 
material. But there are of course practical limitations. 
The most critical of these is the energy resolution, which 
is determined both by the availability of the instrumen- 
tation necessary to energy resolve the photons, and by 
the availability of tunable photon sources of sufficient in- 
tensity. 

In order to tune the incident photon energy to a par- 
ticular edge, a tunable x-ray photon source is essential. 
This can be achieved with synchrotron radiation sources 
and their increase in brilliance over the past decades has 
been many orders of magnitudes in the 10^-10^ eV x- 
ray regime. The next generation photon sources include 
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Hard x-ray regime: Cu K-edge Soft x-ray regime: Cu L-edge 
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FIG. 3 Development of RIXS resolution, count rate and re- 
sulting orders of magnitude increase in the figure of merit 
(count rate/resolution) in the hard x-ray regime, as measured 
by the 6 eV feature in CuGeOs at the Cu K-edge (8990 eV). 



x-ray free electron lasers (FELs), which are coming on 
line at the time of writing. The peak brilliance of these 
sources is again orders of magnitude larger than that of 
the third generation synchrotrons and it is likely that 
these sources will provide further advances, particularly 
for time-resolved experiments. Such future developments 
will be touched upon in Section [yl] 

This vast increase in photon flux has been matched by 
advances in the RIXS instrumentation: the monochroma- 
tors, analyzers and spectrometers. The resulting increase 
in resolution of RIXS experiments over time, as measured 
for instance at the hard x-ray Cu K-edges and soft x-ray 
Cu L3, is shown in Figs. [3] and |4j respectively. We will 
briefly discuss the current state-of-the-art in Sec. |ll] 

In concert with the great progress in the RIXS ex- 
periments over the past decades, there has been a simi- 
lar rapid advance in the theoretical understanding of the 
scattering process and of the dynamic correlation func- 
tions that the technique probes. Taken together, the the- 
oretical and experimental advances have driven an enor- 
mous increase in the number of RlXS-related publica- 
tions and citations to these publications, as shown in 
Fig. [5] 

It seems likely that this strong growth will continue. 
First, because of the ongoing push to better energy reso- 
lutions. Second, and perhaps more importantly, because 
there are a multitude of different x-ray absorption edges, 
in particular for the heavier elements in the periodic ta- 
ble, and each one of these can, in principle, be exploited 
for RIXS measurements. The bulk of RIXS data so far 
has been collected at transition metal and oxygen edges. 
This is motivated by the intense scientific interest in 
strongly correlated transition metal oxides such as the 
high-Tc cuprate superconductors and the colossal magne- 



FIG. 4 Progress in soft x-ray RIXS resolution at the Cu L- 
edge at 931 eV (a) ( [Ichikawa et al.\ |1996|), BLBB @ Pho- 
ton Factory (b) 1511-3 @ MAX II (|Duda et aZ.[ |2000b), (c) 
AXES Q IDOS, ESR FUGhiringhem ei al.\ |2004p (d) AXES 
@ ID08, ESRF(|B raicovich et ai, 2009Tr(ejSAXES @ SLS 
( Ghiringhelli et al. i^2010^ . Courtesy of G. Ghiringhelli and L. 
Braicovich. 
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FIG. 5 Left: Number of publications on RIXS in five year 
periods. Right: Number of citations to RIXS publications in 
five year periods (Source: ISI Web of Knowledge) . 



toresistance manganites. This focus on transition metal 
oxides is an accident of history. It has been very ben- 
eficial to the field, driving advances in instrumentation 
and theory at the relevant edges, but there is clearly a 
huge potential for growth as interest moves on to other 
materials and other fields. 



3. Probing Elementary Excitations with RIXS 

The elementary excitations of a material determine 
many of its important physical properties, including 
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transport properties and its response to external per- 
turbations. Understanding the excitation spectrum of 
a system is understanding the system. 

In this respect strongly correlated electron materials, 
e.g. transition metal oxides, are of special interest be- 
cause the low energy electronic properties are determined 
by high energy electron-electron interactions (energies 
on the order of eV's). From these strong interactions 
and correlations a set of quantum many-body problems 
emerge, the understanding of which lies at the heart of 
present day condensed matter physics. Most often this 
many-body physics is captured in model Hamiltonians, 
the exact parameters of which must be determined ex- 
perimentally. RIXS, along with other spectroscopic tech- 
niques, can play an important role there, even if it is as 
a spectroscopic technique applicable to many other ma- 
terials and of course not limited to correlated systems. 

In the following, we discuss the relevant excitation en- 
ergy and momentum scale on which RIXS can probe the 
excitation spectrum of a solid. We then briefly introduce 
the kinds of elementary excitations that are accessible to 
RIXS. 



a. Excitation Energy and Momentum Scale As is shown in 
Fig. [6j the elementary excitation spectrum in solids spans 
the range from plasmons and charge transfer excitations 
at a few eV, determining for instance optical properties, 
through excitons, dd excitations and magnons down to 
phonons at the meV scale. In principle, RIXS can mea- 
sure the momentum-dependence of the excitation energy 
of all these modes - i.e. their dispersion - because the 
photon transfers momentum as well as energy to the ma- 
terial under study. 

This is unusual if one is accustomed to optical light 



scattering, such as Raman scattering (Devereaux and 



Hackl 2007). Photons in the visible range of the spec- 



trum with an energy of a few eV carry negligible momen- 
tum compared to the quasi-momentum of the elementary 
excitations of a solid (fig. 1). A photon of 2 eV has a mo- 
mentum of roughly hq ^ 10^^^ kg m/s, or a wavevector 
q = 10~^ whereas elementary excitations in a crys- 
tal with a lattice constant of say 3 A have wavevectors 
up to g = 27r/3 « 2 A~^. On this scale optical light scat- 
tering is in essence a zero momentum probe. To measure 
the dispersion of elementary excitations for momenta in 
a sizable portion of a typical Brillouin zone, x-rays with 
energy on the order of 1 keV or more are needed, corre- 
sponding to for instance the Cu L-edge. 



b. Overview of elementary excitations In this paragraph 
we briefly discuss the different elementary excitations ac- 
cessible to RIXS ~ a detailed review is provided in Chap- 
ter IVl 
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FIG. 6 Different elementary excitations in condensed mat- 
ter systems and their approximate energy scales in strongly 
correlated electron materials such as transition metal oxides. 



Plasmons. Collective density oscillations of an electron 
gas are referred to as plasmons. They can be observed by 
inelastic x-ray scattering (IXS) or by optical probes since 
they occur at finite energy for q=0. Plasmon-like excita- 



tions were also observed early on in RIXS (Isaacs et al. 



1996), but their resonant enhancement with respect to 
IXS is weak, and little work has been done since, see 
Section |V31 

Charge-transfer excitations. Charge transport in a 
condensed matter system is determined by the energet- 
ics of moving electrons from one site to another. In 
a transition-metal oxide, there are two relevant energy 
scales for this process. The first is the energy associated 
with an electron hopping from a ligand site to a metal 
site. This is known as the charge transfer energy. A, 
where A = E{d'^~^^L) — E{d'^), and L represents a hole 
on the ligand site. The second energy scale is the energy, 
U, associated with moving a d electron from one metal 
site to another where U = E{d"-+'^) + E{d"-'^) - 2E{d"). 
Strongly correlated insulators may be classified by which 
of these two energies is the larger (Zaanen et al. 1985| ). 
If f/ > A, then the gap is of the charge transfer type 
and the system is said to be a charge-transfer insulator. 
Conversely, ii U < A, then the gap is controlled by the 
d — d Coulomb energy and the system is said to be a 
Mott-Hubbard insulator. 

The bulk of the interesting transition metal oxide com- 
pounds, including the cuprates, nickelates and mangan- 
ites are all in the charge transfer limit. This means the 
lowest lying excitations across the optical gap are charge 
transfer excitations and therefore these are of central im- 
portance in these materials. Key questions include the 
size of the gap (typically on the order of a few eV) and 
the nature of the excitations - do they form bound ex- 
citon states? Are these localized or can they propagate 
through the lattice? What are their lifetimes, symme- 
tries and temperature dependence, etc. While some stud- 
ies have been performed using other techniques, notably 
EELS and optical conductivity measurements, RIXS of- 
fers a powerful probe for many of these questions and has 
been applied extensively. We review the experimental re- 
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suits and their theoretical interpretation in section V.B 



Crystal Field and Orbital excitations. Many strongly 
correlated systems exhibit an orbital degree of freedom, 
that is, the valence electrons can occupy different sets of 
orbitals. Orbitally active ions are also magnetic: they 
have a partially filled outer shell. This orbital degree 
of freedom determines many physical properties of the 
solid, both directly, and also indirectly because the or- 
bitals couple to other degrees of freedom. For instance, 
the orbital's charge distribution couples to the lattice, 
and according to the Goodenough-Kanamori rules for su- 
perexchange the orbital order also determines the spin- 
spin interactions. The nature of the orbital degree of 
freedom - the orbital ground state and its excitations 
are an important aspect of strongly correlated systems. 

In many Mott insulators this orbital physics is gov- 
erned by the crystal field: the levels of the orbitally active 
ion are split and the orbital ground state is uniquely de- 
termined by local, single-ion considerations. The orbital 
excitations from this ground state are transitions between 
the crystal field levels. Crystal field transitions between 
different d orbitals are called dd excitations. Such excita- 
tions are currently routinely seen by RIXS and are now 
well understood. 

In other cases the crystal field does not split the levels 
of the outer shell very much, leaving an orbital (quasi- 
)degeneracy in the ground state. This local low-energy 
degree of freedom can couple to orbital degrees of freedom 
on neighboring sites by superexchange processes, and in 
this way collective orbital excitations can emerge. The 
quanta of these collective modes are called orbitons, in 
analogy to spin waves and magnons. Definitive proof 
of the existence of orbitons remains elusive - RIXS 
is contributing significantly to the search for orbitons. 
The distinctive advantages of RIXS over other spectro- 
scopies when probing orbital excitations raise promising 
prospects, see Section [V.D.3| 

Magnetic excitations. Magnetism and long-range mag- 
netic ordering are arguably the best known and most 
studied consequences of the electron-electron interactions 
in solids. When usual magnetic order sets in, be it either 
of ferro-, ferri- or antiferromagnetic type, the global spin 
rotation symmetry in the material is broken. As a re- 
sult characteristic collective magnetic excitations emerge. 
The resulting low energy quasiparticles, the magnons, 
and the interactions between them, determine all low 
temperature magnetic properties. Magnon energies can 
extend up to ~ 0.3 eV (e.g. in cuprates) and their mo- 
menta up to ~ 1 A~^. Recently magnon dispersions have 
been measured for the first time at the Cu L-edge on thin 
films of La2Cu04 (Braicovich et al. 2010). In K-edge 



RIXS bi-magnon excitations and their dispersions have 



also been observed ( Hill et al. 2008 1 



spin system, or by the frustration of magnetic interac- 
tions between the spins, can result in the formation of 
spin-liquid ground states. Spin liquids potentially have 
elusive properties such as high-temperature superconduc- 
tivity or topological order, which one is only beginning to 
explore and understand. Some of the more exotic mag- 
netic excitations that emerge from these ground states, 
such as spinons and triplons can also be observed by 
RIXS dSchlappa "al| [2009 ) . In Section 



we review 



A melting of the long-range ordering, for instance 
through an increase in quantum fiuctuations as a result of 
the introduction of mobile charge carriers in a localized 



both the experimental and theoretical status of magnetic 
RIXS in detail. 



A. The RIXS Process 

The microscopic picture of the resonant inelastic x- 
ray scattering process is most easily explained in terms 
of an example. We will choose a copper-oxide mate- 
rial as a typical example, but it should be stressed once 
more that the focus of RIXS on transition metal oxides is 
rather an accident of history. In a copper-oxide material, 
one can tune the incoming photon energy to resonate 
with the copper K, L or M absorption edges, where 
in each case the incident photon promotes a different 
type of core electron into an empty valence shell, see 
Figs. [7|and [sj The electronic configuration of Cu^+ is 
ls'^2s'^2p^5s'^3p^3d^, with the partially filled 3d valence 
shell characteristic of transition metal ions. The copper 
K-edge transition Is — 4p, is around 9000 eV and in 
the hard x-ray regime. The L2,3-edge 2p — > 3d 900 
eV) and M2,3-edge 3p — )■ 3d (~ 80 eV) are soft x-ray 
transitions. Alternatively, by tuning to the Oxygen K- 
edge, one can choose to promote an O Is to an empty 2p 
valence state, which takes ^ 500 eV. 

After absorbing a soft or hard x-ray photon, the sys- 
tem is in a highly energetic, unstable state: a hole deep 
in the electronic core is present. The system quickly de- 
cays from this intermediate state, typically within 1-2 
femtoseconds. Decay is possible in a number of ways, 
for instance via an Auger process, where an electron fills 
the core-hole while simultaneously emitting another elec- 
tron. This non-radiative decay channel is not relevant for 
RIXS, which instead is governed by fluorescent decay, in 
which the empty core-state is filled by an electron and at 
the same time a photon is emitted. 

There are two different scattering mechanisms by 
which the energy and momentum of the emitted photon 
can change from the incident one. These are known as 
direct and indirect RIXS. The distinction between these 
two is discussed below. 



1. Direct and Indirect RIXS 

Resonant inelastic x-ray scattering processes are clas- 



sified as either direct or indirect (Van den Brink and 
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FIG. 7 In a direct RIXS process the incoming x-rays excite an 
electron from a deep- lying core level into the empty valence. 
The empty core state is then filled by an electron from the 
occupied states under the emission of an x-ray. This RIXS 
process creates a valence excitation with momentum Kk' — Kk 
and energy hw]^ ~ hwy^i . 



FIG. 8 In an indirect RIXS process, an electron is excited 
from a deep-lying core level into the valence shell. Excitations 
are created through the Coulomb interaction Uc between the 
core hole (and in some cases the excited electron) and the 
valence electrons. 



van Veenendaal 2005 2006). This distinction is use- 



ful because the cross-sections for each are quite differ- 
ent. When direct scattering is allowed, it is the domi- 
nant inelastic scattering channel, with indirect processes 
contributing only in higher order. In contrast, for the 
large class of experiments for which direct scattering is 
forbidden, RIXS relies exclusively on indirect scattering 
channels. 



lence and conduction states directly. Although the direct 
transitions into the valence shell dominate the spectral 
line shape, the spectral weight can be affected by inter- 
actions in the intermediate-state driven by, for example, 
the strong core-hole potential. 



a. Direct RIXS For direct RIXS, the incoming photon 
promotes a core-electron to an empty valence band state, 
see Fig. [7] Subsequently an electron from a different state 
in the valence band decays and annihilates the core-hole. 

The net result is a final state with an electron-hole 
excitation, since an electron was created in an empty 
valence band state and a hole in the filled valence band, 
see Fig. [7] The electron-hole excitation can propagate 
through the material, carrying momentum hq and energy 
huj. Momentum and energy conservation require that 
q = k' — k and lu = — Wk, where hk. (/ik') and hiO]^ 
(^Wk') are the momentum and energy of the incoming 
(outgoing) photon, respectively. 

For direct RIXS to occur, both photoelectric transi- 
tions - the initial one from core to valence state and 
succeeding one from conduction state to fill the core hole 
- must be possible. These transitions can for instance be 
an initial dipolar transition of Is — >■ 2p followed by the 
decay of another electron in the 2p band from 2p — > Is, 
in for example wide-band gap insulators. This happens 
for instance at the K-edge of oxygen, carbon and silicon. 
At transition metal i-edges, dipole transitions give rise 
to direct RIXS via 2p — > 3d absorption and subsequent 
3d — )■ 2p decay. In all these cases, RIXS probes the va- 



b. Indirect RIXS The indirect RIXS process is slightly 
more complicated. For pure indirect RIXS to occur, pho- 
toelectric transitions from the core-state to conduction- 
band states must be weak. Instead, the incoming pho- 
ton promotes a core-electron into an empty state sev- 
eral electron volts above the Fermi level. Subsequently 
the electron from this same state decays again to fill the 
core-hole, see Fig. [8] The most studied example is RIXS 
at the transition- metal i^-edges (Is — > 4p). Obviously, 
in the absence of any additional interaction, no inelastic 
scattering would be observed. But in the intermediate 
state a core-hole is present, which exerts a strong po- 
tential on the 3c? valence electrons, that therefore tend 
to screen the core-hole. The core-hole potential scatters 
these valence electrons, thereby creating electron-hole ex- 
citations in the valence band. After the 4p — >■ Is decay, 
the electron-hole excitations are then left behind in the 
system. 

Indirect RIXS is thus due to shakeup excitations cre- 
ated by the intermediate state core-hole. The fact that 
close to the absorption edge the Is core-hole and 4p elec- 
tron bind together to form an exciton does not change 
this picture conceptually - in this case, one may think of 
the valence electrons as scattering off this exciton. 
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B. Relation to other spectroscopies 

1. RIXS, (N)IXS, XAS, XES and EELS 

Resonant inelastic x-ray scattering is sometimes also 
referred to as resonant x-ray emission spectroscopy, res- 
onant x-ray fluorescence spectroscopy, and resonant Ra- 
man spectroscopy. These different names all indicate es- 
sentially the same process. Throughout this review, we 
will only use the term resonant inelastic x-ray scattering 
(RIXS). 

RIXS is a two-step process involving a photoelectric 
transition from a deep-lying core level into the valence 
shell, with the material absorbing an x-ray photon. After 
evolving for a very short period in this highly-excited 
state, the system decays under the emission of an x-ray 
photon. 

The term inelastic x-ray scattering (IXS) is reserved 
for processes where the photon scatters inelastically by 
interacting with the charge density of the system. IXS 
measures the dynamic charge structure factor Sq{uj). To 
emphasize the difference with RIXS, IXS is sometimes 
also referred to as non-resonant IXS (NIXS). 

X-ray absorption spectroscopy (XAS) is the first step 
of the RIXS process and thus RIXS is related to it. A 
common way to measure XAS is to study the decay prod- 
ucts of the core hole that the x-ray has created, either by 
measuring the electron yield from a variety of Auger and 
higher-order processes (known as electron yield) or by 
measuring the radiative decay (fluorescence yield). The 
total fluorescence yield corresponds approximately to the 
integration of all possible RIXS processes. 

The term X-ray emission spectroscopy (XES) is gener- 
ally used for the situation where one studies the radiative 
decay of a core hole created by X-rays or electrons with 
an incident energy far away from an absorption edge. 
However, even with the incident energy not tuned to a 
resonance, its theoretical description is much closer to 
RIXS than to IXS, with the modification that the ini- 
tial excitation is not into the valence shell, but into the 
vacuum continuum at much higher energy. 

Electron-energy-loss spectroscopy (EELS), is compara- 
ble to RIXS and NIXS in the sense that in the inelas- 
tic electron scattering process charge neutral excitations 
are created. The EELS cross section is determined by 
the charge structure fact or Sgjuj) of the system under 
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FIG. 9 Two classes of inelastic energy-loss features in the 
RIXS spectrum of CUB2O4 (iHancock et al.\\2009\. The RIXS 



intensity is represented on a color scale versus transferred en- 
ergy (energy loss) hco and incident energy huj^. The zero-loss 
line is the vertical line at zero transferred energy. X-ray Ra- 
man features are parallel to that line and tend to resonate 
strongly at a specific incident energy, for instance the huj =2 
eV loss feature resonating at ficjk =8992 eV. Emission lines 
appear as diagonal features in this fku - hun^ plot because in 
this case energy of the emitted photon hcou' is roughly con- 
stant. 



fects in the beam. It has the advantage that it is very 
sensitive for lower momentum transfers, around q = 
but its intensity rapidly decreases for large q. Further, 
at large q, multiple scattering effects become increasingly 
important, making it hard to interpret the spectra above 
fl - 0.5 - 1.0 A-i. Since EELS is measured in transmis- 
sion it requires thin samples and measurements in the 
presence of magnetic or electric fields are not possible due 
to their detrimental effect on the electron beam. There 
are no such restrictions for x-ray scattering and EELS 
therefore provides spectroscopic information that is com- 
plementary to IXS/RIXS. Comparing spectra measured 
by these different techniques can be a fruitful method 
both the physical properties of the materials investigated 
and to deepen our understanding of the techniques them- 
selves. 



study ( Schnatterly , 1979 Schuster et al. 2009) and is, 2. Raman and Emission peaks in RIXS spectra 



therefore, closely related to NIXS. When dealing with 
excitations in the valence shell, the RIXS cross section 
can, within certain theoretical approximations, also be 
related to the structure factor, see Section [lV.C[ It should 
be stressed, however, that to what extent, and under 
what circumstances, the RIXS response is proportional 
to S'q(w) is still an open question and an active subject 
of theoretical and experimental investigations. 

The EELS intensity is limited due to space-charge ef- 



Two classes of inelastic energy-loss features are dis- 
tinguished in RIXS, depending on how their energy huj 
changes upon changing the incident energy /icjk. For one 
set of RIXS features, the energy Hoj lost by the photon 
and transferred to the material is independent of the in- 
coming photon energy fiwk and they appear as vertical 
lines in a plot of transferred energy hio versus the incom- 
ing energy ttwu, see Fig. [9] This implies that when the 
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incident photon energy is increased, the energy of the 
outgoing photons /lajk' increases by the same amount. 
This is similar to Stokes phonon losses measured in op- 
tical Raman scattering. Using this analogy one can refer 
to these energy-loss peaks as resonant x-ray Raman fea- 
tures. 

The other class of loss features is characterized by the 
energy of the outgoing photons /iWk' being independent of 
the incident energy fi.Wk- This happens when the absorp- 
tion process promotes a deep core-electron to the vacuum 
continuum far above the Fermi level and subsequently an 
electron from, for example, a shallow core-level decays ra- 
diatively, filling the deep core and emitting a photon. The 
energy of the x-ray photon emitted in such a core-core 
transition is constant, i.e. independent of the incident 
energy Swk. Thus in this situation, when the incident 
energy fiwk increases, the energy lost by the photon hio 
increases by the same amount. As the presence of such 
loss features does not depend on the incident energy be- 
ing close to an absorption edge or a resonance, one can 
refer to such peaks in the RIXS spectrum as x-ray emis- 
sion or x-ray fluorescence features. 

In a plot of the transferred energy huj versus the in- 
coming energy /iWk, x-ray emission features move along 
the diagonal, see Fig. |9] Note that both Raman and 
emission features can appear close in energy. One should 
also consider that for incident photon energies below the 
edge, emission lines occur at constant energy loss. Thus, 
a full description of the RIXS cross section should include 
both phenomena. 



C. Aim of this review 

Here, we review the field of resonant inelastic x-ray 
scattering from elementary excitations in solids, with a 
strong focus on strongly correlated materials. This im- 
plies we do not include RIXS on gas phase atoms or 



molecules, reviewed in Rubensson (2000) and Salek et al. 



(2003), liquids, or non-crystalline and biologically rele- 



vant transition metal complexes ( Glatzel and Bergmann 



2005). Elementary excitations here are understood as 



excitations of electronic, lattice, orbital or magnetic ori- 
gin with an energy of up to a few electron volts above 
the Fermi level. Such low-energy excitations have be- 
come accessible experimentally in the past decade due to 
the progress in instrumentation and x-ray sources. We 
discuss both these experimental developments and the 
simultaneous theoretical developments. Pre-2001 RIXS 



work has been reviewed by Kotani and Shin (2001). 



D. Structure of this Review 

In the following section, we survey the experimental 
methods and foundations of RIXS, including the practi- 
calities of the method, and instrumentation used in soft 



and hard x-ray RIXS. Sections |III| and IV provide a self 



Other reference material includes the books of ISchiilkel 



(2007) and De Groot and Kotani (2008). 



contained review of the theory of RIXS, including effec- 
tive theories and approximate descriptions that are used 
to calculate the complex RIXS process. 

Readers not interested in the experimental methods 
or theoretical background can directly proceed to Sec- 
tion |Vj where the experimental and theoretical investi- 
gations of the past decade on the different elementary 
excitations measured by RIXS are reviewed side by side. 
These include plasmons, inter-band and charge transfer 
excitations, charge density wave, crystal field, orbital ex- 
citations, orbitons and finally magnetic excitations. A 
summary and outlook is presented in Section |VI[ An ap- 
pendix on spherical tensors and a table of symbols are 
included at the end of this review. 



II. EXPERIMENTAL METHODS AND FOUNDATIONS 
A. Practicalities of the method 

1. Choice of Absorption Edge 

The first choice for an experimentalist interested in 
doing a RIXS experiment is which absorption edge to 
work at. The most common edges for transition-metal 
compounds are the transition-metal L- and -RT-edges; less 
work has been done on transition-metal M-edges and lig- 
and edges, such as the oxygen K edge. Although at first, 
one might assume that there should be little difference, 
it turns out that different valence excitations have signif- 
icantly different spectral weights at different absorption 
edges. 

The transition- metal L- {2p 3d) and M- (3p — ?► 3d) 
edges are in many ways comparable. The dipolar transi- 
tions directly involve the 3c? orbitals. This allows access 
to all 3d orbitals providing the experimentalist with the 
opportunity to study dd excitations and orbital effects. 
In addition, since spin is not a good quantum number in 
the intermediate state, one is also able to create single- 
magnon excitations. Since the excited electron effectively 
screens the core hole, there is relatively little shake up 
apart from the excitations made by the dipolar transi- 
tions. 

At the transition-metal K-edge (Is — >■ 3d/ 4p), direct 
transitions into the 3d are quadrupolar and therefore rel- 
atively weak. The major RIXS intensity is related to the 
Is — >■ 4p transitions. Since the excited Ap electron poorly 
screens the core hole, significant shake-up processes oc- 
cur due to the transient presence of a Coulomb potential 
related to the Is hole and 4p electron. These can in- 
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elude transitions aeross the charge-transfer gap. Since 
x-rays in the i^T-edge absorption region carry significant 
momentum, these charge excitations can be studied as a 
function of both energy and momentum. Since no spin- 
orbit couphng is present in the intermediate state, spin is 
a good quantum and one cannot observe single magnon 
excitations, although spin-conserving excitations, such as 
bi-magnons are certainly possible. The ligand edge is 
used less often and provides complementary information 
to L- and M-edge RIXS. 



2. Choice of Incident Energy 

After deciding upon an excitation edge, one has the 
choice of incident energy. In general, in an experimen- 
tal run, scans are made for a range of difi'erent energies. 
The intensities of the features can strongly depend on 
the incoming photon energy. At an L-edge, probing the 
main absorption line usually leads to transitions that are 
metal centered. When using energies in the satellite re- 
gion above the main line, the charge transfer features are 

a few eV 



emphasized. Figure 10 shows this for NiO 
above the absorption maximum, the inelastic processes 
are dominated by CT excitations. At the K edge, one 
can distinguish several regions. In the weak pre-edge re- 
gion, direct excitations are made into the ScZ shell, mak- 
ing it comparable to RIXS at the L- and M-edges. Here, 
dd excitations are emphasized. For the main edge, the 
RIXS spectrum is dominated by charge-transfer excita- 
tions. The most efficient screening mechanisms appear at 
the lowest energies. For example, for cuprates, the most 
efficient screening mechanisms are excitations across the 
charge-transfer gap. Scanning further above the edge 
into regions where the core hole is not well screened, in- 
creases the intensities of the unscreened final states at 
higher energies. 



3. Relative Intensity of Elastic Line 

One of the most striking features of RIXS is that the 
size of the elastic line relative to the inelastic scattering is 
very different at the L- or M-edges . For K- and M- 
edge experiments, the elastic scattering is so large that it 
is typically the tails of this peak that determine whether 
a specific low energy feature can be seen, and not the 
resolution Full Width at Half Maximum (FWHM) - al- 
though a better resolution narrows the elastic peak. In 
contrast, for L-edge experiments, the inelastic scattering 
typically dominates the elastic. It is therefore important 
to understand the elastic scattering, and to the extent 
possible seek to minimize it in order to carry out a suc- 
cessful experiment. 

There are several contributions to the elastic scatter- 
ing. The first of these are phonons. While strictly inelas- 
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FIG. 10 Relative intensity of charge transfer spectral weight 
in NiO (defined as all intensity between 5 eV and 20 eV energy 
loss) compared to the whole RIXS spectrum. The charge 
transfer processes dominate for incident energies more than 
4 eV above the peak in the absorption. Figure taken from 
Ghiringhelli et al. ( |GhiringheUi et al.\ |2005| ) 



tic processes, scattering from phonons is generally unre- 
solved and therefore contributes to the apparent hcu — 
signal. For single crystals of high quality, this may con- 
tribute as much as 50% of the elastic intensity at room 
temperature and may be minimized by cooling to as low 
a temperature as possible. The overall intensity of the 
phonon scattering varies as (q • ^)^, where S is the dis- 
placement vector for the atomic motion, and q is the total 
momentum transfer. As such it increases rapidly at the 
higher momentum transfers where K-edge data tends to 
be taken. In contrast, the smaller q at L- and M-edges 
may be one factor why the elastic response is smaller in 
these experiments. 

Other contributions to the elastic scattering include 
diffuse scattering due to strain or defects in the crys- 
tal. Diffuse scattering is spread throughout reciprocal 
space, peaking in the vicinity of Bragg peaks in charac- 
teristic patterns and is commonly known as Huang dif- 
fuse scattering. This scattering also scales with and 
so becomes an increasing problem at higher momentum 
transfers. The contribution from such scattering is thus 
minimized by choosing crystals of as high a crystalline 
quality as possible and by working at low q. 

A similar source of elastic scattering is the crystalline 
surface quality. A rough surface will produce diffuse, 
static scattering spread throughout reciprocal space. For 
rough surfaces, the integrated intensity of this scattering 
drops off rapidly away from Bragg peaks. Very smooth 
surfaces will give rise to rods of scattering connecting 
the Bragg peaks along directions normal to the surface, 
known as crystal truncation rods. These may be avoided 
by moving away in reciprocal space in a direction parallel 
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to the surface. 

All the contributions listed above are seen in the non- 
resonant scattering channel, and as such have e • e po- 
larization dependence, where e(e ) is the incident (out- 
going) photon polarization. As such, these contributions 
may be minimized by working with a scattering geom- 
etry for which the incident photon is polarized in the 
plane of scattering (tt incident) and the scattering angle 
is 90°. Such scattering geometries are used routinely in 
soft x-ray experiments working at the L- and Af-edges, 
but until recently have been less commonly used for K- 
edges experiments, and this is another reason why the 
K-edge elastic intensity can appear large. Of course such 
requirements are at odds with the desire to move around 
in reciprocal space and investigate the dispersion of a 
given excitation. There are certain cases - for example, 
one or two dimensional materials, where the dispersion in 
one or two directions may be negligible and one may take 
advantage of this fact to maintain the 20 = 90° condition 
whilst still varying the momentum transfer in a relevant 
direction. This has proved of great utility in studying 
such low dimensional materials at both L- and X-edges. 

Finally, when considering the relative strength of the 
elastic and inelastic contributions, one must consider the 
size of the resonant enhancement for the inelastic chan- 
nel. For the L-edge experiments, it is typically, though 
not always, a direct RIXS process, and as such the in- 
elastic scattering is relatively strong. Further, the res- 
onances are large. This latter fact is evidenced by the 
strong "white line" features seen at L-edges in absorp- 
tion data and is contrast with the very weak peaks, or 
absence of peaks in M- and K-edge XAS spectra. As a 
result the resonant enhancements are significantly larger 
for the L-edge experiments. The net result of these two 
effects is that the inelastic scattering is much stronger at 
an L-edge than a M- or K-edge. 



4. Momentum Dependence 

Generally, RIXS spectra on elementary electronic ex- 
citations show a dependence on the momentum that the 
photon transfers to the solid. This experimental observa- 
tion concurs with the theoretical observations of Platz- 



that the photon momentum is within the first Brillouin 



zone). In studies of La2Cu04, Kim et al. (20071 showed 



man and Isaacs (Platzman and Isaacs 1998) who argued 



that while the matrix elements involved in the resonant 
scattering process are all highly local, the coherent sum 
of such transitions on different sites would lead to over- 
all crystal momentum conservation, and further, that 
momentum-breaking processes such as the creation of 
phonons would be small (by a factor wd/F, where ujd 
is the Debye frequency and F the core-hole lifetime). 

In exploring the momentum dependence of the inelastic 
scattering in a K-edge RIXS experiment, one has a choice 
of which Brillouin zone in which to take the data (for TM 
L-edge and Af-edge RIXS experiments, the energy is such 



that to a first approximation the spectra did not depend 
on which Brillouin zone they were measured in. They 
therefore concluded that the momentum dependence ob- 
served in RIXS experiments arises from the crystal lat- 
tice, i.e from electronic states obeying Bloch's theorem 
and respecting the periodicity q = k -I- nG, where q is 
the total momentum, k the reduced wave- vector and G a 
reciprocal lattice vectoi[^ The reason that the observed 
RIXS spectra depend on k but not on nG is that the 
local core-hole interaction is the main excitation mech- 
anism, and this has the translational symmetry of the 
lattice. This point will be discussed in detail in Section 

A detailed study of the Brillouin Zone dependence of 
the zone center scattering in Nd2Cu04 was carried out by 



Chabot-Couture et al. (20101. After careful normaliza- 



tion of their data, the authors conclude that there was a 
small Brillouin zone dependent component of about 10% 
of the inelastic signal. Its intensity is best described by 
an empirical form; cos'^{9 — 77), where d is half the scat- 
tering angle and r\ the deviation of the scattering vector 
from the surface normal (in this case along the c-axis 



of the crystal). Chabot-Couture et al. (2010) argue that 



the functional form originates from a photon polarization 
dependence of the cross-section - not the momentum- 
dependence of the excitation spectrum of the solid. This 
implies it is a correction to the 4p o,s a spectator approx- 
imation to K-edge RIXS discussed in section [TV.CI 



5. Normalization of data 

During the RIXS scattering process, the incident and 
outgoing photons are attenuated by Auger emission and 
core-level fluorescence as they propagate through the 
sample. This attenuation, called self absorption, de- 
pends strongly on the scattering geometry used, and de- 
termines the illuminated sample volume. Correcting for 
this extrinsic effect is thus important in order to com- 
pare spectra collected in different scattering geometries. 



In Ref. (Chabot-Couture et al. 



2010[ ), the authors de- 
scribe how an extrinsic 20-25 % attenuation (or enhance- 
ment) of the measured Cu K-edge RIXS spectra can be 
corrected for by using the fluorescence signal of a known 
emission line (measured in the same scattering geome- 
try) and Cu K-edge X-ray absorption spectra. Similar 
considerations also apply for L- and Af-edge experiments 
and in general care should be taken to properly account 
for self absorption. 



^ Note that contrary to the present convention many papers use 
Q and q as the total and reduced momentum, respectively. 
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FIG. 11 Comparison of RIXS spectra on La2Cu04 taken at 
different total q positions but with an equivalent reduced wave 
vector K for three different q's corresponding to k = (0,0), 
(7r,0) and (tt, tt), from top to bottom. The two-dimensional 
reciprocal space net is shown in the top panel. The same 
symbols are used both to denote the q position in reciprocal 
space and to show the spectra. All data were obtained with 
hLUk=8991 eV (iKim et aZ.l|2007l. 
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FIG. 12 (Color online) The schematic layout of the SAXES 
spectrometer, at the Swiss Light Source. It is based on a 
variable-line-spacing spherical grating and a two dimensional 
position sensitive CCD detector. Panel A: the main optical 
parameter symbols are defined. Panel B: the main mechan- 
ical components are sketched together with their degrees of 
freedom. Figure from Ghiringhelli et al. ( [Ghiringhelli et aL\ 



B. Instrumentation 



2006b I 



A RIXS instrument may be divided into two parts, the 
beamline, which is required to deliver a well-collimated, 
highly monochromatic, and often focused, beam onto the 
sample, and the spectrometer, which must provide the 
sample environment, manipulate the scattering angles to 
vary the momentum transfer, collect the scattered radi- 
ation over an appropriate solid angle, and energy an- 
alyze the scattered radiation. All of these are highly 
non-trivial tasks and it has taken a great deal of work 
performed over a large number of years by a large num- 
ber of people to reach the current state-of-the-art. In 
parallel to the instrumentation development there has 
been a tremendous development of synchrotron sources 
in terms of ever improving brightness in the relevant en- 
ergy ranges. The combined effect is that there has been 
dramatic improvement in energy resolution that can be 
achieved in momentum-resolved RIXS instruments and 
still have sufficient count rates for experiments to be per- 
formed in a finite amount of time. 

Hard x-ray instruments are invariably built around 
Bragg crystal optics to achieve the very high energy res- 
olution required, ultimately benefitting from the very 
highly ordered crystalline lattice of (typically) Si or Ge 
that is well matched to the x-ray wavelengths being uti- 
lized. Such schemes do not work well in the soft x-ray 



regime, where the long wavelengths prohibit the use of 
crystal optics. Instead artificial periodic structures are 
utilized - diffraction gratings. This difference means a 
hard x-ray RIXS spectrometer is fundamentally differ- 
ent from a soft RIXS instrument. In the following two 
sections, an example of each is presented. 



1. A soft RIXS instrument 

As mentioned above, diffraction gratings lie at the 
heart of soft x-ray beamlines and spectrometers. These 
are lined gratings that are typically operated in graz- 
ing incidence. The lines act as a series of slits and the 
light diffracts from these. The resolution is determined 
by the number of lines per mm. These are not par- 
ticularly efficient optics, and offer small collecting an- 
gles when used as analyzers. Therefore bright x-ray 
sources, synchrotrons and in the coming years soft x-ray 
free-electron lasers, are a pre-requisite for high-resolution 
RIXS. Pioneers in the construction of such instruments 



were 


Callcott et al. 


(1986 


) and 


Nordgren et al. ( 


1989 


) and 


others (Chuanga et al. 2005 Cocco et al. 2001 


2004 


Dallera et al. 1996| Hague et al.\ 


2005 Shin et al. 


1995 
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FIG. 13 The mechanical design of SAXES. The sample is 
placed at the far right. Diffraction grating is in the vacuum 
chamber approximately Im from the sample. The CCD is 
inclined at an angle of 20° and is 5.1 m from the sample. 
The whole spectrometer is supported by a single steel girder 
and can rotate around a vertical axis about the sample on air 
cushions. This allows the momentum transfer to be varied. 



Figure from Ghiringhelli et al. ( 2006b I 



Tokushima et al. 2006 ) 



As an example, we discuss here the instrument at the 
Swiss Light Source. It consists of a beamline, known as 
ADRESS for "ADvanced RESonant Spectrscopies" , and 
a spectrometer, known as SAXES for "superadvanced x- 
ray emission spectrometer" . The beamline is based on 
a design that utilizes a plane grating monochromator 
(PGM) diffraction grating with a line spacing of up to 
4200 mm^^, together with coUimating mirrors upstream 
of the PGM and focusing mirrors downstream of the 
PGM. (Another common choice is to combine these latter 
two optical elements into a single element - a spherical 
grating). Operating in the energy range 400-1800 eV, 
ADRESS has a number of key features that make it well 
optimized for RIXS experiments. In particular, it deliv- 
ers a high flux: - lO^^ ph/s/0.01% BW at 1 keV, with 
a variable incident polarization in a small spot (8.5 /im 
X 100 fim). This latter is essential in allowing high en- 
ergy resolution in the spectrometer, as we shall see below. 
Further, the energy resolution is extremely good, with a 
resolving power of over 30,000 at 1 keV (i.e. an energy 
resolution of 30 meV) fPSH |20T0 1 . 

The spectrometer, SAXES, must collect the scattered 
radiation from the sample and provide both momentum 
and energy resolution. Its design and performance is de- 
scribed in detail in 'Ghiringhelli et al. (2006b I. The basic 
design is shown in figure [T2[ 

The spectrometer is designed around a single optical 
element, a spherical variable line space (VLS) grating 
which collects the scattered radiation and disperses it 
onto a CCD camera. The choice of a single grating - 
as opposed to a plane grating and a concave collection 
mirror - simplifies the minimization of higher order op- 
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FIG. 14 The calculated resolving power and energy resolution 
of SAXES. Top panel: The expected energy resolution for 
various openings of the entrance slit. Bottom panel: The 
corresponding values of the resolving power, E/AE. Figure 
from iGhiringhelh et al.\ (|2006b|) 



tical aberrations at the price of reduced collection angle. 
The overall energy resolution of the instrument is largely 
determined by the line spacing of the grating and the 
spatial resolution of the CCD camera. In recent years 
there has been steady improvement in the manufactur- 
ing technology for making these gratings. The grating in 
SAXES has an average line spacing ao=3200 mm~^ and 
a fixed radius of curvature oi R = 58.55 m. The grating 
angle and position, a and ri respectively, are optimized 
at each energy to minimize coma aberration. Slits 5 mm 
downstream of the sample determine the illumination of 
the grating by the scattered beam and may be aperturcd 
down to improve the resolution. The grating is estimated 
to be 3-4% efficient, operating at an incident angle of 
88.3°. 

The dispersed light is then allowed to propagate in vac- 
uum until it is incident on the CCD camera. In principle, 
the energy resolution is then determined by the distance 



the light is propagated (r2 = 4.03 m, figure 12 ) and the 
pixel size on the camera. In fact, the physics of the CCD 
means that while the actual pixel size is 13.5 /zm x 13.5 
/im, a single photon will appear to be observed by sev- 
eral pixels due to charge leakage between neighboring 
pixels (Ghiringhelli et al. 2006b I. The camera is there- 



fore operated at a grazing incidence (20°) to match the 
effective pixel size to the charge distribution caused by a 
single photon. Finally, it is cooled by liquid nitrogen to 
— 100 °C to reduce the dark current to negligible levels. 

The whole spectrometer (grating plus CCD) is then 
mounted on a large (5 m) girder as shown in figure 13 
This allows control of the vertical and transverse posi- 
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FIG. 15 Schematic of optical components in the 30 IDB in- 
elastic x-ray scattering beamline at the Advanced Photon 
Source, USA. Figure courtesy of D. Casa (APS). 



tion of the spectrometer together with its roll and yaw, 
to correctly position the spectrometer about the scat- 
tered beam. Finally, the whole girder is able to rotate in 
the horizontal plane about the sample position - moving 
on air cushions over a granite substrate - to control the 
scattered momentum transfer. 

The calculated performance of this instrument is 
shown in figure [14] as a function of photon energy and 
for various openings of the entrance slit (i.e. the slit just 
downstream of the sample). At the smallest slit settings 
(5 /im), it is seen that a resolving power of better than 
10,000 is expected throughout the whole energy range. 
One can see that the effective pixel size from the detector 
is the largest single contribution to the energy resolution 
suggesting that further improvements in such spectrome- 
ters will require associated advances in detector technol- 
ogy. Overall, this instrument is, when coupled with the 
ADRESS beamline, capable of sub-lOQ meV resolutions. 
As shall be seen in later sections, this has led to dramatic 
advances in the experimental results. 



2. A hard RIXS instrument 

RIXS instruments in the hard x-ray regime are built 
around Bragg crystal optics and as such, are fundamen- 
tally different from soft RIXS instruments, such as the 
one discussed above. Here, we describe the so-called 
MERIX instrument (for Medium Energy Resolution In- 
elastic Xray) which operates on the 30 ID beamline at 
the Advanced Photon Source, USA. Other, related in- 
struments include an earlier version of MERIX designed 
by the same team and in operation at 9ID, APS (Hill 



at medium energy resolutions (~ 70 meV) for incident 
photon energies from 4.9 keV to 15 keV. It uses two lin- 
early polarized 2.5 m long 30 mm period devices as a ra- 
diation source. The primary monochromator comprises 
two water-cooled diamond (111) crystals in a fixed exit 
geometry. This is capable of handling the heat load from 
the undulators and provides an incident band pass of 
^0.6 eV at 10 keV. This beam is then incident on a sec- 
ondary monochromator comprised of 4 asymmetrically 
cut Si(400) crystals. These are arranged in two pairs to 
provide a fixed output beam with zero offset from the 
low resolution beam. This arrangement allows the ex- 
perimenter to quickly switch from the high-band pass, 
high flux beam from the diamond monochromator to the 
narrow-band pass, lower flux configuration, a feature that 
is frequently of great value. The design of the Si(400) 
high-resolution monochromator is such that it provides a 
more or less constant band pass of ~70 meV over the 5-15 
keV range of the beamline, with an average spectral effi- 
ciency (the number of incident photons that lie within the 
required band pass that are passed by the optic) of 70 %. 
Following the high-resolution monochromator, there are 
two plane mirrors arranged in a Kirkpatrick-Baez con- 
figuration to provide horizontal and vertical focusing of 
the monochromatic beam. The focused spot size is 45 
x6/im^ (H X V) and the incident flux on the sample is 
2 X 10^^ photons/s/lOOmeV at 9 keV. A schematic of the 



beamline layout is shown in flgure 15 



et al. 2007) 



The beamline is designed to provide high incident flux 



The MERIX spectrometer itself is shown in flgure 16 
It features a 2 m detector arm which is able to move 
in both the vertical and horizontal plane. This allows 
the experimenter to vary the incident polarization with 
respect to the momentum transfer (cr incident or tt in- 
cident geometries). Further, fixing the scattering angle 
at 90° in the horizontal plane with horizontally polar- 
ized incident light minimizes the intensity of the elastic 
scattering due to the e • e dependence of most sources 
of such scattering. The analyzer is situated inside a tank 
filled with He gas. This minimizes absorption and back- 
ground scattering due to windows that would otherwise 
have to be present. The analyzer may be moved along a 
linear track to ensure it is placed on the Rowland circle. 
Common distances used are 1 and 2 m. 

The analyzer is the heart of the instrument. In hard 
RIXS instruments, these are again based on Bragg op- 
tics. A single crystal is positioned to intercept the scat- 
tered radiation, and Bragg reflect it onto a detector. The 
scattered energy is then scanned by changing the Bragg 
angle of the crystal. In order to maximize the scattered 
intensity (at the expense of momentum space resolution) 
a spherically-bent analyzer is used to intercept a large 
solid angle. At MERIX, analyzers are typically 10 cm in 
diameter and intercept ~10 mSr at 1 m. These are oper- 
ated in the so-called Rowland circle geometry, such that 
the sample, analyzer and detector are all mounted on a 
circle of diameter R. When this radius is equal to radius 
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been used (following their highly successful use for very 




FIG. 16 The MERIX spectrometer for hard x-ray inelastic 
x-ray scattering at 30 IDB, Advanced Photon Source, USA. 
The sample is located at the center of the instrument and 
the beam is incident from the right hand side of the figure. 
The red cones show the path of the scattered light onto the 
analyzer and back onto the detector. There are three degrees 
of freedom on the sample and the detector arm may move 
in both the vertical and horizontal plane. The detector is 
typically a position sensitive detector. Figure courtesy of S. 
Coburn (BNL). 



of curvature of the analyzer, one has the so-called Johan 
geometry, which is an approximation to the ideal figure 



(Johansson geometry). This is shown in figure 15 1. For 



the very center of the analyzer, positioned right on the 
Rowland circle, all appropriate scattering conditions are 
met. For the parts of the analyzer away from this point, 
aberrations are introduced. 

The energy resolution of such analyzers has four con- 
tributions; i) the intrinsic Darwin width of the Bragg 
reflection being utilized - this is a function of the Debye- 
Waller factor, the structure factor and the Bragg angle 
itself - and the strain introduced by the bending which 
broadens this resolution according to the equations of 
dynamical diffraction theory, ii) The contribution in the 
vertical and horizontal directions, due to the Johan er- 
ror discussed above. This increases quadratically with 
r/i?, where r is the radius of the analyzer, and decreases 
with increasing Bragg angle, iii) The source size result- 
ing from the finite beam footprint on the sample. This 
increases the energy resolution, linearly increasing with 
the angular width of the source as seen by a point on the 
analyzer and decreasing with increasing Bragg angle. 

The introduction of strain due to the bending of the 
analyzer causes the largest contribution to the energy res- 
olution and therefore in recent years, diced crystals have 



high-resolution, meV non-resonant spectroscopy (Ver- 
beni et aL\ [2005 )). The dicing however, leads to the 



fourth contribution to the energy resolution; iv) The fi- 
nite block size, due to the angular width of that block as 
seen by the source point. This is linear in block size and 
decreases with increasing Bragg angle. Because the fi- 
nite size contributions, and the Johan corrections all scale 
with cot 6b, these are minimized when at back-scattering, 
Og = 90°, thus diced analyzers work extremely well in 
these conditions. 

Unfortunately, in a RIXS experiment, and in contrast 
to non-resonant IXS, the incident energy is chosen by 
the resonance process. One must then find a set of Bragg 
planes in a suitable analyzer material for which the Bragg 
angle at this energy is close to = 90°. This is by 
no means a trivial task. It turns out that at the Cu K- 
edge, Ge(733) has 9b — 86.8° and this is sufficiently close 
to 90° that relatively good energy resolutions have been 
achieved. However, for other transition metal i^-edges, 
neither Si nor Ge has any appropriate reflections and 
resolutions are compromised. Further, even at the Cu 
K-edge, dicing still limits the ultimate resolution achiev- 
able. Therefore in the very latest generation instruments, 
the traditional point detector is replaced with a position 
sensitive detector (jHuotari et al. 



2006 2005) 



As shown in figure 15 a), and as first emphasized 



by Huotari et al. (2006 2005), the scattered radiation 



is dispersed from the Bragg reflection from a single block 
with an energy-angle correlation given by Bragg's law. 
With a point detector, this spread is integrated over and 
gives rise to the broadened energy resolution (point iv) 
above). However, with a position senstive detector, this 
information can be recovered, removing this broaden- 
ing effect. Further, there is the added advantage that 
these energies - corresponding to different energy trans- 
fers, may be counted simultaneously, greatly increasing 
the throughput of the instrument. The block size no 
longer contributes to the energy resolution, but rather 
now determines the range of energy transfers that may be 
measured simultaneously. Such dispersive schemes, echo- 
ing those of the soft RIXS instrument discussed above, 
are the latest instrumentation development, improving 
both resolution and count rate (see figure|3|. With such a 
scheme, MERIX has achieved 83 meV FWHM at the Cu 
K-edge (overall resolution), of which the biggest contri- 
bution is the monochromator. Further, good resolutions 
(< 100 meV) may be achieved at the other transition 
metal iiT-edges using this scheme and Si or Ge reflec- 
tions. In addition, groups are now exploring alternate 
analyzer materials, including sapphire (AI2O3), quartz 
(Si02) and LiNbOa, which offer many more reflections 
and the chance to move closer to back scattering and 
thus improved resolution. While a great deal of devel- 
opment will be required for each of these, it seems likely 
that the resolution of hard RIXS will continue to improve 
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in the coming years. 



III. THEORY OF RIXS: I. FUNDAMENTALS 

In RIXS, the soUd is taken from a ground state with 
energy to a final- state with excitations and an energy 
Ef. The energy and momentum of the excitation is de- 
termined by the difference in photon energy hoji^ — hui^i 
and momentum Kk' — Kk, respectively. The RIXS inten- 
sity can in general be written in terms of a scattering 
amplitude as 

/(w,k,k',e,e') = ^| J"/g(k, k', e, e', Wk)! 
/ 

X S{Ef + hioi^i — Eg — huk) 

where the delta function enforces energy conservation 
and the amplitude J-"/g(k, k', e, e', Wk) reflects which ex- 
citations are probed and how, for instance, the spectral 
weights of final state excitations depend on the polar- 
ization vectors, e and e' of the incoming and outgoing 
x-rays, respectively. The following sections derive the 
RIXS scattering amplitude and demonstrate how it can 
be broken down into separate pieces. 



First, we need to derive a general expression for the 
RIXS scattering amplitude. Section III. A looks at the 
interaction between photons and matter. RIXS refers 
to the process where the material first absorbs a pho- 
ton. The systems is then in a short-lived intermediate 
state, from which it relaxes radiatively. In an experi- 
ment, one studies the x-rays emitted in this decay pro- 
cess. This two-step process cannot be described simply 
by using Fermi's Golden Rule, but requires a higher- 
order treatment, known as the Kramers-Heisenberg equa- 



scattering amplitude is a linear combination of these fun- 
damental scattering amplitudes weighted by the angular 
functions T^(k, k', e, e'). 

The next step is to understand the fundamental scat- 



tion (Kramers and Heisenberg 1925). Since the absorp- 



tion and emission are single-photon processes, the inter- 
actions between the x-rays and the material are domi- 
nated by the terms in the cross-section proportional to 
p • A, where p is the momentum of the electrons in the 
material and A is the vector potential of the photon. 
The interaction between the x-rays and the material de- 
pends on external quantities, such as wavevector k and 
polarization vectors e of the x-rays, and operators, such 
as p and r. As a result the electronic transitions are 
intermingled. Sections III. A. 2 and III. A. 3 describe how 
the scattering amplitude can be split into an angular and 
polarization dependence T'^(k, k', e, e') related to the ex- 
perimental geometry and spectral functions F^(fc, k', Wk) 
that measure the properties of the material, see Table 
|Tj This separation can be done exactly. It is important 
to note that there are only a finite number of fundamen- 
tal scattering amplitudes F^{k, fc', Wk) and that the RIXS 



tering amplitudes. This is done in Section IV This can 
be done numerically, as is outlined in Section IV. A In ad- 
dition, several authors have used approximation schemes 
in order to provide more insight into the scattering am- 
plitude. Generally, the approximations involve the prop- 
agation of the system in the time between the absorption 
and emission processes. The schemes generally allow the 
separation of the fundamental scattering amplitudes into 
a resonance function P(a;k,i^k') and an effective transi- 
tion between ground and final states {f\W^\g), see Ta- 
ble |l] The resonance function gives the strength of the 
fundamental scattering amplitude which is a combina- 
tion of radial matrix elements of the transition operators 
and energy denominators that describe the resonant ef- 
fect as a function of Wk. The effective transition opera- 
tors create excitations in the valence shell similar to an 
optical excitation. In certain cases, these operators can 
also be related to correlation functions such as the dy- 
namic structure factor. The approximations depend on 
the RIXS process. Direct RIXS is approximated by using 
a fast-collision approximation, see Section [IV. B[ Indirect 
RIXS can be approached via perturbative methods or 
an ultra-short core-hole lifetime expansion, see Section 

HvTi] 



A. Interaction of light and matter 

To develop the theory of RIXS, we first need to de- 
rive the Hamiltonian that describes the interaction of 
the incident X-ray beam with the electrons in the sam- 
ple. The interaction terms in this Hamiltonian are small, 
controlled by the dimensionless fine structure constant 
a — e^/4:TT€oHc « 1/137, with e — \e\ the magnitude 
of the elementary charge and eq the permittivity of free 
space. Therefore they can be treated as a perturbation 
to the terms in the Hamiltonian that describe the system 
under study. To second order in such a perturbation the- 
ory, we obtain the Kramers-Heisenberg formula, which 
describes RIXS very well. We need to go to second or- 
der because two interactions are needed: one to create 
the core hole, and one for the subsequent radiative de- 
excitation. 



1. Kramers-Heisenberg cross section 

The incident X-rays are described by an electromag- 
netic field with vector potential A(r,t). The coupling 
between such a field and electrons is given by the the- 
ory of quantum electrodynamics. It is common to start 
from the exactly solvable case of a single electron without 
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TABLE I In the theory of RIXS, the scattering amphtude Tjg occurring in the Kramers- Heisenberg equation is separated 
into several pieces. One can split off the angular and polarization dependence T"^, leaving fundamental scattering amplitudes 
F^. Several approximation schemes then break down these scattering amplitudes into a resonance function P and effective 
transition operators . 



Factorizing the Kramers-Heisenberg Scattering Amplitude Tjg 



scattering amplitude 
exact J/g(k, k', e, e', cjk) 



fundamental scattering amplitudes x angular and polarization dependence 
F-(fc,fc',a;0 ■ T-(k, k', e, e') 



approximate 



fundamental 
scattering amplitude 



effective 

resonance function x transition operators 



potentials (A and electric potential (/)(r,t)). Then, the 
potentials are (perturbatively) introduced and one takes 
two limits. The first of these is that the electrons travel 
at speeds, v, small compared to the speed of light. This 
is a good approximation even for, for instance, copper Is 
core electrons, where we estimate v ~ hZ/mao « 0.21c 
with Z the atomic number for copper and ao the Bohr 
radius. At first glance, v/c might not appear small in this 
latter case, but 7 = l/-y/l — v'^/c^ « 1.02 and relativistic 
effects are still small. The second limit is that the poten- 
tials related to both the electrons and photons in the sys- 
tem are small compared to twice the mass of the electron: 
e0/2mc^, e |A| /2toc ^ 1 (m is the electron mass). Al- 
though the intrinsic potentials of materials diverge close 
to the nuclei, they may be treated consistently within 



the whole procedure for Z <C 137 (see Messiah (1962), 
page 948). Photon potentials at existing x-ray sources 
satisfy these limits. However, in future very strongly fo- 
cussed X-ray Free Electron Lasers, the electric field of 
the photon is projected to exceed 10^^ V/m, which gives 
e |A| ^ 2mc at a photon energy of ~ 8 keV and these 
approximations are no longer valid. However, for the 
purposes of this review, such effects are neglected and 
the formalism is developed for non-relativistic electrons 
in small potentials. 

In these limits, one obtains for a system with N elec- 



trons, in SI u nits (see pages 944-947 in [Messiah] ([1962| , 
pages 85-88 in|Sakurai| ((l967[)) 



N r 



(pi +eA{Yi)) 
2m 



2m 



■B(r,)- 



eh 



2(2r 



(Ti ■ E(r,) X (p, + eA{vi)) ~ (p, + eA(rj)) x E(r,) 



+ 



8(TOc)^eo 



Coulomb 



^f^^H. ( oLa^E + ^ ) (1) 



where Pi, and cr,; are, respectively, the momentum 
and position operators and the Pauli matrices acting on 
electron i. A(r) is the vector potential, E(r) = — Vi/) — 
dA/dt, the electric field, and B(r) = Vx A, the magnetic 
field, a^l annihilates (creates) a photon in the mode with 



wave vector k and polarization vector e. The second 
term yields the Zeeman splitting, and the third includes 
spin-orbit coupling. The interaction of electrons with 
an external electric potential and with other electrons 
and nuclei in the sample (including the Darwin term) 
are all described by ffcouiomb- The vector potential can 
be expanded in plane waves as 



Mr) - E 



), (2) 



where V is the volume of the system. 

In order to derive the photon scattering cross section 
one splits the Hamiltonian H into an electron-photon in- 
teraction part, H' , and the remaining terms, Hq, which 
describe the electron and photon dynamics in the ab- 
sence of electron-photon interactions. H' is then treated 
as a perturbation to Hq. To calculate the RIXS cross 
section in this perturbation scheme, it is assumed that 
there is a single photon in the initial state with momen- 
tum hk., energy /iwk and polarization e that is scattered 
to {fik' , huk' , e') in the final state. Photon scattering then 
induces a change in the material from ground state \g) to 
final state |/), with energies Eg and Ef respectively. In 
the process, the photon loses momentum hq = Kk ~ Kk' 
and energy hiu = fi-ajk — huju_' to the sample. Fermi's 
Golden Rule to second order gives the transition rate for 
this process: 

27r ■ 



E 

f 



[f\H' 



S) 



E 



{f\H'\n) {n\H'\g) 



E„ — En 



(3) 



where the initial state 
In) and the final state 



If) = 



l5;ke) 
l/;k'e 



intermediate state 
I are eigenstates of 



Hq with energies Eg = Eg + huj^, £"„ and Ef = Ef + huji^^i, 
respectively. The first order amplitude in general dom- 
inates the second order, but when the incoming X-rays 
are in resonance with a specific transition in the material 
{Eg w En), then the second order terms become large. 
The second order amplitude causes resonant scattering, 
while the first order yields non-resonant scattering. 
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In order to derive H' it is useful to classify the terms of 
Eq. ([I]) by powers of A. Terms of H that are quadratic in 
A are the only ones to contribute to the first order ampli- 
tude, because they contain terms proportional to a\^i^,a-^^^ 
and aj^^aj^/g,. To be specific, the quadratic contribution 
from the first term of H gives rise to non-resonant scat- 
tering, while the third term of H yields magnetic non- 
resonant scattering. Although both appear in the first 
order scattering amplitude, they in principle also con- 
tribute to the second order, but we neglect these pro- 
cesses because they are of order a^/'^. 

The interaction terms linear in A do not contribute to 
the first order amplitude, but do contribute to the second 
order. They thus give rise to resonant processes. In the 
following, we neglect such contributions that come from 
the third term of Eq. ([I]), because they are of second 
order in two separate expansions. Firstly, this term of 
H is of second order in the limits discussed above, and 
secondly, it appears in the second order of the scattering 
amplitude. 

Finally, all terms in Eq. ([T]) that are independent of A 
are included in i?o- The relevant remaining terms are 



N r 



(2mc)2 



m 2m 



eh 
2rn 



• X A r, 

at 



(Ti-V X A{Yi) 



(4) 



where the gauge was fixed by choosing V • A(r) = so 
that A • p = p • A. 

The two terms of H' that contribute to the first order 
amplitude are the one proportional to A^ and the cr ■ 
{dA/dt) X A term. The latter is smaller than the former 
by a factor tujj]^(^i)/mc? <C 1, and is therefore neglected. 
The first order term in Eq. ([s]) then becomes 



e 

2^ 



;f|^A2(r,)|g) 



he" 



2mVeo ^WkWk 



= (/iE' 



Iff) 



(5) 

When the incident energy fiwk is much larger than any 
resonance of the material, the scattering amplitude is 
dominated by this channel, which is called Thompson 
scattering. In scattering from a crystal at zero energy 
transfer, this term contributes amongst others to the 
Bragg peaks. It also gives rise to non-resonant inelastic 
scattering. In practice, RIXS spectra show a strong res- 
onance behavior, demonstrating that for these processes, 
it is the second order term that dominates the first or- 
der term. We therefore omit the A^ contribution in the 
following. More details on non-resonant inelastic X-ray 
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when huj]^ matches a resonance energy of the system, 
and the incoming photon is absorbed first in the inter- 
mediate state, creating a core hole. The denominator 



Eg -t- ?icjk — En is then small, greatly enhancing the sec- 
ond order scattering amplitude. We neglect the other, 
off-resonant processes here, though they do give an im- 



portant contribution to non-resonant scattering ( Blume 



1985). The resonant part of the second order amplitude 



IS 



N 



EE 



2m2Veoi/aJki^k' 

(/I e-^i^'"-' {e'* • p, - f • k' X e'*) |n) 
Eg + hujk - En + iTn 

X {n\ e'^ "^ (e ■ Pj + ^cTj • k x e j |g) , 



(6) 



where a lifetime broadening r„ is introduced for the 
intermediate states. This accounts for the many non- 
radiative interaction terms that are not included in H' 
(for example Auger decay) , which make the intermediate 
states very short lived. 

Resonant scattering can thus occur via a magnetic and 
a non-magnetic term. An estimate shows that the latter 
dominates. The size of localized Is copper core orbitals 
is roughly a^jZ k, 0.018 A so that for 10 keV photons 
the exponential e*'' '' is close to unity and can be ex- 
panded. The non-magnetic term can induce a dipole 
transition of order |p| ~ hZja^ ^ 5.9 • 10"^'^ kg m/s, 
whereas the magnetic term gives a dipole transition of 
order (k • r);i|k| /2 - 2.5 • lO^^s kg m/s. We thus ig- 
nore the magnetic term here, and the relevant transition 
operator for the RIXS cross section is 



(7) 



where a prefactor has been introduced for convenience in 
the following expressions. 

The double-differential cross section /(w, k, k', e, e') is 
now obtained by multiplying by the density of photon 
states in the sohd angle (Kl (= V k'^rf|k'|(jl7/(27r)^), and 
dividing by the incident flux c/V ( [Blume 1985| Sakurai 
19671 ISchiilkel [20071): 




/(a;,k,k',e,e') 



2 2 



k-Wk^ |J"/g(k,k',e, e',a;k,a;k' 



X 8(Eg - Ef + huj), 



(8) 



1 

where the classical electron radius = -t^ — The 
scattering amplitude at zero temperature is given by 



7"/g(k,k',e,e',Wk,Wk') = 



iflV'^H {n\V\g) 
E„ + fiWk - En+ iTn 



(9) 



where the prime in V indicates it refers to transitions 
related to the outgoing X-rays. Eqs. ([S]) and ([£]) are 
referred to as the Kramers-Heisenberg equations, which 
are generally used to calculate the RIXS cross section. 
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Alternatively, we can rewrite the denominator for the 
intermediate-states in terms of a Green's function, which 
is also referred to as the intermediate-state propagator, 
which describes the system in the presence of a core hole: 



G(^k) 



1 



\n){n\ 



Z\<i- H ^ Z-k~ En 

where \n) forms a complete basis set and 



Zk = E„+ huJk 



(10) 



(11) 



where T is taken to be independent of the intermediate 
states. The quantity Zk is the energy of the initial state 
combined with the finite lifetime of the core hole. In the 
following we will often suppress the explicit label k of 
Zk and denote it simply by z, with an implicit incident 
energy dependence. With the core-hole propagator G 
and transition operators D in place, the RIXS scattering 
amplitude Tfg finally reduces to the elegant expression 



fa 



{f\V'^Giz^)V\g). 



(12) 



2. Scattering amplitude in dipole approximation 



In the previous section, Eqs. ^ and ([£]) give the 
Kramers- Heisenberg expression for RIXS. The next step 
is to separate the part pertaining to the geometry of the 
experiment from the fundamental scattering amplitudes 
that relate to the physical properties of the system, see 
Table [l] In addition, better-defined transition operators 
will be obtained. Due to the complexity of the multipole 
expansion, we first give a derivation in the dipole limit 
allowing the reader to better follow the arguments. In 
the next section, we present the higher order transitions. 

In the dipole limit, one asumcs that e**''"' = gik-Ri 
where Ri indicates the position of the ion to which the 
electron i is bound. Note that is not an operator. 
This has as a result that the electronic transitions are 
due to the momentum operator p and Eq. [7] becomes 



V = e T> with D = 



N 

\ ^ „ik R, 



Pi, 



(13) 



Generally, the matrix elements are expressed in terms of 
the position operator r. For example, in the absorption 
step, one can write 



zk-R, 2 



N 



oik-Ri 



J2j—iEn-Eg){n\r, \g) 

N 

J2e^^-^^ {n\r,\g} 



where huju = En — Eg. The operator thus reduces to 
the well-known dipole operator D = J2i=i^^^ ^'^i H^^l 
causes electronic transitions. 

The next step is to separate the part that pertains to 
the geometry of the experiment (the polarization vectors 
e' and e) from the physical properties of the system. Ul- 
timately, our interest lies in the spectral functions of a 
material. The geometry is chosen in an optimal way to 
measure them. Using spherical-tensor algebra, see Eq. 



( A3 ) in the Appendix, we can rewrite the scattering am- 



plitude, Eq. (12), remaining in the dipole limit, Eq. (13 1, 

as 



^f9 = Et^l'^ii^ 



(/|[e'*,6]-.[Dt,G(zk)D]n5) 



using the shorthand [h ■ ■ ■ In] = (2^i + l) • • • (2^„ + l); 
is a normalization constant, and [ , ]^ is a tensor product, 
see Appendix. Since the tensor product couples tensors 
of rank 1 (the polarization vectors and the position vector 
r), the rank x of the tensor products can assume the val- 
ues 0, 1, and 2. The fundamental scattering amplitudes 
are given by 



F-(zk) = (/|[Dt,G(zk)Df 1.9) 



(14) 



For each value of x, there are 2x + 1 components 
with q = —X, —X -|- 1, • • • , X. Note that, whereas there is 
an infinite number of different scattering amplitudes, for 
dipole transitions, there are only nine fundamental ones 
(3x3 = 1-1-3-1-5 = 9) . All the other possible scattering 
amplitudes are combinations of these fundamental scat- 
tering amplitudes with a weighting determined by the 
angular dependence; 



T-(e,e') = Nr 



(15) 



which again has nine components ^^(e, e'). For x = 0, 1, 
the angular dependence is given by the inner product, 
r^(e, e') = le'* • e, and the outer product, T^(e,e') = 
l(e'* X e)c of the polarization vectors, respectively. The 
total scattering amplitude in the dipole limit can now be 
written as 



^/,(e,e',L.k) = ^T-(e,e')-F^(zk) 



(16) 



a;=0 



The spectra for different x and q are combinations of 
the spectra for different polarizations. Usually, the scat- 
tering amplitudes are calculated in terms of the com- 
ponents Da of the dipole operator, where a = 1,0,-1 
in spherical symmetry or a = x,y,z in Cartesian coor- 
dinates. The spectra for different polarizations are then 
combined to form the fundamental scattering amplitudes. 
This can be compared with x-ray absorption. The circu- 
lar dichroic spectrum (the x = I fundamental spectrum 
for x-ray absorption) is usually calculated by subtracting 
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the spectra for left and right circularly polarized light 
{a = ±1). The scattering amplitudes in terms of the 
components of the dipole operator are given by 



{f\DlGiz^)D^ \g) 

{f\Dl\n) {n\D^\g) 
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(17) 



Note that again there are only nine spectra and this is 
just a representation of the nine fundamental spectra in 
a different basis. The simplest scattering amplitude is 
the isotropic one given by x = 0. The tensor contain- 
ing the isotropic scattering amplitudes F° has only one 
component F^, 
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(18) 



which is just a sum of all the different polarization com- 
ponents. For the expressions in spherical symmetry, note 
that, since r|^, = (—1)" ri _a', there is no net transfer of 
angular momentum to the system for the isotropic scat- 
tering amplitude. Since the angular dependence is given 
by T° — e'* ■ e, the isotropic contribution to the spectral 
line shape is removed in many experiments by the use 
of a 90° scattering condition with the incoming polariza- 
tion vector in the scattering plane (tt polarized). This 
makes the incoming polarization vector perpendicular to 
both possible outgoing polarization vectors and therefore 
e'* • e = 0. In addition, this has the advantage that it 
strongly reduces the non-resonant term from the ex- 
perimental RIXS data (which has the same polarization 
dependence). This contributes mostly to the elastic line 
and is frequently the major experimental impediment to 
measuring low-energy excitations. 

Tensors of rank x = 1 have three components. For 
example, the q — component is given by 



xy 



p 



(19) 



For resonant elastic x-ray scattering, the Fq scattering 
amplitude is the one that gives rise to, amongst others, 
magnetic scattering. The angular dependence for x = 1, 
is given by an outer product = e'* x e. 

At this point it is useful to make a comparison with 
X-ray absorption (XAS) and resonant X-ray scattering 
(RXS), which are determined by the scattering amplitude 
J-gg where 

-?"xAs(e,Wk) -ilm[J'gg(e, e,Wk)] (20) 
/Rxs(e,e',wk) = |Jgg(e,e',Wk)P, (21) 

where for X-ray absorption, there is only a polarization 
vector for the incident X-rays, and e' = e. Since for XAS 
and RXS the "final" state is equivalent to the ground 
state in the scattering amplitude (|/) — \g)), an addi- 
tional restriction is imposed on the scattering. In many 
symmetries, this means that only the q = component 



contributes, reducing the scattering amplitude determin- 
ing XAS and RXS to 



^gg{u;^) = Y.TSie,e')FS{z^) 



(22) 



x=0 



This implies that of the 3x3 = 9 components in the full 
scattering amplitude, only 3 components, corresponding 
to a; = 0, 1, 2 and q — 0, remain. For X-ray absorp- 
tion, these correspond to the well-known isotropic, circu- 
lar dichroic, and linear dichoic spectra, respectively. 



3. Scattering amplitude for a multipole expansion 

We next generalize the ideas from the previous section 
to include the different types of multipoles arising from 
the p • A interaction in Eq. Q. Since the dipolar and 
quadrupolar transitions in RIXS are predominantly exci- 
tations from a localized core hole into the valence states, 
the common approach is to expand the plane wave in 
the vector potential, see Eq. ([2|, around the site where 
the absorption takes place. Essentially, one is using an 
approximation of the type e*'^ ''' = H-ik-r; but in spher- 
ical harmonics. In the previous section, we treated the 
case that e'^ '^' ^ 1. The plane wave can be expanded 
in terms of spherical harmonics Yim{0,(p) and spherical 
Bessel functions ji ( Varshalovich et al. 1988| ) 



oo / 
1=0 m=-l 



In order to arrive at the standard transition operators, it 
makes sense, at this point, to rewrite the above equation 
in terms of spherical tensors, see Appendix. A common 
set of tensors are the normalized spherical harmonics, 
which we write as the tensor f with components fm = 
^/4TT/[l]Yim{0r,(Pr)- In addition, r^) = r'fC'. Note that 
i-(o) — I Pqj. spherical harmonic tensors of order 1 = 1, 
the superscript is dropped r = r^'^. This allows us to 
rewrite the expansion as 



(23) 



/=o 



As in the previous section, we want to separate the mo- 
mentum and polarization vectors of the x-rays (the geom- 
etry of the experiment) from the transitions in the mate- 
rial under consideration (the fundamental spectra). This 
can be done by recoupling the different tensors. Recou- 



of the tensors using Eq. ( 
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) leads 
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to 
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^kfc'[P.,rf)]^-[e,kW]^(24) 



where the approximation ji{kr) = {krY for fcr ^ 1 
has been used, with the double factorial III = l{l — 2) ■ ■ ■ . 
Note that the operators acting on the electrons, namely 
the momentum and position r^, are coupled together 
to form an effective operator D'-^ of rank L — I — 1,1,1 + 
1. The quantities related to the photons, namely, the 
wavevector k and polarization e also form a tensor of rank 
L. Let us first consider the transition operators in Eq. 



(24) by introducing the transition operators (Van Vee- 



nendaal and Benoist 1998 1 



D- = ^f:[p,:,rf)]^ 
imuj ^ — ' 



(25) 



with piL{k) = [lL]nli^k'/[l]l\ = l,fc/2,fc/6 for IL = 
01, 11, 12, respectively. The values of I and L give rise to 
the usual dipolar {IL = 01), magnetic dipolar (11), and 
quadrupolar (12) transition operators. For I — 0, one has 
rf'^ = 1 and the operator simplifies to D"^ = J2i 'Pi/i'muj, 



which is equivalent to the dipole operator in Eq. ( 13 1 in 
the previous Section. In cartesian coordinates, the op- 
erator [ri,pi]^ = Li and D^^ = ^f^X^i^i, with oq the 
Bohr radius and the angular momentum given in h. The 
orbital moment forms, together with the Zeeman term in 
Eq. ([I]), the magnetic dipole transition. Since magnetic 
dipole transitions are of the order of = 1/137^, i.e. 
about five orders of magnitude, smaller than the electric 
dipole transitions of the same wavelength, they will be 
neglected in the remainder of the paper. The next opera- 
tor is D^^ = I X^ili ''i^'' which is the electric quadrupole 
operator. 

In the remainder, we limit ourselves to the electric L- 
pole transitions, and we can drop the I = L — 1 from 
the expressions, i.e. D^^^'-^ — > D^. The transition 
operators are then V ^ [e,k(^^i)]-^ • with = 

Phik) 'Yl!i=i ''i^' with L = 1, 2 for dipolar and quadrupo- 
lar transitions, respectively. The relative strengths of 
the components of the multipole transition operators r^^'' 
depend on the direction of polarization and wavevector 
through [e, k*^-^"^)]-^. These reduce to e and [e,k]^, for 
the electric dipolar and quadrupolar transitions, respec- 
tively. As discussed in the previous section, the part that 
depends on the geometry of the experiment and the fun- 
damental spectra (9 and 25 for dipolar and quadrupo- 
lar transitions, respectively) that describe the physical 
properties of the system can be separated exactly. This 
can again be achieved by applying a recoupling using 



Eq. (A3 1 on the scattering amplitude, which can then be 



rewritten as 



2L 



Ffg{\<i, k', e, e', Wk) - ^'^"(k, k', e, e') ■ F^-(fc, k' , uj^). 



x=0 



Neglecting interference effects between different multi- 
poles, the scattering amplitude for a particular multipole 



is given by (Van Veenendaal and Benoist 1998) 



F^-(fc,fc',u;k)=E 



[(/|(D^)t|n),(n|D^ 



hujlr + En 



IL ^ 

which has angular dependence 

T^-(k,k',e,e') = [x]ni^J[e'^k'(^-l)]^[e,k(^-l)]^]^ 

The above equations give an exact separation of the 
Kramers-Heisenberg expression for RIXS into an angu- 
lar dependence and fundamental scattering amplitude 
achieving the first step in Table |T] 

B. Direct and Indirect RIXS 



In the previous section the Kramers-Heisenberg expres- 
sion for the RIXS scattering amplitude J^fg , Eq. |9j was 
derived and re-expressed as a product of a photon ab- 
sorption operator T), the intermediate state propagator 
G and a photon emission operator 2?^, sandwiched be- 
tween the RIXS final and ground state: 



{f\D^GV\g) 



(26) 



The presence of the intermediate state propagator is 
what makes the theory of RIXS complicated - and in- 
teresting. The propagator G is defined in terms of 
the inverse of the total Hamiltonian H of the system, 
G{z-k) = (zk — H)~^ , where the operator H naturally di- 
vides into the ground state Hamiltonian (governing 
the quantum system without a core-hole) and the core- 
hole Hamiltonian He perturbing the system after photon 
absorption: H = Hq + Hq ■ It should be noted that even 
if one commonly refers to He as the core-hole Hamilto- 
nian, it also includes the interaction between the electron 
excited into the conduction band and the rest of the ma- 
terial. As core-hole and excited electron together form 
an exciton, their separate effects on the system cannot, 
in principle, be disentangled . 

At this point it is useful to separate the full propagator 
G into an unperturbed one Go — (zk — Ho)~^ and a term 
that contains the core- hole Hamiltonian He, using the 
identity G = Gq + GqHcG. This also separates the RIXS 
amplitude into two parts, which define direct and indirect 



RIXS (Van den Brink and van Veenendaal 2006): 



■direct 
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(27) 
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Intermediate State Propagator 


Definition 


RIXS amplitude Jig = {±\V'^GV\g) 
Intermediate State Propagator G(-Sk) = ^ — ' ^ ~ -^c 
Ho ground state, He core-hole Hamiltonian 


Exact Propagator 


G = Go + GqHcG 

^ \ 


Direct RIXS 


Indirect RIXS 

jr^nMrect ^ (f | JC^^GP | g) 




Approximation 
to propagator 


Fast Collision'"' 
G(zk) -i> 1/zk 


Perturbation Expansion'""^' Ultrashort Core-hole Life-time^'^^-'' 
GqHcG — > GqHcGo GqHcG — > GqHcGc 


Valence excitation 


caused by transition operators D 


caused by core-hole interaction He 



TABLE II Theoretical approach to the intermediate state propagator, classifying direct and indirect RIXS processes and 
common approximations to the propagatoiQ 

°( |Luo et ai ! || l993'). ''('Platzman and Isaacs' '1998), "(Doring et aZ.''2004'), (Abbamontc et ai." 1999a'), (N omura and I garashT 2004'), 
f ' Nomura a nd Igarashi, ,2005J, ^(iMarkiewicz and Bansil, ,2006| ), "( |Van den Brink and van Veenendaal, |2005| , 'i jVan den Brink and van] 
Veenendaall 1200611, HAment et al||2007||. 



and 



-pindirect 



{f\V^GoHcGV\g) 



(28) 



Note that this definition of direct /indirect RIXS, based 
on the Kramers-Heisenberg expression, is exact. 

For the direct RIXS amplitude, the core-hole does not 
play a role - the photon absorption and emission matrix 
elements determine which electronic transitions are al- 
lowed. The physical picture that arises for direct RIXS 
is that an incoming photon promotes a core-electron to 
an empty valence state and subsequently an electron from 
a different state in the valence band decays, annihilating 
the core-hole, see Fig. [7j Thus for direct RIXS to oc- 
cur, both photoelectric transitions - the initial one from 
core to valence state and the succeeding one from valence 
state to fill the core hole - must be possible. These tran- 
sitions can, for example, be an initial dipolar transition 
of Is 2p followed by the decay of another electron in 
the 2p band from 2p — >■ Is. This happens at the X-edge 
of oxygen, carbon and silicon. In addition, at transi- 
tion metal (TM) L-edges, dipole transitions causing di- 
rect RIXS are possible via 2p — >■ 3d and 3c? 2p dipolar 
transitions. In all these cases RIXS probes the valence 
and conduction states directly. 

For indirect RIXS, the scattering amplitude depends 
critically on the perturbing core-hole Hamiltonian - with- 
out it the indirect scattering amplitude vanishes. In 
general the scattering, j:™direct ^ arises from the com- 
bined influence of He and transition matrix elements T). 
Most often for indirect RIXS, V/V'^ create/annihilate an 
electron in the same state, far above the Fermi level. 
For instance at the TM if-edge, the Is Ap process 
creates/annihilates an electron in Ap states electonvolts 
above the TM 3d valence shell. The delocalized Ap 



electron can then be approximated as being a spectator 
because (Coulomb) interactions involving the localized 
core-hole are usually much stronger and dominate the 
scattering cross section. 

It should be noted that if scattering is direct, as for 
instance at TM L-edges, indirect processes can also con- 
tribute to the total scattering amplitude. However, as 
indirect scattering arises in this case as a higher order 
process, it is normally weaker than the leading order di- 
rect scattering amplitude. Conversely, in case of indirect 
RIXS, direct processes are absent by definition. 



IV. THEORY OF RIXS: II. CALCULATING THE 
CROSS-SECTION 

Computing the exact RIXS response of a material us- 
ing the full propagator is in general a daunting if not im- 
possible task. This section describes the different approx- 
imate methods used to evaluate the Kramers-Heisenberg 
scattering amplitude for a specific system. 

The most straightforward approach is to evaluate the 



exact scattering amplitude Eq. ( 12 ) and intensity Eq. (18 



numerically for a simplified system. One can for instance 
neglect or simplify correlation and core-hole interaction 
effects and resort to an effective single-particle descrip- 
tion, for instance with band structure methods. Another 
possibility is to keep the correlation effects but reduce 
the number of degrees of freedom by considering instead 
of the full macroscopic system only a small cluster of 



atoms (Kotani and Shin 2001). We will briefly discuss 



such cluster methods in Section ITV. Al 

A complementary approach is to approximate the 
intermediate-state core-hole propagator. At first, this 
appears less appealing than a direct evaluation of the 
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RIXS cross section. However, one should note that, in 
particular for strongly correlated systems, numerical cal- 
culations can only be done for small systems, such as an 
ion, a cluster of atoms, or an impurity model. Often, one 
is also interested in the coupling of RIXS to collective 
excitations, not present in these small systems. There- 
fore, effective methods that project out the intermedi- 
ate states are of crucial importance for understanding 
the RIXS cross section. These approaches provide effec- 
tive response functions for RIXS which can be related 
to response functions for other spectroscopic techniques. 
As a next step one can then compute, in cases where it 
is feasible, the effective RIXS response function exactly. 
Depending on whether one considers direct or indirect 



RIXS, one uses Eqs. (27) or (28), as a starting point to 
derive an effective expression for the RIXS cross section. 

For direct RIXS the transition operators are key to 
the scattering process. One can highlight this physics by 
neglecting the dependence on He of the full intermedi- 
ate state propagator, so that G — >■ 1/zk - this amounts 
to neglecting the indirect contributions to the scattering 
process. This is known as the fast collision approxima- 
tion (Luo et al. 1993 Van Veenendaal 2006). This ap- 



proximation allows the reduction of the full direct RIXS 
amplitude to a rather elementary response function of 
effective scattering operators. This is a great advantage 
and simplification, even if evaluating this response func- 
tion for a real (strongly correlated) material still presents 
a challenge in itself - as is evaluating any other response 
function for a correlated electron system. These effective 
theories for direct RIXS are reviewed in Section lIV.BI 

In the case of indirect RIXS the scattering is due to 
the intermediate state core-hole potential. The approx- 
imation closest at hand is to treat this perturbing po- 
tential as being weak, expand G in a perturbation se- 



ries and retain only the leading terms ( Abbamonte et al. 



1999a'; 'Boring et al, 2004; MarkicwiczandBansil, 20W. 
Nomura and Igarashi, ,2004, ,2005, ,Platzman and Isaacs^ 



1998 Takahashi et al. 2007). To deal with the strong 



core-hole interaction, the Ultrashort Core-hole Lifetime 
(UCL) expansion of the core-hole propagator was devel- 



oped (|Ament et al. 


2007 


Van den Brink and van Veenen- 


daal 2005 


2006 


), which treats the core- hole potential as 



the dominant energy scale. Section [IV. C| presents a sur- 
vey of the different effective theories for indirect RIXS. 



A. Numerical Evaluation of Kramers-Heisenberg 

There are several numerical approaches used to evalu- 
ate the Kramers-Heisenberg equation directly. Although 
there are many different methods to calculating the elec- 
tronic structure, the general approach is always the 
same. First, one finds the eigenstates and eigenvalues 
for ground, intermediate, and final states. Subsequently, 
one calculates the relevant electric-multipole matrix ele- 



ments between the core levels and the valence shell for 
the appropriate scattering condition, and evaluates Eq. 
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1. Single particle approach 

A common approach in the description of the electronic 
structure of solids is the use of an independent-particle 
approximation. After applying a Local-Density Approxi- 
mation, the electronic structure can be calculated ab ini- 
tio. Numerical methods based on this approximation can 
be successfully applied to systems where the interaction 
between the electrons can be described by an effective 
potential. For RIXS, this approach is particularly useful 
in the study of semiconductors and wide-band insulators 
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Jia et al. 
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1994 Kokko et al. 
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Here, the absorption process excites an electron from 
a deep-lying core level into the bands above the Fermi 
level. Varying the energy of the incoming photon en- 
ergy over the absorption edge essentially probes the un- 
occupied density of states. The de-excitation removes an 
electron from the states below the Fermi level. The final 
states have an excited electron-hole pair. Since the trans- 
ferred energy and momentum are known, the RIXS spec- 
trum is essentially a momentum-resolved convolution of 
the projected occupied and unoccupied density of states. 
For relatively shallow core levels such as the Is level of 
Carbon, the momentum transfer is very small and the 
created electron-hole pairs have g = 0. 

This approach entirely neglects any interaction that 
can occur in the intermediate state between the core- 
hole and valence-shell electrons. It is well known that for 
semiconductors or large band-gap insulators, the attrac- 
tive Coulomb interaction between the core hole and the 
valence electrons gives rise to an the intermediate state 
with an exciton consisting of the core hole and the elec- 
tron excited into the conduction band. Since this is a two- 
particle problem, it can be solved exactly using Bethe- 



Salpeter type equations (Shirley 1998 Van Veenendaal 



and Carra 1997). In a metallic system, the absence of 
a gap leads to the creation of low-energy electron-hole 
pair excitations that screen the core-hole potential. This 
response leads to X-ray edge singularities in the RIXS 



hams 1974) 



spectral line shape ( Ahn et al. 2009 Nozieres and Abra- 



2. Exact diagonalization and other cluster methods 

In recent years, the majority of RIXS experiments 
have been performed on materials that display strong 
Coulomb interactions between the electrons. For these 
systems, ab initio methods based on the assumption of 
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independent particles are not applicable and one gener- 
ally resorts to model Hamiltonians. In these calculations 
one defines an appropriate basis set. Subsequently one 
sets up the Hamiltonian using parameters obtained from 
ab initio methods, Hartree-Fock type calculations, or in 
a more heuristic fashion. The Hamiltonian is then diago- 
nalized for ground, initial, and final states and the RIXS 



cross section is calculated using Eq. (17) 



Exact diagonalization methods have the advantage of 
being able to treat many-body effects since the eigen- 
states are not limited to a single Slater determinant. The 
major disadvantage is that the basis set increases very 
rapidly since all possible permutations of electron in the 
basis set need to be considered. Despite these hmitations, 
these approaches have been successful in describing the 
RIXS cross section, in particular for the L and M edges 
in transition-metal compounds. This has been described 



extensively in Refs. (|Agui et at. 
eTaT, 2005 



2005 



2008 : |De Groot et al.[ 

2004; GhiringhcUi et aZ.V2005V2008' 
2000: Huotari et at 



2009 



Chiuzbaian 



1998; Ferriani et al.[ 



2006b; Harada et aV 
_ 2008 ; Magnuson et al. , 2002a b ; Mat- 

subara et a/.||2005||Van Veenendaal and Fedro(|2004]fVeF 



nay and Delley 2009) and reviewed in Refs. (De Groot 



and Kotani 2008) and (Kotani and Shin 2001) 



Since the many-body basis set that one can handle nu- 
merically is limited, it is essential to decide which states, 
interactions and electronic configurations are to be taken 
into account, and conversely, which ones can safely be 
neglected. 

In direct RIXS, after a np ^ 3d (with n ~ 2, 3) dipo- 
lar transition, the excited electron is strongly bound to 
the np core hole. Therefore the most important con- 
figurations are the metal-centered ones. The process can 
then be described as n/3d™ np^Sd^+i ^ rip*'(3d™)*, 
where m are the number of 3d electrons and (3^™)* in- 
dicates an excited 3d state. Since L and M RIXS in 
transition-metal compounds is mainly a direct process, 
many excited states can be accessed. 

In order to interpret the dd type excitations that are 
typically observed at these edges, it is important that the 
Hamiltonian include a Coulomb interaction that includes 
higher-order terms. This interaction splits the 3d states 
into Coulomb multiplets that can be separated by tenths 
of eV to several eV (Griffith 1961). A well-known effect 



of these higher-order terms in the Coulomb interaction is 
Hund's rule that determines the ground state for a free 
ion. To mimic the effect of a solid, a crystal field can be 
introduced, which further splits the Coulomb multiplets. 
An additional interaction that needs to be included is the 
spin-orbit coupling, in particular since at the L edges in 
transition-metal compounds, the 2p spin-orbit coupling 
is by far the strongest interaction in the intermediate 
state, splitting the absorption edge into two distinct fea- 
tures. Since the direct RIXS process creates excitations 
involving all 3d states, it is in addition very important to 
properly include dipole/quadrupole transitions, since po- 



larization effects are crucial in the interpretation of the 
RIXS at the L and M edges in transition-metal com- 
pounds. 

A next step in including solid-state effects is to take 
into account the ligands, such as oxygen, surrounding 
the ion where the scattering takes place explicitly in 
the calculation. Usually, a metal ion is surrounded by 
4-6 ligands. The Hamiltonian then needs to include 
the hybridization between the metal and the ligands 
and the on-site energy of the ligands. An important 
parameter in this respect is the charge-transfer energy 
A = i;(3d™+iL) - £;(3d'"), where L indicates a hole 
on the ligands, that gives the energy needed to transfer 
an electron from the ligands to the metal site. This ap- 
proach is known as (small) cluster calculations, referring 
to the cluster of atoms taken into account in the exact 
diagonalization. Instead of including only the surround- 
ing ligands, one can also include the bands corresponding 
to the ligands. This is known as a single-ion Anderson 
impurity model. However, often it is more important to 
include additional metal sites in the small cluster cal- 
culations, especially when considering strongly covalent 
and/or doped materials. 

For indirect RIXS, other considerations are often more 
important in limiting the size of the basis set. Here, va- 
lence excitations are shake-up structures related to the 
screening of the strong transient core-hole potential in 
the intermediate state. Since the screening processes of- 
ten only involve the 3d states closest to the Fermi level, 
not all 3d orbitals are necessarily relevant. On the other 
hand, the spatial extent of the cluster becomes increas- 
ingly important, since screening charge is pulled in, not 
only from neighboring ligands, but also from neighboring 
metal sites. Furthermore, since the excited electron can 
usually be considered a spectator, the transfer of angu- 
lar momentum to the valence band is much smaller and 
polarization effects are less important. 

Therefore, the appropriate reduction of the basis set 
depends strongly on the type of final state excitations 
that one considers. Let us illustrate this with an ex- 
ample: consider a divalent copper compound, which has 
holes in the 3dx2-y2 orbital. In describing the L and 
M edge RIXS it is important to include all the 3d or- 
bitals and their appropriate crystal fields, since the final 
states consist primarily of local dd transitions. In ad- 
dition, spin flips can occur due to the strong core-hole 
spin-orbit coupling in the intermediate state. In addi- 
tion, since this is a direct process, polarization effects are 
important in reproducing the correct RIXS intensities. 
At the K edge, on the other hand, the screening dynam- 
ics in the intermediate state primarily involve the Zd^^^-yi 
orbitals. If the calculation is restricted to a small cluster 



containing only one copper ion (Hill et al. 1998), only 
high-energy bonding-antibonding between the 3d'" and 
g^m+i^^ configurations on a local copper plaquette are 
included. This misses out the important charge-transfer 



25 



excitations between different copper plaquettes. In order 
to describe tfiese and obtain insigfit into tlieir dispersion 
additional simplifications can be made by reducing the 
copper-oxygen plaquette to an effective single site sys- 
tem. This leads to the well-known single-band Hubbard 
model, where a number of exact diagonalization studies 



for RIXS have been carried out (Tsutsui et al. 19991. 



B. Effective theory for direct RIXS 

The assumption in direct RIXS is that the excitations 
are predominantly created by the radiative transitions. 
Essentially, this approach was followed in papers on sys- 



tems with wide bands such as boron nitride (Jia et al. 



1996), diamond (Ma et al. 19921, and graphite (Carlisle 



et al. 1995). In these systems, the intermediate state 
propagator (or, equivalently, the x-ray absorption spec- 
trum) can be described as the ground state plus an elec- 
tron excited into the unoccupied states. The RIXS spec- 
trum is then interpreted as the creation of electron-hole 
pairs and effectively a (projected) joint density of states 
due to the two consecutive radiative transitions. The po- 
larization dependence of the direct RIXS process shows 



up nicely for graphite ( Carlisle et al. 1995 1 in the change 



in intensity as a function of the angle of the polariza- 
tion of the incoming light with respect to the graphite 
planes. However, for strongly correlated systems, such as 
transition-metal compounds, such an intuitive approach 
to the intermediate states often fails due to intermedi- 
ate state interactions. However, at the transition-metal 
L and M edges, the assumption that the excitations are, 
to first order, due to the dipolar transitions is still valid. 
In order to understand what valence excitations can be 
created by the radiative process, one can apply the UCL 
approximation to lowest order (also known as the fast- 
collision approximation). In this limit, the intermediate- 
state propagator is replaced by a resonance factor. 



1. Fast collision approximation and effective transition 
operators 

To obtain effective valence operators for the L and M 
edge, the dipole operator is expressed in second quanti- 
zation ( Judd 1967 ) . For a transition between a core level 
and the valence shell with angular momentum c and /, 



respectively, we have 



ncC,nL 



(29) 



where the ^/\cl] has been kept separate for use in the 
derivation of sum rules and effective operators. The re- 



duced matrix elements are contained in 

/ L c 



P' 



= PL-lx{k){-^) 







xriciL/ drr'+''Rni{r)RnjSr), (30) 

where the coefficient naix is defined in the Appendix; 
The term in parentheses is a 3j symbol, which is propor- 
tional to the Clebsch-Gordan coefficients. The integral is 
over the radial parts of the wavef unctions. Although it is 
important in determining the strength of the transition, 
in general, it only scales the spectrum and we can often 
take Pni,ncc = 1- The essential physics is contained in 
the operator 



E 



(31) 



where Cma- annihilates an electron in the core level with 
orbital component m and spin a and ll^,^ creates an elec- 
tron in the valence shell with quantum number m' and 
the same spin. The coefficients are given by 



w, 



ll'L 
mm' M 



I' 



L I 
M m 



"ll'L- 



(32) 



The 3j symbol imposes the well-known selection rule 
m' = m -\- M since the sum of the elements in the lower 
row has to equal zero. 

In direct RIXS, the spectral line shape is mainly de- 
termined by the transition operators, see Eq. (27 1. In 



order to obtain effective operators that give the com- 
bined action of the absorption and emission process, we 
need to approximate the intermediate-state propagator. 
In order to show the procedure, let us first consider a 
single edge, before considering the situation where the 
edge is split by the core-hole spin-orbit coupling. This 



is usually done using a fast-collision approximation ( Luo 



et al. 1993 Van Veenendaal 2006 ) 



f;^,,(z) = (/I [Dt, G(zk)D]^ \g) = G(z^) (/I [Dt,D]- \g) , 

where G{z)^) — (zk — E)^^ . This approximate interme- 
diate state propagator is then effectively the resonance 
function P(a;k), see Table |lj describing the strength of 
the scattering amplitude as a function of the incoming 
photon energy. The dipole operator can be written as 



where the reduced matrix element, 
which only scales the spectrum, has been taken unity, 
Pn^Cjui = 1- Since neither the ground state nor the final 
states contain a core hole, we can effectively remove the 
creation and annihilation operators related to the core 
states (Luo et al. 1993 Van Veenendaal 2006). Using 



spherical tensor algebra, it can be shown (Luo et al. 
2U06| 



1993 Van Veenendaal 



G(zk)[d](/|[(w^'^)t,W^'^] 
G(zk)[c](/|«;f 



15) 



(33) 
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where the usual situation I = c + L has been taken; w°'^ 



is a tensor containing the operators w'^^ with M 
1, • • • , —L. Equation (33 1 is the central expression for the 



fundamental scattering amplitude in the absence of core- 
hole spin-orbit coupling. Note that the combined action 
of the two operators w'^'-^ representing the absorption 
and emission between the core level and the valence shell 
has been replaced by one effective hole operator 



- Ill 



E 



"^mm' M'"in' af-rna 



(34) 



acting only on the valence shell. The coefficients are 
again given by Eq. (32). For a; = 0, the operator 
= rih is the number of holes operator. Within an 
ionic framework, this only contributes to the elastic line, 
although it can give rise to inelastic satellite structures. 
For X = 1,2, the operator is the angular momentum 
w''^ = L/Z and the quadrupolar moment, respectively. 
The transitions made by these operators will be discussed 
in the Section on dd excitations. 

In many cases, the core-hole spin-orbit coupling cannot 
be neglected. For example at the L23-edges (2p — )■ id) in 
transition-metal compounds, one has two distinct edges 
corresponding to a total angular momentum — c± ^ — 
|, I of the core hole with angular momentum c—1. The 
same evaluation of the combined action of the transition 
operators for absorption and emission can be done at a 
particular spin-orbit split edge. Restricting ourselves to 
the common situation where I = c + Q, the expressions 
can be reduced to 



, W-' 



l\9) 



TABLE III Orbitals for the s, p, and d orbitals with / = 
0, 1, 2, respectively. The value oijj, — l,l—l,--- , —I. The sym- 
metry notations in Oh (octahedral) and are also given. 



orbital 


/ 




Oh 


Dai, 


s 








Ai 


^1 


Px 


1 


1 


Ti 


E 


Pz 


1 





Ti 


A2 


Py 


1 


-1 


Ti 


E 




2 


2 


E 


Bi 


dzx 


2 


1 


T2 


E 




2 





E 


Ai 


dyz 


2 


-1 


T2 


E 


dxy 


2 


-2 


T2 


B2 


dependence and can 


assume 


the values and 1, cor- 



responding to spin-independent and spin-dependent op- 
erators respectively; x = z — + l relates to the 
orbital part of the operator. First, note that for the 
spin-independent operators = and that when 
adding the transition operators from the different edges, 

i^q^{j^) +Y-q^{j~) = [c]iq°^, the rcsuh in the absence 
of spin-orbit coupling is reproduced. Therefore, the first 
term on the right-hand side in Eq. (36 1 gives, at each 
edge, the transition as if there was no spin-orbit interac- 
tion. Note that the intensity at each edge is proportional 
to the multiplicity of the core level -I- ^ . The additional 
operators are spin-dependent. Most notable are the spin 
v"" = 2S and the spin-orbit coupling vl'^° = fL'S. The 
effect of the spin-dependent operators will be discussed 
in the Section on magnetic excitations. 



G,(zk)(/|Z, {f)\9), 



(35) 



where the fast-collision approximation is now done for 
a particular spin-orbit split edge Gj{z\^) = (zk — Ej)~^, 
where Ej is an average energy for each edge with j = 
Z ± ^. The rank of the total tensor is now z (reserving 
the X for the orbital part of the operators). The operator 
is now 



- llz . I , 
Vq U^) = U 



± 



1 



hq 



iz+l){z+2) 



z+l,lz 



(36) 



with = c ± Since the spin is no longer a good 
quantum number at a particular spin-orbit split edge, 
the operators are now spin dependent 



r,xvz _ 



(37) 



where S° = 1 and — S. The II in the superscript 
has been assumed implicitly to avoid crowding of in- 
dices. The rank z = 0, • • • , 2L is still related to the 
rank of the transition operators; y indicates the spin 



2. Selection rules 



Equation (33) shows that, within the fast-collision ap- 
proximation, the combined action of the transition op- 
erators for absorption and decay can be written as an 
effective transition operator in the valence shell. The 
rank of the transition operator is x. For dipole opera- 
tors, we have x = 0,1,2. The presence of the x = I 
term might seem surprising since this is comparable to 
an effective dipole operator (due to r, a tensor of rank 
1) in the valence shell, and it is well known that dipole 
transitions within one shell are not allowed. However, 
such transitions arise from the reduced matrix element 



in Eq. ( 30 ) and not from the transition operator in Eq. 
(32). Since the reduced matrix element in RIXS is the 
product of the reduced matrix elements for absorption 
and emission, the x = 1 term is allowed for transition- 
metals at the transition-metal L and M edges. For the 
components, we have the well-known selection rule 



m = m + q. 



(38) 



27 



When performing RIXS at a particular spin-orbit split 
edge the effective operators in the valence shell can be 
written in terms of spin-dependent operators, see Eq. 
(35). The selection rule is then effectively modified to 



(Van Veenendaal 2006) 



(39) 



When X — 0, there is no net transferred angular momen- 
tum to the system {q ~ 0). In the absence of spin flips, 
the effective operator is the number of holes Wq°° = Uh 
which only contributes to the elastic line. However, when 
considering spin-dependent operators, inelastic features 
can be observed. The spin-dependent x = operator is 
the spin-orbit coupling = (2/?)L • S. Inelastic scat- 
tering in spherical symmetry is due to the L±S:f: term, 
where the change in angular momentum is compensated 
by the change in spin. 

Many materials are more appropiately described in 
a lower symmetry than spherical. Often, the basis set 
used for octahedral symmetry is used for transition-metal 
compounds. Unfortunately, the standard notation is not 
very conducive for the definition of selection rules. A 
more appropiate notation is given by the use of tesseral 
harmonics 



3. Sum rules 

Sum rules relating the integrated spectral weight to 
ground-state expectation values of such quantities as 



the orbital and spin moment (Carra et al. 1993 Thole 



et al. 1992) have played a major role in x-ray absorp- 



tion and x-ray magnetic circular dichroism studies. The 
derivation of such rules is significantly more complex for 
higher-order spectroscopies, due to the presence of the 
intermediate-state propagator, and approximations are 
necessary. Using the fast-collision approximation, sum 
rules can be derived in the situation where the decay is 



between two core levels ( Carra ei aZ. 1995 Van Veenen- 



daal et al. 1996a). Since the decay process can then be 



decoupled from the excitation process, the sum rules re- 
late the integrated RIXS intensities to expectation values 



of one-particle operators (Carra et al. 1995 Van Vee- 



nendaal et al. 1996a). When the decay involves the 



valence shell, the RIXS sum rules involve two-particle 
operators which are difficult to interpret. When the po- 
larization of the outgoing x-rays is essentially isotropic, 
numerical evaluation of the integrated intensities shows 



(Van Veenendaal et al. 1996b) that the sum rules are 



close to those for x-ray absorption (Carra et al. 1993 



Thole et al. 1992). 



where, with Sp 




'^V' (r). 



(40) 



±1 



1 



(i+«)/2 and B,. 



for 
= 0, 



0, A, 



for = 



and > 0, respectively. Roman and Greek letters are 
used for the components of the angular momentum in 
for spherical and tesseral harmonics, respectively. The 
tesseral harmonics correspond to the well-known orbitals, 
see Table 



III 



The operators can be transformed to 
VFJ'^ in a tesseral basis via a unitary transformation, 
see Appendix. The selection rule for transition between 



2010) 



shells I and I' in Eq. ( 39 1 then becomes ( Van Veenendaal 



(41) 



respectively, and ti = ±1 for Z is even/odd. In addition, 
the transition to /x' = is not allowed when Sij,s^ti+i'+x < 
0. A notable example of the difference between is the 
transition operator w^'^^ — —Lq/1 = —L^/l. In spherical 
symmetry, this gives a transition m' = to, i.e. the orbital 
does not change. However, for tesseral harmonics, one 
has ^' = — /i, and therefore scattering — O xy and 
zx O yz. Obviously, the selection rules should give the 
same results as those that can be obtained from general 
group-theoretical considerations, see TablelTTT] Examples 



of the use of these selection rules are given in Section V.D 



C. Effective theory for indirect RIXS 

In the previous section, we have seen that for direct 
RIXS the lowest term in the UCL expansion (also known 
as the fast-collision approximation), can be written in 
terms of effective transition operators corresponding to 
the states into which the absorption takes places. How- 
ever, often for indirect RIXS these terms do not con- 
tribute to the RIXS cross section. For example, i^T-edge 
RIXS is dominated by excitations into the transition- 
metal 4p states. Since the Ap states are usually almost 
completely empty, the effective operators for direct RIXS 
only contribute to the elastic line, where the effective 
transition operator creates an electron in the valence shell 
in the excitation step and annihilates it again in the emis- 
sion process. 

Experimentally, however, RIXS is observed at the K- 
edge. Particularly prominent are the charge-transfer type 
excitations that will be discussed in Section IV. Bl Also 
the excitation of dd transitions and magnons have been 
observed, see Sections |V.D| and |V.E| The general con- 
sensus is that these excitations are created through the 
interaction between the valence shell and the Is-Ap exci- 
tation created in the absorption process. Most work has 
focused on the interaction with the potential of the Is 
core hole, which is known to be of the order of 6-8 eV. 
This potential can be written as 



Ho = 



kk'q/KTCT' 



/Sk'CT'rfkMO-, (42) 
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where fj, sums over the different orbitals. The poten- 
tial can in principle contain exchange terms, but these 
are negligible at the K-edge. The transient presence of 
this potential in the intermediate state leads to strong 
screening dynamics in the valence shell giving rise to the 
final-state excitations. This Section discusses some of 
the methods used to describe the excitations created by 
interactions in the intermediate state. 




Is/Vls 




1. Momentum dependence for indirect RIXS 

Recognizing that for indirect RIXS the core-hole dom- 
inates the scattering process has an important conse- 
quence for the momentum dependence. In the hard X-ray 
regime photons have a momentum q that can span sev- 
eral Brillouin zones because it is larger than the recipro- 
cal lattice vectors G. The photon momentum reduced to 
the first Brillouin zone is, by definition, k = q— nG. The 
translational invariance and localized nature of the core 



potential in Eq. ( 42 1 imply that the momentum depen- 
dence of RIXS is determined by the reduced momentum 
K. It will only weakly depend on nG as in reality a fi- 
nite, but small, length-scale is associated with the core 
potential. RIXS spectra will therefore appear practically 
identical in different Brillouin zones. This is confirmed 



experimentally by Kim et al. ( 2007 ) . The weak variations 



found in Chabot-Couture et al. (2010), are attributed by 



the authors to polarization effects. 

This is remarkable because in IXS the total momen- 
tum q determines the scattering amplitude. The reason 
for this is that in IXS q enters directly into the transi- 
tion matrix elements, which in RIXS are dominated by 
dipolar transitions for which e*'^ '" = 1 and are therefore 
independent of q. 

In the following, we review how in certain limits the 
indirect RIXS amplitude can be related to the dynamic 
electronic structure factor Sk{uj), which is directly mea- 
sured by IXS. The important difference is thus that IXS 
measures S'q(w) and RIXS is, in these cases, related to 
S^{uj). 



2. Perturbative approach 

The most straight-forward approach to include effects 
of the interaction He between the core hole and the 
valence shell is the use of perturbation theory. This 



amounts to replacing G by Go in Eq. (28) (jAbbamonte 


et al. 


1999a 


|D6ring et al. 


2004 Nomura and Igarashi 


2004 




20051 


Platzman anc 


Isaacs 1998 Semba et al. 


2008 


so that 



-r-indirect 



{f\V^G,HcGoV\9), 



(43) 



which is also referred to as the Born approximation and 
shown in terms of Feynman diagram expansion in Fig. 17 



FIG. 17 Feynman diagram for the transition probability in 
an indirect RIXS process in the Born approximation. Greens 
functions for Cu Is, 4p, and 3d electrons correspond to the 
solid lines labeled Is, 4p, and 3d, respectively. The wavy 
and broken lines represent the photon propagator and core 
hole potential il7is,3d, respectively. The shaded triangle is the 
effective scattering vertex renormalized interactions between 
the valence electrons in the 3d-shell ( [Nomura and Igarashi| 
^0051. 



For dipolar Is ^ Ap transitions at the K-edge, we have 



2? = ^^^lV4pE.ka£ 



with I 



x,y,z. 



and 



^is ip the reduced matrix element containing the integral 
over the radial parts of the wavefunction, see Eq. ( 30 ) . 



In indirect RIXS, one considers the case where the Is- 
4p exciton created in the absorption step is annihilated in 
the emission process. Since there is a momentum trans- 
fer q from the photons to the system, this implies that 
the momentum of the ls-4p exciton must have changed 
in the intermediate state. This can only be a result of 
interactions of the Is-Ap exciton with the valence shell. 
If the dominant interaction is the Coulomb interaction 
of the core hole with the valence shell, then the isotropic 



scattering amplitude can be rewritten as (Doring et al. 



2004 Schiilke 2007) 



-T-indirect 

■^fg 



(q,c^) = P(wk,Wk')r(e,e')(/|Pql5), (44) 



Note that all the operators involving the Is and 4p states 
have been removed from the expression. The density 
operator is 



= E4 



(45) 



kcr 



The resonance behavior is determined by the resonant 
function 



3(-Pls,4p)^ 



s.3d 



(46) 



using Eq. (30) and the fact that hu — hojk ~ huj^' — Ej. 
The resonant function is more complex than for direct 
RIXS reflecting the fact that this is a higher-order exci- 
tation. The polarization dependence requires some care- 
ful consideration. In the situation, where the Ap electron 
is a spectator, an electron is excited into the Ap band 
with momentum k and band index n and subsequently 
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removed from the same state, 

= ^ E Ceo.{0\p^c.'\Et){Et\pijO) (47) 



with a ~ x,y,z. In the atomic hmit, this expression 
reduces to e'* • e, since the orbital of the 4p electron is 
unchanged in the intermediate state. In the absence of 
band effects, a change in the polarization therefore im- 
plies that the angular momentum of the valence elec- 
trons has changed. However, the 4p states form wide 
bands that are mixtures of the different 4p orbitals, and 
these local-symmetry arguments only apply at the F- 
point. Therefore, the use of a scattering condition where 
the incoming polarization vector is perpendicular to the 
outgoing polarization vectors does not necessarily imply 
that a symmetry change has to occur for the valence elec- 
trons. 

The essential physics of the material is contained in 
the fundamental scattering amplitude 



FfM = (/|Pq|5)- 



(48) 



This quantity is directly related to the dynamic structure 
factor, through 



^q(c.)-^|F/3p^(i?^-M 
/ 

— Im(Q|p_n 



iO- 



:Pql.9), (49) 



which corresponds to the bubble in the Feynman diagram 
in Fig. [it] It should be noted that RIXS measures a pro- 
jected S'q(a;), meaning that pq contains only djj^^q ^dko- 
terms. This is a direct result of the fact that the core-hole 
Coulomb interaction does not scatter between different 
orbitals. This is different from IXS, where in principle the 
photon can induce a direct transition from the d states to 
the ligands. This does not imply that RIXS does not cre- 
ate charge-transfer excitations, since the charge-transfer 
states also have d character. In (Igarashi et al. 2006) 



multiple scattering corrections to the Born approxima- 
tion are also considered on the basis of a Keldysh Green's 
function formalism. It was found that multiple scatter- 
ing effects lead to small modifications in the shape of the 
RIXS spectrum, which partly justifies the Born approx- 
imation for wide gap insulators such as La2Cu04 and 



NiO (Takahashi et al. 2007). 



The scattering amplitudes for direct and indirect RIXS 
in Eqs. (lisl) and (iMl) look similar. For direct RIXS, the 



detailed dependence on the polarization is given in the 
fundamental scattering amphtude FJ^ ^. When includ- 
ing the polarization dependence for indirect RIXS, one 
obtains a similar fundamental scattering amplitude. The 
quantity Ffg in Eq. ( 48 ) then reduces to Fj^ q and cor- 
responds to the isotropic term. In the absence of interac- 
tions causing a transfer of angular momentum between 



the Ap and the valence shell, the indirect RIXS ampli- 
tude is simply proportional to pq. In terms of tensors, 
pq = <0(q). 

A significant difference between the two processes is 
that when direct excitations are made into the valence 
shell {e.g. ip/Sp — > 3d), the effect of the operator 
^ddo _ jg relatively small, since the excited elec- 
tron screens the 2p/3p core-hole potential very well. The 
isotropic contribution then mainly contributes to the 
elastic line. In indirect RIXS, the excited delocalized 4p 
electron does not screen the Is core-hole potential very 
well. This produces instead in the intermediate state 
appreciable screening dynamics of the valence electrons 
- the reason why the pq response generates significant 
inelastic scattering intensity for indirect RIXS. 



3. Ultrashort Core-hole Lifetime expansion 

The potential that the core-hole exerts on the valence 
electrons is strong, the attraction Uis,3d between a Is 
core-hole and 3d electron is typically ~ 6 — 8 eV, which 
is of the same order as the d — d Coulomb interaction 
UzdM that appears in Hubbard-like models. Treating 
such a strong interaction as a weak perturbation ren- 
ders a perturbation expansion uncontrolled. To deal with 
the strong core-hole interaction, the Ultrashort Core- 



hole Lifetime (UCL) expansion was developed in ( Ament 
et al.\ |2007| [Van den Brink and van Veenendaal[ |2005[ 
2006), which treats the core-hole potential as the domi- 



nating energy scale. 

The UCL relies on three observations. First for most 
RIXS intermediate states, the core hole lifetime broad- 
ening is quite large: typically F is of the order of 1 eV. 
This yields a time scale r = 1/2F = 4 fs. Only dur- 
ing this ultrashort time, the system is perturbed by the 
core hole. Many elementary excitations have an intrinsic 
time scale that is much larger than 4 fs. This intrinsic 
timescale is the fundamental oscillation period, related to 
the inverse frequency lj of an excitation with energy hw. 
For example, phonons have a typical energy scale up to 
100 meV, and magnons up to 250 meV, thus correspond- 
ing to timescales almost an order of magnitude larger. 
Even low energy electronic valence band excitations can 
be within this range. 

The resulting physical picture of a RIXS process in- 
volving low-energy excitations is therefore that the dy- 
namics in the intermediate state are limited because of 
lack of time, provided that the excitation time scale is 
not decreased significantly by the core hole. The sec- 
ond observation is that the core-hole potential can to 
good approximation be treated as a local potential, i.e. 
its dominating effect is to perturb electrons on the same 
atom where the core hole resides. Finally, the core-hole is 
considered to be immobile, which is a reliable assumption 
for the deep core-states such as Cu Is. 
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The calculation of the indirect RIXS amplitude within 



the UCL expansion by ( 


Ament et al. 


2007 




Van den 


Brink and van Veenendaal 


2005 


2006 


) is based on a 



series expansion of the Kramers-Heisenberg equation, 
(Eq. l9|. But a Green's function approach is equally vi- 



able, which then starts by inserting in Eq. ( 28 1 the iden- 
tity G^Gc + Gc'HoG 



■indirect 
is 



- {f\V^GoHcGc{l + HoG)V\g) 



where the Green functions, Go = (zk — Hq) 



(50) 

Gc = 



{zk—Hc)^^, G = {zk — H)~^, correspond to the Hamilto- 
nian of the unperturbed system Hq, the valence electron 
- core- hole interaction Hq, and the total Hamiltonian 
H = Hq + He- The UCL is best illustrated by con- 
sidering the core-hole Hamiltonian He = Uc X^i PiPfi 
where Uc = Ui^^^d and pY {pI'^) are the density opera- 
tors counting the number of Is core-holes (3c? electrons) 
at site i. The simplest system one can consider is where 
one allows the 3d states to be only occupied by either 
or by 1 electron, for instance due to strong correla- 
tion effects in the id shell. As there is only one localized 
core-hole present in the intermediate state. He then has 
the interesting property H]~, = U^ ^Hg for any integer 



I > (Van den Brink and van Veenendaal 2005 20061, 



which implies that He is either or Uc- This directly 
implies the relation HcGc = Hc{zk — Uc)~^- One now 
obtains for the indirect RIXS amplitude 



pnMrect ^ 



He 



Zk - Uc 



a + HoG)V\g), (51) 



Note that this expression is exact, but of course specific 
for the present form of the core-hole potential; gener- 



alized forms are given in (j 


\ment et al. 


Brink and van Veenendaal 


2005| 2006 



2007 Van den 



which include 

the spin and possible orbital degrees of freedom of the 3c? 
electrons. 

In the leading order of the UCL expansion one retains 



in Eq. (51 1 the first order term in He so that 



■jrindirect (f | HcV |g) , \lf\^d\x 

^ {z^-^){z^~Uc) = ^^"^ '^^ ' 

(52) 

where the resonance function 

P(^k,^k') = (PiV4p)'t^c((zk-w)(zk-C/c))"' (53) 

was introduced, where Pis,Ap is the Is — )■ 4p dipole tran- 
sition amplitude. The generic shape of the resonance 
function depends on the form of the core-hole potential. 
It is remarkable that the RIXS amplitude found in lead- 
ing order of the strong coupling UCL is directly related 



to the dynamic structure factor S'q(w) of Eq. (49), which 
is a situation very similar to the weak coupling pertur- 
bative approach, see Eq. (|48|. In fact for t/c the 



strong coupling UCL resonance function reduces to the 
perturbative one. This result has important implications 



for the interpretation of RIXS spectra since this approach 
then suggests that with proper handling of the prefactor, 
RIXS can be considered as a weak probe that measures 

S'q(cj). 

The sub-leading contributions to the indirect UCL 
scattering amplitude of Eq. (51) are of the type 
HcHqHc- Such terms a priori cannot be reduced to 
a response of Pq because Hq and He do not commute. 
Physically this term corresponds to an electron (or hole) 
hopping onto the core-hole site in the intermediate state. 
Denoting the hopping amplitude as t, these contribu- 
tions to the scattering amplitude is down by a factor 
t/{zk — Uc) with respect to the leading term. When go- 
ing off-resonance corrections to the UCL expansion thus 
become progressively smaller. On resonance these terms 
constitute contributions to the RIXS intensity of the or- 
der of {t/r)'^, which are thus governed by Uc and the 
inverse core-hole life-time F. Corrections to the UCL are 
thus smaller for shorter-lived core-holes. In cuprates, for 
instance the effective 3d valence bandwidth t « 0.4 eV 
and such corrections are expected to be moderate. Using 
for a specific system the commutation relation Hq and 
He such a higher order term can be calculated explicitly 
and again be cast in the form of a product of a resonance 
function and a generalized charge response function. 

The observation that within the UCL the RIXS cross 
section can be factored into a resonant pre-factor and the 
dynamic structure factor, Sq{uj) was tested experimen- 



tally (Kim et al. 



Kim 



2009a| , see also Section |V.B.3 
et al. (2009a I reporting an empirical comparison of Cu 
K-edge indirect RIXS spectra, taken at the Brillouin- 
zone center, with optical dielectric loss functions mea- 
sured in a number of copper oxides: Bi2Cu04, CuGeOa, 
Sr2Cu304Cl2, La2Cu04, and Sr2Cu02Cl2. Analyz- 
ing both incident and scattered-photon resonances |Kim 



et al. (2009a) extracted an incident-energy- independent 
response function. The overall spectral features of the in- 
direct resonant inelastic x-ray scattering response func- 
tion were found to be in a reasonable agreement with 
the optical dielectric loss function over a wide energy 



range. In the case of Bi2Cu04 and CuGeOa (Kim et al. 



2009a) observe that the incident-energy- independent re- 
sponse function, Sq^Q{uj), matches very well with the di- 
electric loss function, — Im(l/e(w)) measured with spec- 
troscopic ellipsometry, suggesting that the local core-hole 
approximation treatment of the UCL works well in these 
more localized electron systems. Corner-sharing two- 
dimensional copper oxides exhibit more complex exci- 
tation features to those observed in the dielectric loss 
functions, likely related to non-local core-hole screening 
efi^ects. 

The UCL expansion describes the RIXS cross section 
in the limits of small and large core-hole potential, re- 
spectively. In the intermediate region, one has to resort 
to numerical calculations (Ahn et al. 2009). In the dy- 



namics structure factor, excitations are created via pq. 
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implying that electron and holes are excited in an equiv- 
alent fashion. When dynamical effects are strong in the 
intermediate state, this can change and an asymmetry in 



daal and Carra 1997) that the Coulomb interaction be- 



the excitation of electron and holes can occur ( Ahn et al. 



2009). Since the screening electron is strongly bound to 



the core hole in the intermediate state, it is more likely to 
be scattered to higher lying states. The hole excitations 
on the other hand can delocalize and have a tendency to 
be closer to the Fermi level. 

Besides charge excitations also magnetic and orbital 
excitations were studied with the UCL. Theoretically the 
two-magnon response of antiferromagnetic La2Cu04 was 



tween the Is core hole and the excited electron creates an 
intermediate-state exciton. The result is that, although 
the final states are equivalent, the spectral weight is re- 
distributed as a function of the incoming photon's energy. 
Scattering of electron-hole pairs by long-range Coulomb 
interactions leads to collective excitations known as plas- 
mons. Plasmonlike excitations were also observed early 



on in RIXS (Isaacs et al. 19961. Whereas the measure- 



ments of RIXS on wide-band insulators is important, and 



covered in more detail by Kotani and Shin ( jKotani and 
Shin 2001), RIXS experiments in recent years have fo- 



et al. 2008a), agreeing nicely with experiment (Hill et al. 



calculated within the UCL (Van den Brink 2007 Forte cused more on strongly-correlated electron systems 



2008). Collective orbital excitations were investigated 



YTiOa (Ament and KhaliuUin 20101 and compared to 



theoretically for LaMnOs (Forte et al. 2008b I and for B. Charge Transfer Excitations 



experiments on titanates (Ulrich et al. 2009). These in- 



direct RIXS investigations on magnetic and orbital ex- 
citations will be reviewed in Sections |V.E| and |V.D.3[ 
respectively. 



V. ELEMENTARY EXCITATIONS MEASURED BY RIXS 

In this Section, the elementary excitations, as out- 
lined briefly in the introduction, are discussed in detail. 
The Sections start with high-energy excitations, such as 
plasmons, interband and charge-transfer transitions, via 
intermediate-energy transitions, such as dd excitations, 
to low-energy excitations, such as orbitons, magnons, and 
phonons. 



A. Plasmons, Interband Excitations 

A significant portion of the early RIXS work con- 
cerned direct RIXS in wide-band insulators, such as Si 



1994 



(Rubensson et al. 1990), diamond (Johnson and Ma 



Ma et al. 1992), graphite (Carlisle et al. 



BN (Jia et al. 1996), B2O3 (O'Brien et al. 19931, hex 



aborides (Denlinger et al. 
[2003'), and GaN (|Strocov et al. 



1995), 



2002[ ) MgBa ^Zhang et~ar. 



2004 2005). In these 



systems, the RIXS process is very intuitive, as described 
in these papers and more conceptually in Refs. (Ma 



1994 Ma and Blume, 1995). The absorption step ex- 



cites an electron into the empty conduction band and 
the deexcitation removes an electron from the occupied 
valence band. The final excitation is an electron-hole 
pair created across the gap. One of the important con- 
clusions of these experiments was to show the feasibility 
of momentum resolution. Generally, these RIXS experi- 
ments are rather well understood within an independent- 
particle framework. However, on closer inspection, the 
presence of intermediate-state relaxation can be ob- 
served. In graphite, it was demonstrated (Van Veenen- 



To date, charge transfer excitations have been studied 
at metal K- L- and M- edges and oxygen K-edges. For 
the first three cases, a core-hole is created on the metal 
site, and charge moves in from the ligand site to screen 
the core hole, which then decays. Following the decay, 
there is a finite probability that some rearrangement of 
the charge remains - creating a charge transfer excita- 
tion. Thus, this is an indirect RIXS process. This case is 
illustrated in Fig. [18] drawn for the Cu K-edge. Similar 
diagrams may be drawn for the L and M -edges, where 
the core hole is respectively in the 2p or 3p shell. At the 
oxygen K-edge, the core hole is created on the oxygen 
site and the photoelectron annihilates an oxygen 2p hole, 
which is strongly hybridized with a metal d orbital. Then 
a different O 2p electron decays, filling the O Is core hole 
and leaving the system in an excited state. Thus this is 
a direct RIXS process. 

In the study of charge transfer excitations, RIXS of- 
fers the ability to measure their momentum dependence 
and to carry out studies under extreme conditions, such 



as pressure (Rueff and Shukla 2010). These advantages, 



combined with the fact that the energy scales are well 
matched to the achievable experimental resolutions mean 
that a large amount of experimental work has been per- 
formed utilizing RIXS. Here we review that work, focus- 
ing on what has been learned about the character of these 
excitations in the various compounds studied to date. 



1. Early work 

We first briefly mention the early work in this field, de- 
fined here as work carried out before 2000, and typically 
with lower resolution than the later work, to be discussed 
in more detail below. 



a. Metal K-edge The first work on charge transfer exci- 
tations at a metal K-edge was the seminal work of Kao 
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FIG. 18 Indirect RIXS scattering process for observing charge 
transfer excitations at a metal isT-edge, illustrated for the case 
of a copper oxide . Figure from ( Hill et al. 1998 1 
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and co-workers (Kao et al. 19961) who studied them in 
NiO. This was notable for being the first K-edge RIXS 
experiment and was largely responsible for opening up 
the field to higher energy x-rays. Platzman and Isaacs 
followed this work with a study of NiSii.sSeo.s (Platzman 



and Isaacs 1998 ) 



Following this work, attention immediately switched to 
the cuprates and there were a series of studies of increas- 
ing sophistication on a variety of cuprates. The first of 
these was work on Nd2Cu04 (Hill et al. 1998") in which 
the 6 eV antibonding charge transfer excitation was ob- 
served and explained with a simple Anderson impurity 
model. Abbamonte et al. then carried out work on 



This higher-resolution work revealed the charge transfer 
gap at 2eV and hinted at its q-dependence. 

The 6 eV feature is a bonding-antibonding transi- 
tion that can be understood within the limit of a small 
cluster. For a transition-metal ion surrounded by lig- 
ands, the ground state for n electrons (with n = 9 for 
divalent copper compounds) is given approximately by 
a|d") -I- /Jjd^+^L), where L stands for a hole on the lig- 
and. For positive charge-transfer energies, one has a > j3. 
In the presence of a core hole c, the configuration is mod- 
ified to a'|cd"+^L) -I- /3'|c(i"), when considering a partic- 
ular spin-orbit split edge. After emission of a photon, 
the system can return to the ground state. However, 
there is also a finite probability that the antibonding 
state /3|c?") — a|c?"+^L) is reached. The energy differ- 
ence between these two states is about 6-8 eV. This is 
called a local screening process since the screening charge 
comes from the ligands surrounding the site with the core 



hole ( 


Van Veenendaal et al. 


1993 


Van Veenendaal and 


Sawatzky 


1993 





However, there is an alternative longer-range screening 
process, which is often more effective, in particular in 
cuprates that have a small charge-transfer energy. In 



La2Cu04 and Sr2Cu02Cl2 (Abbamonte et al. 1999a). iRicel 11988) 



FIG. 19 Experimental Ni L3 RIXS spectra measured with 750 
meV combined resolution. The lower left panel shows the Ni 
L3 absorption spectrum measured with the same incident en- 
ergy resolution and the labels of the excitation energies used 
for the RIXS spectra. The polarization of the incident pho- 
tons was V-pol, i.e. perpendicular to the scattering plane. 



Figure taken from ( Ghiringhelli et al. 2005 I 



1993 Van Veenendaal and Sawatzky 1993), a screening 



electron comes from a ligand surrounding neighboring 
copper sites. The intermediate state is jcd""'"^) on the 
site with the core hole. An electron is then removed from 
a neighboring site. For cuprates, this creates a two-hole 



state known as a Zhang- Rice singlet (ZRS) (Zhang and 
In the final state 



one is then left, with a 
transition across the charge-transfer gap, which is about 
2-3 eV in cuprates. The filled 3^^° state and the ZRS 
can move independently giving rise to dispersive features. 
This latter fact was more fully explored in Ca2Cu02Cl2 



(Hasan et al. 2000) who observed strong dispersion of 
the 2 eV feature along the (tt, tt) direction and weaker 
dispersion along the (tt, 0) direction. This was modeled 
with a t,t ,t — U model calculated on a 4x4 cluster. 

Early works exploring the polarization dependence of 
the charge-transfer states included Hamalainen et al. 



bamonte et al. 1999b) 



(Hamalainen et al. 2000) and Abbamonte et al. (Ab- 



this nonlocal screening process (Van Veenendaal et al. 



b. Metal L-edge Charge-transfer peaks are generally rel- 
atively weak at the L edges. The main absorption feature 
has a configuration a"|cd"+i) + /3"|cd"+2L). Note that 
the electron excited into the 3d shell directly screens the 
core hole and there is no need for screening electrons from 
neighboring transition-metal sites. Therefore, nonlocal 
screening effects are expected to be small and one ex- 
pects mainly transitions between bonding and antibond- 
ing states of the transition-metal site and its surrounding 
ligands. 
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FIG. 20 The RIXS and absorption spectra calculated with 
the Anderson impurity model at the same excitation energies 
used in figure [191 Colors and symbols: red dashed lines are V- 
pol, black solid lines are H-pol. The spectra were broadened 
with a 500 meV FWHM Guassian in the final state. Figure 
taken from ( [Ghiringhelli et al.\ |2005| l 



These bonding-antibonding charge-transfer features 
become more prominent by tuning the incoming photon 
energy to the weak satellite features above the main ab- 
sorption edge. These satellites have wavefunctions of the 
form a"\cd"^'^■^L)—|3"\cd^'^^) with a higher Ugand charac- 
ter and therefore the overlap with final states that have 
a higher ligand character is larger. The satellites were 



already observed early on in FeTi03 and KMn04 (Bu- 
torin et al. 1997) with approximately 1 eV resolution. 



Higher resolution (400 meV) experiments were done on 
Ti02 ([Harada et al\\2Qm\, see Fig. I29l The antibonding 



states are at 14 eV energy loss. The antibonding states 
are usually at higher energies for early transition-metal 
compounds reflecting the strong coupling between the lig- 
ands and the mainly empty 3(i states. Other prominent 
ligand features, often described as non-bonding states, 
are observed around 5-10 eV. For Ti02, the intermedi- 
ate states are given by a"\cd^j) + P"\cd^d^i L^,) , where 
the symmetry of the Zd orbitals is introduced. When 
the electron is removed from the orbital d^, the system 
returns again to the ground state symmetry. However, 
when the deexcitation removes the electron from the d 
orbital, one is left with a configuraton \cd^L^, 
that the hole on the ligands no longer couples to the Zd 
states, since the d^ electron has a different symmetry. 
Therefore, since there is no coupling to the d° state, the 



Note 



features between 5-10 eV are nonbonding states (Harada 



et al. 2000) 



Another model strongly-correlated system is NiO, 
which has a 3d^ ground state. Several experiments 
have been performed that show dd excitations and the 
presence of bonding-antibonding charge-transfer features 



around 6-8 eV, in particular above the main absorption 
edge, see Fig. 19 (Ghiringhelli et al. 2005 Ishii et al. 



2001 Magnuson et al. 2002b). These RIXS data were 



also modelled using a NiOg cluster. Note, the presence of 
polarization dependence, not only in the dd excitations, 
but also in the satellites. The good agreement between 
theory and experiment, see Fig. 20 demonstrates the 
current understanding of the RIXS process at transition- 
metal L edges is in good shape. 



c. Metal M-edge Amongst the very first RIXS work on 
charge transfer excitations, Butorin et al. studied Ce02 
and UO3 at the Ce and U M-edges (Butorin et al. 



1996b). Using the Anderson Impurity model, these au- 



thors showed that this technique was sensitive to charge 
transfer excitations and that reasonable parameters could 
reproduce the experimental results. 



The oxygen K edge has been used 
In 2002, Harada et al. used O K- 



d. Oxygen K-edge 
only sporadically, 
edge RIXS to study SrzCuOzClz ([Harada et al|[2002|. 



These authors used the polarization dependence of the 
observed scattering, together with theoretical expecta- 
tions for the polarization of the dd excitations to con- 
clude that the feature at 2 eV was a ZRS CT excita- 
tion. In addition, they observe double magnon features 



around 0.5 eV. Duda et al. (2006) measured the oxy- 



gen K edge of NiO. They observed the presence of fea- 
tures related to dd excitations (~1 eV), nonlocal screen- 
ing ( Van Veenendaal and Sawatzky 1993 ) (~4-5 eV) , and 
local bonding-antibonding transitions (above 5 eV). In 
addition, they claim the presence of double-singlet exci- 
tations, where the two holes are transferred between two 
neighboring nickel sites creating two singlet excitations, 
which are about 2Jh, with Jh the Hund's exchange split- 
ting, higher in energy than the ground state that con- 
sists of two triplet states. The possibility of measuring 
momentum dependence at the Oxygen K edge has also 



been suggested (Okada and Kotani 2007). Learmonth 



et al. { 2009 ) used oxygen K-edge RIXS to probe the elec- 
tronic structure of the band insulator a-MoOs and ob- 
serve that when the incident energy is resonantly tuned 
to the TT* states at the Fermi level, momentum selectivity 
allows the observation of states only near certain sym- 
metry points reflecting the highly anisotropic electronic 
structure of a-MoOa. 



2. Momentum Dependence 

a. OD Cuprates Materials in which the CUO4 plaquettes 
are isolated from each other by virtue of the crystal struc- 
ture have been termed zero dimensional. Such systems 
are interesting since they represent perhaps the simplest 
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incarnation of the transition metal oxide physics - and 
one that is particularly amenable to cluster model-type 
calculations. In the context of RIXS experiments, they 
can also therefore provide insight into the resonant cross- 
section. 



Kim et al., studied Bi2Cu04 at the Cu K-edge (Kim 



et al. 2009a). In this paper, the authors analyzed their 



data by dividing out a resonant prefactor, following the 
prescription of Ref. (Ament et al. 2007) with ujres de- 



fined by the absorption spectrum and T and \U\, the 
core-hole lifetime and core-hole potential, respectively, 
viewed as fitting parameters. In practice, such a scaling 
approach is very similar to that first taken by Abbamonte 



\U\=1 eV. dLu et al. 
r = 3.2 eV and \U\ 



et al. (Abbamonte et al. 1999a I, and by Doring et al. 



( [Poring et al.\ |2004[) with U\=0, or Lu et al. who took 
2006). Kim et al. found that for 
the data taken at differ- 



=5.7 eV, 

ent incident energies for Bi2Cu04 collapsed onto a single 
curve, and further that this curve closely resembled the 
loss function as independently obtained from ellipsom- 
etry measurements on the same sample. This demon- 
strates empirically, the utility of the UCL approach, 
discussed above, at least for zero dimensional systems. 
In the same paper, these authors look at a number of 
other cuprates. They find similar agreement between 
the underlying RIXS response (i.e. without the resonant 
prefactor) and the loss function in most cases, though 
this agreement is better in the ID edge-sharing cuprate 
CuGeOa, than in two-dimensional cuprates. Such edge 
sharing cuprates are relatively localized systems because 
of the very small overlap of the oxygen orbitals in neigh- 
boring plaquettes as a result of the 90° bonds and this 
suggests that the UCL works best in localized systems. 
This is to be expected because of the local approximation 
to the core hole in the UCL. For the more two dimen- 
sional compounds studied in the same paper (La2Cu04, 
Sr2Cu304Cl2 and Sr2Cu02Cl2) the agreement was less 
good and additional features were seen in the reduced 
RIXS spectrum that were not present in the loss func- 
tion, suggesting that non-local features in the interme- 
diate state were important in determining the lineshape 
and that were not captured in the UCL. 

Hancock et al. studied a different zero-dimensional 
cuprate, CUB2O4 (IHancock et al.\ 120101 120091). These 



authors observed several features in the Cu K-edge RIXS 
spectrum for this material, specifically peaks at 1.9 eV, 
6.8 eV and 10.9 eV. They observed no momentum depen- 
dence in any of these features, consistent with the zero 
dimensionality of the system. They argued that the 1.9 
eV feature, which is resolution limited in energy, was not 
a charge-transfer excitation, but rather a. d ~ d excita- 
tion. The higher dimensional cuprates do see a charge 
transfer excitation in this same energy range (~2 eV), 
but in those cases it is attributed to a non-local CT ex- 
citation (i.e. where the ligand hole is shared between 
neighboring plaquettes and is not on the same site as 



the filled d shell). Such an excitation cannot exist in 
these zero dimensional cuprates, and only the Molecu- 
lar Orbital (MO) antibonding excitation is seen - where 
the ligand hole and filled Cu d shell are on the same 
plaquette. This occurs at 6.8 eV in this system. The au- 
thors suggest that the dd excitation channel could arise 
from symmetry breaking interaction with the 4p photo- 
electron, or the participation of phonons. They further 
note that the lineshape (Gaussian) and incident energy 
dependence of the energy loss and the width of the MO 
excitation are difficult to understand in the context of an 



electron only model for the RIXS cross-section ( Hancock 



et al. 2010 2009). They argue that the origin of these 



features lies in electron-phonon coupling during and af- 
ter the scattering process and construct a simple Frank- 
Condon model; as the charge is rearranged in the inter- 
mediate and final states, the potential landscape for the 
lattice is changed and the lattice is abruptly out of equi- 
librium, and thus contains many phonons. Their model 
reproduces the salient features qualitatively, though not 
quantitatively. They find that the coupling shifts the 
MO excitation energy upwards, relative to the bare en- 
ergy, by an amount hujphjg^f, where LOph is a typical 
phonon energy (taken to be 75 nieV in this work), and 
7ff- / is a m easure of the electron-phonon coupling ( Han- 



cock et al. 2010). They further argue that this shift may 



explain the anomalous dependence of the MO energy on 



the Cu-0 bond length observed by Kim al. (Kim et al. 



2004c I 



b. ID Cuprates The one dimensional cuprates have at- 
tracted attention again because of their relationship to 
the two dimensional superconducting cuprates, but also 
because interesting new behavior is expected in one 
dimension and because theoretical approaches tend to 
be more reliable in one dimension where, for example, 
larger clusters may be calculated. The one dimensional 
cuprates studied to date may be classified according to 
the connectivity of the CUO4 plaquettes. Either these 
are connected along their edges, so-called edge-sharing 
networks, in which neighboring plaquettes share two oxy- 
gens, or at their corners (neighboring plaquettes share a 
single oxygen), so-called corner-sharing networks. 

1. Edge-sharing ID cuprates The two prototypi- 
cal edge-sharing ID cuprates are Li2Cu02 and CuGeOa. 
The charge transfer excitations in these systems are of 
interest because of the contrasts one may draw with 
the corner-sharing 2D networks of the high-Tc cuprates. 
In these latter systems, the ZRS excitation is the basic 
building block of many effective theories of the electronic 
structure. It is stabilized by the antiferromagnetic ex- 
change between the Cu and the O and, as discussed, 
can propagate to neighboring plaquettes. In the edge- 
sharing network the Cu-O-Cu bond angle is close to 90°. 
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According to the Goodenough-Kanamori- Anderson rules 
this should make the nearest-neighbor copper-copper ex- 
change slightly ferromagnetic. Thus the dynamics of this 
excitation should be quite different in this geometry. 

Li2Cu02 and CuGeOs differ slightly in their the Cu-0- 
Cu bond angle. It is 94° in LiaCuOs and 98° in CuGeOs. 
This is significant because higher order processes such as 
Cu-O-O-Cu hopping contribute antiferromagnetic terms 
to the magnetic exchange and this makes the net ex- 
change sensitive to slight variations in this angle. Both 
have been studied extensively using RIXS. Cu K-edge 



RIXS has been carried out on both Li2Cu02 (Kim et al. 



2004b 


Qian et al. 


2005) and CuGeOa (Hasan et al. 


2003 


IQian et al. 2005| Suga et aL[ |2005| |Vernay et al. 


2008 


] and for Li2Cu02, both K-eAge (|Bondino et al. 


2007 


1 


Duda et al.\ 


2000a Learmonth et al. 2007 land 



Cu L-edge RIXS (Learmonth et al. 2007 (has been per 



formed. The gross features are similar in all these studies 
across these two materials. Specifically, an excitation is 
seen around 1.7 (2.2) eV energy loss, another around 3.8 
(5.4) eV and a third at 6.5 (7.6) eV, where the numbers 
refer to CuGeOs (Li2Cu02). There is in fact a variety of 
interpretations of these features. 

The 2 eV feature is attributed to a dd 



20041 |Kim TTd] |2004b[ [Learmonth TTaT. 



Duda et al. 



2007[ |Suga| 



et al. 2005), while others attribute it to a d-intersite 



ZRS excitation (Qian et al. 2005 Vernay et al. 20081. 



There is general agreement that the feature near 3.8 eV in 
CuGeOa is due to a non-local charge transfer excitation 
with the hole weight on the neighboring plaquettes and 



that the 6.8 eV feature is a local CT excitation (Bondino 



et al. 


2007 


Suga et al. 


2005 



Vernay et al. 2008). In 



Li2Cu02, Kim et al. attribute the 5.5 eV feature to a 
local CT excitation on the basis of an agreement with 



a calculation by Weht and Pickett (Weht and Pickett 



1998 ) , an assignment agreed to by Vernay et al. on the 



basis of a CuaOg cluster calculation (Vernay et al. 2008). 
The feature at 7.6 eV in this material, which resonates 
at the poorly screened intermediate state is also believed 



to be a CT excitation (Kim et al. 2004b) although this 



has not been reproduced to date in a calculation. 

A number of authors looked at the momentum depen- 
dence of these various features with Cu iC-edge RIXS 



Kim et al. 


2004b 


Qian et al. 


2005 


Suga et al. 


2005 



with the momentum transfer along the chain direction. 
These studies observed very little, if any, dispersion - 
consistent with the expected isolated nature of the pla- 
quettes. Kim et al. did observe some small dispersion of 
the 5 eV and 7.5 eV feature in the direction perpendic- 
ular to the CuO chains in Li2Cu02 with a bandwidth of 
^100 meV, suggesting that there is significant interchain 



couphng (Kim et al. 2004b) 



2. Corner-sharing ID cuprates In corner- 
sharing ID cuprates, neighboring CUO4 plaquettes share 
a single oxygen and the Cu-O-Cu bond is 180°. As a 




FIG. 21 (color online). Pseudocolor plot of Cu Tf-edge RIXS 
intensity as a function of momentum and energy transfer 
(fiajk=8982 eV)as measured in the ID corner-sharing cuprate, 
SrCu02. The highly dispersive feature has been interpreted 
as a holon-anti-holon continuum. Figure taken from ( |Kim| 
et aZ.l|2004a|) 



result the hopping along the chain is relatively large and 
quite different behavior is expected, and observed, in 
the charge transfer excitations of such materials. Two 
materials have been commonly studied in this class of 
materials, SrCu02(Kim et al. 2004a Qian et al. 2005 



Wray et al.\ 



et al. 2006 



2007) and Sr2Cu3 (Hasan et al. 2002 Seo 



Suga et al. 2005). In the latter, there is a 



single CuO chain per unit cell, in the former there are 
two essentially independent such chains, connected to 
each other in an edge sharing geometry. 

In general, the results from the different experiments 
for these two materials are in broad agreement. A sinu- 
soidally dispersing feature is observed between ~2 and 4 
eV, with a maximum energy around reduced momentum 
transfer q=7r and a higher energy, non-dispersive feature 
at ~6 eV, see Fig. [21] This latter is attributed to a lo- 
cal charge transfer (bonding to anti-bonding) excitation. 
The interpretation of the dispersive feature, however, is 
more interesting. In ID, S = \ chains, the low energy ex- 
citations are expected to be in the form of spinous (which 
carry spin and no charge) and holons (which carry charge 
and no spin). 

Hasan et al., noting that the dispersion was larger than 
in the 2D cuprates, and on the basis of comparison with 
exact diagonalization calculations of the single band Hub- 



bard model (Tsutsui et al. 2000), identified the disper 



sion of the onset of the excitation as representing the 



dispersion of holons in this ID material (Hasan et al. 
2002 ) . Kim et al. noted that the continuum starts at a 
much lower energy and has a weaker dispersion than that 
of the peak (Kim et al. 2004a). They carried out dynam- 
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ical density matrix renormalization group calculations of 
the extended single band Hubbard model for a 128 site 
cluster. Parameters were obtained by fitting to optical 
conductivity data on the same material. They found t 
= 0.435 eV, V/t = 1.3, and U/t = 7.8. The results 
of these calculations were consistent with the observed 
RIXS spectra, implying that the dispersive peak repre- 
sents a continuum of spin-less charge ±e collective modes 
(holons and antiholons). In contrast to previous calcu- 
lations (Stephan and Pencl 119961 iTsutsui et all 120001), 



no bound holon-antiholon exciton was observed, even at 
reduced q = n. The onset dispersion is weaker and out 



of phase with the dispsersion of the peak (Kim et al. 



2004a) and this is consistent with low energy field the- 



ory predictions which say that it should follow the single 
spinon dispersion ^Jlsingj + A (Controzzi and Essler 
2^01. 



In contrast to these somewhat exotic pictures of the 
low energy excitations, Suga et al. interpret the same 
feature in Sr2Cu03 as an excitation from the ZRS to 
Upper Hubbard Band (UHB), with a dispersion that re- 
flects the dispersion of these two bands - as calculated in 
a Hartree Fock approximation, with an RPA treatment 



of correlations (Nomura and Igarashi, 2004 Suga et al. 



2005). With the parameter set; t 



x,dp 



y,dp 



-1.4 



eV, tpp — —0.7 eV and Udd = H eV, they were able to 
reproduce the energy position and dispersion of this fea- 
ture, though their calculations predict a sharpening in 
energy of the feature at reduced q=(2n -I- 1)tt not seen 
in the experiments. Other calculations on corner-sharing 
ID cuprates were performed by Okada and Kotani using 



exact diagonalization ( Okada and Kotani 2006 ) 



A third ID corner-sharing compound, namely the so- 
called "telephone number compound" Sri4Cu2404i, and 



variants, has also been studied (Higashiya et al. 2008 



Ishh et al. 2007 Wray et al. 20071. This material is 



"self-doped" that is the parent compound already has 
holes in it, with a nominal Cu^'^^"*" valence. This ma- 
terial has a complicated crystal structure that may be 
thought of as being made up of two corner-sharing chains 
(the "ladder") and an edge-sharing chain, each with dif- 
ferent periodicities (which are incommensurate with each 
other). This fact allows the identification of spectral fea- 
tures with the two structural components based on their 
momentum dependence. In general, all three works see 
very similar phenomenology. 

Wray et al. studied Sri4Cu2404i and compared the 
response to ID SrCu02 and 2D Nd2Cu04. In the 14-24- 
41 compound, they found a broad dispersive peak below 
about 4 eV, which they fit to two broadened Lorentzians, 
and a non-dispersive higher energy peak at 5.3 eV. The 
found the two low energy peaks to disperse with a period- 
icity of that of the ladder, exhibiting a maximum energy 
at ■Kiadder- They noted that this dispersion was inter- 
mediate between the ID SrCu02 and the 2D Nd2Cu04, 
suggesting that this implied a continuous evolution of the 



dynamics on varying the dimensionality of the corner- 
sharing network from ID to 2D. These authors modeled 
this dispersion based on a single band Hubbard model 
in which the dispersion of the quasi-particles is given by 
2tL cos{kcL) ±ij_ in the ladder and 2tch cos(kcch), where 
the subscripts refer to the ladder and chain and t± is 
the hopping across ladder rungs. This gives a continuum 
of excitations and the authors associated the observed 
peaks with the edges of this continuum where intensity 
is expected to become large. 

Higashiya et al. also found two excitations in 
Sri4Cu2404i with the same momentum dependence as 
observed by Wray et aL, one at around 4eV, which they 
treated as a single excitation and associated it with ex- 
citations from the ZRS band into the UHB on the chain, 
and one around 5.5 eV associated with the local charge 
transfer excitation which appears to reside on both the 



ladder and the chain ( Higashiya et al. 2008 ) . (With the 
incident photon energy set to the pre-edge, they also saw 
a feature at 2 eV which they associated with d-d transi- 
tions). They found that the 4 eV feature exhibited signif- 
icant dispersion, with a bandwidth of approximately 1 eV 
and a periodicity commensurate with the ladder, exhibit- 
ing a maximum energy at niadder, consistent with Wray et 
at. The 5.5 eV feature is dispersion-less. These authors 
also carried out a cluster calculation on a CugOig corner- 
sharing ID cluster, with electron correlation treated ex- 



actly (Okada and Kotani 1997, 2006). They found that 



in order to model the data, it was important to include 
a Cu 3d - O 2p intersite interaction, Upd. Without this 
interaction, the calculated energy losses were too low by 
about 1.5 eV (i.e. by Upd). With this addition, they 
were able to reproduce the incident energy dependence 
and momentum dependence observed in the data for the 
4 eV and 5.5 eV features. The full parameter set in their 
model was A = 3.5 eV, (dpa) = -1.5 eV, {ppa)-{ppTr) — 
0.65 eV, Udd = 9-0 eV, Upp = 4.0 e V, Upd = 1.5 eV and 
Q = 7.0 eV ( [Higashiya et al.\ \2008l 

Ishii et al. studied Sri4Cu2404i, La5Sr9Cu2404i, 
and Sr2.5Caii.5Cu2404i. The hole concentration in 
La5Sr9Cu2404i is very small in both the ladder and the 
chain, while for Sr2.5Caii.5Cu2404i the doping in the 
ladder is high enough that the material becomes a super- 



conductor under high pressure ( Ishii et al. 2007 ) . These 



authors also identify the excitation between 2-4 eV as 
being associated with the ladder structural unit. In this 
case, in addition to the momentum dependence, they also 
looked at Ca2Y2Cu50io which lacks the ladders, but is 
otherwise similar. The Cu i^T-edge RIXS spectrum on 
this latter compound lacked the 2-4 eV feature. They 
attribute this feature to excitation from the ZRB to the 



UHB of the ladder ( Ishii et al. 2007 ) . As discussed above, 
this feature is highly dispersive. Ishii et al. found that 
the peak position shifted to higher energies as the dop- 
ing increased, but the dispersion remained roughly the 
same. In addition, they looked at the dispersion along 
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the "rungs" of the ladder (i.e normal to the ladder direc- 
tion) as well as along the ladder itself, finding the gap 
increased to higher energies in this direction. They mod- 
eled their data with exact diagonalization calculations for 
a single band Hubbard model on a ladder. This model 



reproduced the dispersion in both directions ( Ishii et al. 



2007). They attribute the insensitivity of the dispersion 
to doping to the fact that in these ladder compounds the 
antiferromagnetic correlations are left largely unchanged 
by doping - in contrast to the 2D cuprates, for which the 
magnetic correlations are strongly dependent and this re- 
duction in AF correlation has been argued to reduce the 



dispersion of this excitation as seen by RIXS (Tsutsui 



et al. 2003 1 



c. 2D Cuprates In this section, we review the work on 
charge transfer excitations in the 2D cuprates. In doing 
so, we will focus on the features between 2 eV and 4 eV. 
The other prominent CT feature is that around 6 eV, first 



seen in Nd2Cu04 (Hill et al. 1998) and La2Cu04 (Abba 



monte et al.. 1999a I and subsequently in all systems fea- 
turing Cu04 plaquettes. As has already been discussed 
above, this feature is attributed to a local CT excitation 
within the plaquette from the bonding to antibonding or- 
bitals. It's dependence on the Cu-0 bond length is dis- 
cussed in ( Kim et al. 2004c I . It is dispersion- less and is 



correctly predicted by, for example, single ion Anderson 
models. As such, it is relatively uncontroversial. There 
has been much more discussion on the complex starting 
around 2 eV, coincident with the optical gap in these 



materials (Uchida et al. 1991). 



In this context, 2D planar cuprates studied to date 
with if -edge RIXS include La2_a;Sr2,Cu04 (|Collart et al. 



2006 Ellis et al. 


20081 |Kim et aZ.[ |2002| |Lu et al. 


2006 


Wakimoto et al. i 


>005), Nd2-:rCe^Cu04 (Hill et al. 


1 1998 


Ishii et al. 


2005 


b), Ca2Cu02Cl2 (Hasan et al. 20001 


Sr2Cu02Cl2 (Abbamonte et al. 1999a 


, YBa2Cu307-5 


(Ishii et al. 


2005c 


i), and HgBa2Cu04_|_5 


Lu et al. 


2005) 




(0,0) 



(71,0) 



(Jl,JI) 



(0,0) 



FIG. 22 (color online). Dispersion of CT excitations for op- 
timally doped Hgl201 (squares), compared with previous re- 



sults below 6 eV for Ca2Cu02Cl2 (triangles) (Hasan et al. 



20001 and La2Cu04 (circles) ( Kim et al.\ 2002 1 , and new data 
for La2Cu04 (crosses). Figure from Lu et al. ( |2005 I. 



sive features were also seen in the hole-doped cuprate, 



Ndi.85Ceo.i5Cu04 (Ishii et al. 2005b I, though one should 



note that here they are intraband transitions. The 2 eV 
feature is only seen around the origin. This fact, and the 
dispersion of the complex was explained by an exact di- 



agonalization calculation of a 4 x 4 cluster (Ishii et al. 



2005b) 



However, in later work, initially on an optimally doped 



superconductor HgBa2Cu04+5 (Lu et al. 2005), it was 
pointed out that as a function of incident photon energy, 
a multiplet of discrete peaks were observed at approxi- 
mately 2, 4, 5, and 6 eV and that each of these was more 
or less dispersion-less in position (with an upper bound 
of 0.5 eV dispersion for each feature), though the relative 
spectral weight varied with reduced momentum transfer 
in the 2D Brillouin zone. This is shown in figure [22l This 
multiplet of charge transfer excitations was then seen in 
La2Cu04 suggesting that such a picture is common to 



( [Harada et al.\ |2002[) 
Sr2,Cu04, ( Ghiringhelli 



At the same time, soft x-ray work has been done at the 2D cuprates (Lu et al. 2006). They argue against 
the O K-edge on Sr2Cu02Cl2 
and at the Cu L-edges on La2- 



et al. 2004 2007) and Sr2Cu02Cl2, Bi2Sr2CaCu2 08, and 



Nd2-:,Ce^Cu04 (Ghiringhelli et al. 2004) 



the bound exciton picture on the basis of the very small 
dispersion of the 2 eV feature. Still later work, at still 
better resolution (130 meV) mapped out this variation 
in La2Cu04 (Elhs et al. 2008). The authors speculate 



The picture of the charge transfer complex, and its 
momentum and incident energy dependence has slowly 
evolved as the energy resolution of the instruments has 
improved. Early experiments saw this as a single fea- 



that the lowest energy feature, around 2 eV may still be 
a bound exciton, but that its dispersion is significantly 



reduced from the predictions of Zhang and Ng ( )1998| ) be- 
cause coupling to the phonons gives it a larger effective 
ture dispersing (I Abbamonte et al.\ |1999al IHasan et al.\ mass ( Ellis et al. 2008 1 . Excitonic states in this energy 



2000). Subsequently, with an energy resolution of 300 
meV, it was seen as two features, one around 2 eV, the 
other around 4 eV, with different dispersions in-plane and 



region are also found in an exact diagonalization study, 
even if for small cluster sizes the intensity of these fea- 



no dispersion out of the CuO plane ( Kim et al. 2002 1 . 
The authors speculated that the lower energy feature was 
a highly-dispersive, almost-bound exciton predicted by 
Zhang and Ng (Zhang and Ng 19981. Highly disper- 



tures tends to be overestimated (Chen et al. 2010). 

As might be expected, the full range of theoretical op- 
tions discussed in Section [TV] have been applied to calcu- 
lating the RIXS cross-section in 2D cuprates. Amongst 
the first of these was exact diagonalization of small clus- 
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ters (jHamalainen et al. 


2000 


Hasan et al. 


2000 Hill 


et al. 


1998 


Ishii et al. 2005b 


Tsutsui et al.\ 


1999 


2003) 



These calculations were able to reproduce (and predict) 
the general features of the dispersive excitations, assign- 
ing them to transitions from the effective ZRS band to 
the UHB. 

In addition, in work led by Igarashi and co-workers 



( Nomura and Igarashi 


2005 


2006 


Takahashi et al. 


2008 


2009 


) use realistic band structures, include correlation ef- 



fects at the level of RPA and multiple scattering effects 
from the core-hole (not a large correction to the Born 
approximation) to calculate the interband transitions for 
a three band Hubbard model. They find the interband 
transitions, plus correlations, can explain, broadly, the 
features of the lower resolution data, i.e. the two peak 
dispersing structure in La2Cu04, without the need to 
appeal to bound excitons. The fine structure observed 
in later experiments is not seen in these calculations. 
A similar approach is taken by Markiewicz and Bansil 



(Markiewicz and Bansil 2006 Markiewicz et al. 2008) 



who calculated the RIXS response for both hole doped 
and electron doped cuprates. These authors found 5 
"RIXS bands" in the undoped case, from 2 eV to 6 eV 
and mapped out their dispersion. They found that for 
electron doping, the lowest RIXS band was shifted down 
to nearly zero energy at the T point, in agreement with 



the experimental results of ( Ishii et al. 2005b ) 



d. 3D Cuprates Doring et al. (2004) studied the charge 



transfer excitations in the 3D cuprate, CuO, in order 
to investigate their polarization and momentum depen- 
dence. In this material, the CUO4 plaquettes are close 



to rectangular and arranged in two chains ( Doring et al. 



2004). They found no momentum dependence to the 5.4 
eV feature, when either the direction or the magnitude 
of the momentum transfer was varied. This is consistent 
with their assignment of this feature to a strictly local 
charge transfer excitation (Doring et al. 2004). No lower 



energy CT excitations were observed, though the energy 
resolution was only 0.8 eV. Ghiringhelli et al. looked at 



this same compound using Cu L-edge RIXS ( Ghiringhelli 



et al. 2004 ) and saw weak CT excitations out to 5 eV. 



Kim et al. { 2009b ) looked at a different 3D copper com- 
pound, namely CU2O ("cuprite") . In contrast to most 
of the rest of the materials in this review, this material 
is a weakly correlated system and in fact is a direct gap 
semiconductor with well characterized exciton lines, of 
the Mott-Wannier type. The exciton shows up clearly in 
the K-edge RIXS data (taken with an energy resolution 
of 130 meV) and shows strong sinusoidal dispersion with 
reduced momentum transfer, along the T — R direction, 
with a bandwidth of ^0.5 eV. The authors note that if 
one assumes that in a weakly correlated system, RIXS 
measures the joint density of states, then this dispersion 
should reflect the lower bound of possible transitions and 



follow the dispersion of the valence band. This latter has 
been measured by ARPES in this material to be 0.6 eV, 
in good agreement with this interpretation of the RIXS 
result. This same material has also been looked at with 



Cu L-edge RIXS. (Hu et al. 20081. 



e. Nickel-based compounds 

1. NiO NiO is one of the prototypical transition 
metal oxides and has thus received a great deal of at- 
tention over the years. Although simple band theory 
predicts it to be a good conductor, it is in fact an anti- 
ferromagnetic insulator at room temperature with a rel- 
atively large band gap (determined to be about 4 eV 
by optical measurements). Following the early work dis- 
cussed above, this compound has been returned to re- 
peatedly for studies at the O K-edge ()Duda ei al.\ ]2006|, 
the Ni L-edge (Ghiringh elh et al.\ 2005 ; Ishii^ aZ. , 20(311 



Kotani et al.\ |2006j |Magnuson 


et al.^ |2002b( |Matsubara 


et al. 


20051 


and the Ni 7^-edge 


Huotari et al.\ 


2008 


Kao 


et al. 


1996 


Shukla et al. 2003 





At the Ni L-edges, one sees a number of features. With 
the incident photon energy tuned to the peak of the ab- 
sorption, features are seen at around 1.1, 1.6 and 3.0 eV 
energy loss. These are assigned to d — d excitations (in- 
cluding spin-flip processes and non-local pair singlet cre- 
ation) ([Ghiringhelh "oH |2005| [Ishii et aL[ [2001| |Mag- 



nnson et al. 2002b Matsubara et aL 2005 ), and are dis- 



cussed in section V.D If, however, the incident photon 
energy is tuned to features in the absorption spectrum 
corresponding to a d^L configuration, then the charge 
transfer excitations are emphasized, as shown in figure 
[T0| Figure [T9[ shows the evolution of the charge transfer 



complex as a function of incident energy (Ghiringhelli 



et al. 2005). There is broad consensus in understand- 
ing such data, based on a Single Ion Anderson Model 
(SIAM) type approach, that is one Ni site surrounded 
by 6 oxygens, as discussed above in the context of the 
cuprates. Specifically, the model includes parameters for 
the charge transfer energy. A, defined here as the split- 
ting between the 3d* and the 3d^L configurations in the 
absence of hybridization, and a finite bandwidth, W, for 
the oxygen 2p states (approximated with a number of 



evenly spaced discrete levels (Matsubara et al. 2005)). 
Spin-orbit parameters and Slater integrals are taken from 
Hartree-Fock calculations, scaled by an ad-hoc scaling of 
typically 0.85 and the full 3d multiplet is included. In 
contrast to the cuprates, three interaction configurations 
are included, viz: 3d^ 3d^L and 3d^°L^. The 3d wave- 
functions are allowed to expand (contract) due to the 
addition of one 3d electron (presence of the core- hole). 
The RIXS process is then calculated using the Kramers- 
Heisenberg formula. Using parameters that also work for 
photoemission and optical spectroscopy produces broad 



agreement with the experimental data, see figure 20 
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The limitation of the SIAM approach is that it is only 
able to describe strictly local processes. In interpreting 



their O if -edge data, Duda et al ( Duda et al. 2006 1 uti 



lized a NigOig cluster, which they found was the mini- 
mum size that reproduced the excitation at 4.5 eV, ar- 
guing that this was therefore a non-local charge transfer 
(NLCT) excitation, as it is in the cuprates. (They also 
found that the 1.9 eV feature is a non-local double singlet 



creation excitation, see section V.D I 



This cluster is too large to carry out a full exact di- 
agonalization calculation, but by disregarding the d — d 
excitations and the next-nearest neighbor Ni-to-Ni CT 
paths, and restricting the basis set to O 2px.y^z and 
Ni 3d3;2_j^2 3_22_r2 , then the calculation could be carried 
out. Parameters very similar to the Ni L-edge work 
above were obtained. In particular, it was found that 
A = 
A 



3.8ey and [/d<i=7. 2 eV puda et al. 
Z.^eV and Udd=T- 2 eV from the L-edge work 



2006) (c.f. 



( Ghiringhelli et al. 2005)). This consistency is reassur- 
ing and furthermore demonstrates that NiO is indeed in 
the charge transfer limit. 

A holy grail of the field of strongly correlated electron 
systems is to understand the development of the strong 
correlation as one tunes from the well-understood limit 
of good metals to that of the strongly correlated bad 
metals. A step towards this understanding can be taken 
by studying how NiO evolves out of the charge-transfer 
limit as pressure is applied, and the bandwidths driven 
up. The hard x-ray nature of the Ni K-edgc allows the 
use of diamond anvil pressure cells, and thus the pres- 
sure dependence of the charge transfer excitations may 
be studied with this technique. Shukla et al. carried 



out such studies on NiO powder (Shukla et al. 20031. 



Working at room temperature and utilizing a horizontal 
scattering geometry and a scattering angle of 26 = 90° 
to minimize the elastic scattering from the sample and 
from the diamond, these workers were able to observe 
the pressure dependence of the charge transfer complex 
up to 100 GPa, figure [23] They found that the position 
of the features (for example at 5.3 V and 8.5 eV), shifted 
very little as a function of pressure, but that the peaks 
became smeared out as the bands broadened, suggest- 
ing that there would ultimately be a bandwidth-driven 
metal-insulator transition in this material. 



(Takahashi et al. 2007) calculated the K-eAge RIXS 



spectra for NiO, treating the correlations in NiO in 
an Hartree-Fock approach on a multi-orbital tight- 
binding model, augmented with dynamic correlation ef- 
fects within the Random Phase Approximation. The 
core-hole - valence interaction is treated perturbatively, 
within the Born approximation, which is controlled by 
the very large gap of NiO. Two prominent peaks are 
found around 4-6 and 8 eV in the continuum spectra 
which originate from the 3d eg electron densities of states. 
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FIG. 23 Charge transfer excitations in NiO as a function of 
pressure. Data were taken in a powder sample at q = 6 
with the incident energy tuned to the pre-edge resonance. 
Bottom dashed line: nonresonant background. As the pres- 
sure is increased the resonant intensity decreases and the 5.3 
and 8.5 eV structures are seen to smear out due to the in- 
creasing electronic bandwidth. Figure taken from Shukla et 
al. ( [Shukla et dl\ pOOS] ) 



2. Nickelates A natural extension to the NiO work 
is to look at La2Ni04, which is the structural analog 
of the 2D cuprates discussed above. As revealed by Ni 



K-edge work (CoUart et al. 2006), the charge transfer 



excitations behave quite differently in the nickelate than 
in the cuprate. The charge transfer gap is seen around 4 
eV, in good agreement with optical measurements on this 
system. A shoulder is seen on the CT complex, which is 
interpreted as a low lying CT exciton at 3.6 eV. In con- 
trast to the case of the lowest lying CT exciton observed 
in La2Cu04, this exciton is strongly localized in the nick- 
elate, see figure |24j The authors explain this difference 
by noting that the lowest lying excitation in the cuprate 
can form a S=0 singlet bound exciton, involving two sites, 
both of which are spinless and can therefore propagate 
freely through lattice without disturbing the the S=l/2 
antiferromagnetic background. In contrast in the nick- 
elate system, while there can be a S=0 singlet exciton, 
I d^L) , this is on a single plaquette and cannot propagate 
through the S=l antiferromagnetic background without 
disturbing that background. The non-local CT excita- 
tion discussed above for NiO, |d^L; d^) may also form an 
excitonic state, but this is not spinless. Thus on general 
grounds, a different dispersion is expected, and in fact 
observed. 



f. Manganese- based compounds 

1. MnO Early work on MnO was carried out by Bu- 
torin et al. who saw both d — d and CT excitations 
at the Mn L2,3 edges with about 1 eV energy resolu- 
tion (Butorin et al. 1996a). These authors, however, 
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FIG. 25 (Color online) Comparison between the calculated 
Anderson model results and the experimen in the charge 
transfer excitation region as measured by L-edge RIXS in 
MnO. Labels A-C refer to incident energies in the vicinity 
of the peak of the L3 resonance. H is the peak of the L2 



resonance. Figure taken from Ref ( Ghiringhelli et al. 2006a I 



FIG. 24 (a) RIXS dispersion in La2Ni04 along the 
[100] = [010] direction. The spectra have been shifted verti- 
cally, the baseline is indicated for each spectrum. Solid line: 
three-point smoothed guide for the eye; dotted line: BZ edge 
spectrum; open circles label crystal field excitations, (b) Dif- 
ferences between smoothed spectra at various points in the 
BZ and the BZ edge spectrum from (a), (c) Dispersion of the 
exciton. Figure taken from CoUart et al. (Collart et al. 20061 



did not attempt to model the charge transfer excitations. 
Ghiringhelli et al. continued the study of MnO with a 



much improved resolution of 320 meV ( Ghiringhelli et al. 
[2006a| ). With this resolution, the authors found a wealth 
of features in the inelastic spectrum. In particular, a 
number of sharp excitations were observed out to 6 eV, 
attributed to local d — d excitations, and a broad excita- 
tion between 7 eV and 14 eV attributed to charge trans- 
fer excitations. The authors modeled the spectra using 
the SIAM, and a local crystal field model. Both repro- 
duced the d—d excitations well, but as expected, only the 
SIAM approach reproduced the charge transfer excita- 



tions, as shown in figure 25 The parameters used in the 
SIAM model were the charge-transfer energy, A=6.5 eV, 
the on-site Coulomb interaction, Udd='!-'^ eV, the intra- 
atomic core-hole potential, C/dc=8.0 eV, the hybridiza- 
tion strengths, V{eg) — —2V{t2g)~l.^ eV, the crystal 
field splitting IQDq^Q.h eV and the O 2p bandwidth, 
W=3.0 eV. This latter was modeled as 6 discrete states. 
Spin orbit coupling and Slater integrals were calculated 
using atomic Hartree-Fock codes with the latter then be- 
ing scaled by 80%. The results were then broadened to 



account for lifetime effects, experimental resolution and 
the discrete O 2p states. As is shown in the figure, in gen- 
eral the agreement is reasonably good, reproducing the 
position, width, intensity and polarization dependence. 



2. Manganites Agui et al. utilized Mn L-edge 
RIXS to study the charge transfer excitations in the 
doped, bilayered manganite, Lai.2Sri,sMn207, for which 



the Mn ions have a nominal valence of d^-^'^ (Agui et al. 



2005 ). Data were taken at room temperature with the in- 



cident and outgoing wavevectors were in the plane, such 
that the contributions from the 3(ia;2_y2 were emphasized. 
They observed CT features from 5 eV energy loss to 10 
eV energy loss with a tail extending even higher. The 5 
eV feature was interpreted as the Zd^L charge transfer ex- 
citations. These authors obtained good agreement with 
cluster calculations for the Mn'^"'" ion when the full niul- 
tiplet and JT distortion were taken into account. They 
found no need to include any contribution from Mn^+ 
sites to explain the spectra and concluded that the elec- 
tronic organization was one of d^^ + d**^, rather than 
^3.6+ 

Kuepper et al. studied the single layer parent man- 
ganite LaSrMn04, working at the Mn L2.3- and O K- 
edges with 0.6 eV energy resolution. At the Mn L-edge, 
they were able to observe local Mn dd excitations around 
2.2 eV and charge transfer excitations in the range 6- 
9 eV energy loss which were well reproduced by mul- 
tiplet charge transfer calculations, including the details 
of the fine structure. (Kuepper et al. 2006). At the 



O K-edge the energy loss spectra reflects the site spe- 
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cific local partial density of states (LPDOS). By tuning 
the incident energy through the oxygen if -edge, the con- 
tributions from the in-plane (01) and out-of-plane (02) 
inequivalent oxygens may be emphasized. While good 
agreement was seen existing LDA-I-U band structure cal- 



MnO, 



culations (Park 20011 for the out-of-plane oxygens, less 



good agreement was found for the in-plane oxygens, sug- 
gesting that further work is required to understand the 
electronic structure in these materials. 



3. Cubic Manganites Inami et al finami et al. 



2003[ ) carried out the first Jf-edge RIXS study of 
LaMnOs. They observed peaks at 2.5, 8, and 11 eV 
attributing the 2.5-eV peak to a transition from effective 
LHB to UHB across the Mott gap, while the 8- and 11-eV 
peaks were assigned to excitations from the O 2p bands 
to the empty Mn 3d and 4s/4p bands, respectively. These 
authors also measured the momentum and polarization 
dependence of the features. No significant dispersion was 
observed in any of the features, while the 2.5 eV feature 
showed a two-fold periodicity in intensity under rotations 



around the scattering vector (Inami ei aL 2003). Follow- 



ing an approach detailed in (Kondo et al. 2001), these 



authors calculated the inelastic response in the orbitally 
ordered state and found that it had the same weak dis- 
persion and characteristic polarization dependence of the 
2.5-eV peak. 

It should be noted that this picture is not uncontro- 



versial. Semba et al. (2008) also calculated the K-eAge 



RIXS spectrum for LaMnOa. Describing the RIXS inten- 
sity as a product of an incident-photon-dependent fac- 
tor and a density-density correlation function in the 3d 
states, they calculated the incident photon dependence 
from an ah initio band structure calculation of the 4p 
density of states. The correlation function was evalu- 
ated using a multiorbital tight-binding model calculated 
within the Hartree-Fock approximation. Correlation ef- 
fects were included via the random-phase approximation, 
which were shown to strongly modify the spectra. Rea- 
sonable agreement was then obtained with the extant 
experimental data and the authors concluded that the 
features arose from band-to-band transitions to screen 
the core-hole potential - rather than orbital excitations 
created by the off-diagonal part of the 3d-4p Coulomb in- 



teraction as proposed by Inami et al (Inami et al. 2003). 



This debate echoes a similar one concerning the origin 
of the resonant enhancement of elastic scattering at the 
Mn K-edge when looking at orbital wave- vectors. Here it 
was eventually concluded that the 3d-4p interaction was 



too small to give rise to the observed effects ( Shen et al. 



2006D . 

Other studies followed this first experimental work, no- 



tably Ishii et al. ( Ishii et al. 2004 ) who studied the doped 



systems La2_a;Sr2;Mn03 for x=0.2 and 0.4 and found that 
the gap was partially filled by doped holes, but that the 2 
eV feature remained. They also studied the temperature 
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FIG. 26 Temperature dependence of the K-eAge RIXS in 
tensity at fixed energy transfer of 2.25 eV in Lao.8Sro.2Mn03 
The arrow indicates the insulator to ferromagnetic metal tran 
sition temperature (To) (Figure taken from Ref (Ishii et al. 
20041 



dependence of these features and found that the 2 eV 
increased in energy at lower temperatures and that the 
temperature dependence is only observed along the (100) 
direction, and not the (110) direction - using Pbnm- The 
authors speculated that the origin for this anisotropy lay 
in the magnetic interactions and the relative alignment of 
spins in the different directions on crossing the insulator 
to ferromagnetic metal transition. Indeed the RIXS in- 
tensity tracks the ferromagnetic order parameter, figure 

m 

This early temperature dependence study foreshad- 
owed more detailed work by Grenier et al. who 
studied a range of cubic manganites ( Grenier et al.\ 
2005). Specifically, these authors studied the tempera- 



ture dependence of the inelastic scattering at (0,2.3,0) in 
Lao.875Sro.i25Mn03, Pro.6Cao.4Mn03, Lao.yCao.aMnOa 
and Ndo.sSrp.sMnOa. These were chosen because they 
exhibit a variety of different ground states, respectively: 
Ferromagnetic insulator, ferromagnetic metal, antiferro- 
magnetic insulator, and antiferrimagnetic insulator. In 
some cases, exhibiting more than one such phase as a 
function of temperature. What these authors found is 
that the integrated spectral weight from 2eV to 5 eV, 
tracked the magnetic correlations, increasing on enter- 
ing a ferromagnetic phase and decreasing on entering an 
antiferromagnetic phase, regardless of the metallicity (or 
otherwise) of that phase. 

In order to understand these results, the authors 
turned to LDA-I-U calculations of LaMnOa (Grenier 



et al. 2005 1 . They argued that the observed tempera- 



ture dependence arose because, in the small-q region, the 
low-energy features in the RIXS spectra are dominated 
by intersite d-d transitions between Mn atoms (mediated 
via the 0-2p). Ignoring improbable spin-flip excitations, 
such intersite transitions can occur only between neigh- 
bors with the same spin alignment, since no states exist 
at low-energy in the opposite spin channel. Thus the 
more pairs of ferromagnetically aligned neighbors in the 
system, the stronger the low-energy spectra will be, con- 
sistent with the experiments. 
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Working at the Mn L-edge, Butorin et al. showed 
that the AIM and full multiplet theory could explain 
the RIXS profiles in thin films of Lao.sCao.sMnOs as 
long as the calculations were performed for Mn^+ ions 

Related work in 



and not Mn'*+ ( [Butorin etd] |2002[ ). 
cubic manganites includes that in Lai 
Lai-a;Baa;Mni_yMy03, with TM=Co,Ni (|Bondino et al. 



cNaajMnOa and 



2004 Kuepper et al. 2000 1 




2 4 6 

Energy (eV) 



g. Other Compounds Charge transfer excitations have 
also been studied in a range of other oxides in recent 



years. Examples include CoO at the Co L-edges (Magnu- 



son et al. 2002b ) and L and M-edges ( Chiuzbaian et al. 
2008). Here as in the materials discussed above 
SIAM was found to describe the data well. In ref. ( Chi- 



the 



uzbaian et al. 2008 lit was further shown that the same 
set of parameters could describe both data sets providing 
reassurance as to their validity. Charge transfer excita- 
tions, at 6 eV, have also been studied in NaV205 at the 



V L-edges, see e.g. ref (Zhang and Callcott 2006 Zhang 



et al. 2002 ) . In this latter work, the angular dependence 
was also probed. It was found that the CT excitation 
had very little angular dependence, though the ratio of 
the CT to d — d excitations changed, with the d — d be- 
coming substantially larger as the incidence angle was 
increased. This was attributed to the change in the inci- 
dent polarization varying the excitation into the different 



d orbitals (Zhang and Callcott 2006). 



Ag ui et al. studied F eTiOg ("Ilmenite") at the Ti L 
edge ( |Agui et al\ |2009[ ). Here, Fe is divalent , d^ 
and Ti is tetravalent Ti^+ ,d°. Both are octahedrally 
coordinated and occupy neighboring sites in the basal 
plane. Importantly, the neighboring octahedra are face- 
sharing. In addition to the usual charge transfer features 
seen between 5-10 eV energy loss, new features were seen 



at 2.5 eV and 4.5 eV that were absent in Ti02 (Matsub- 



ara et al. 20001 - and also absent in single cluster (TiOe) 
calculations. These authors therefore carried out calcu- 
lations for a double cluster (FeTiOg), i.e. the central Ti, 
surrounded by its oxygen neighbors and a neigboring Fe 
and its ligands. Such a double cluster was able to re- 
produce the new features in the spectra and they were 
therefore attributed to intermetallic charge transfer ex- 
citations (i.e. between the Fe and the Ti). The authors 
state that these excitations are not mediated by the oxy- 
gen, but are rather "direct" intermetallic charge transfer 
excitations. 



3. Incident Energy Dependence 

By its very nature, RIXS is characterized by a very 
strong incident energy dependence to the inelastic re- 
sponse. Here, by incident energy dependence, we are not 
simply referring to which edge is being utilized, but also 



FIG. 27 The result of dividing the energy denominator, equa- 
tion |54| out from a series of incident energy spectra for 
the molecular orbital charge transfer excitation at 6 eV in 
La2Cu04. The thick line is a Gaussian fit ( [Abbamonte et al\ 
1999al). 



to the position within the edge. As discussed in section 
[n] the incident energy determines not only the overall 
size of the signal, but also the relative weight of the var- 
ious excitations that are observed. This makes it hard a 
priori to determine the response function of the, for ex- 
ample, 3d electrons from a series of measurements when 
that response function is apparently strongly dependent 
on the energy at which it was measured. 

Focusing here on indirect RIXS measured at the Cu K- 
edge, there have been a number of attempts to "divide 
out" the incident energy dependence and experimentally 
derive something that might represent a characteristic re- 
sponse function, see Section [IV.CI Fundamentally, these 
attempts are based on the assumption - motivated by 
various theoretical approaches - that the RIXS response 
may be factored into a resonant prefactor and a reduced 
response function. The resonant prefactor is generally 
taken to be a variant of 



{uJk' - UJres)^ + Ffc' 



(54) 



Where w^esi ^k.k' and Uc are potentially experimen- 
tally determined parameters representing the resonance 
energy, the core-hole lifetimes (in principle allowed to be 
different) and the core-hole potential, respectively. Ex- 
pressions similar to this have been derived a number of 



times from different perspectives (see e.g. (Abbamonte 


et al. 


1999a; 'Ament et al. 


2007| |Van den Brink and van 


Veenendaalj |2006f jDoring et al. 


2004 


1) 



The first attempt to recover the response function was 
due to Abbamonte et al. ( 1999a), who motivated such an 



expression (albeit with Uc=0 and F = Ffc = Ffe/) based 
on a third-order shake-up description of RIXS ( Abba- 
monte et al.\ |1999a| Platzman and Isaacs 1998). They 
then took ujres, F as fitting parameters, fitting the data 
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to find the best fit values. Dividing each individual spec- 
trum by the appropriate resonant prefactor then gave the 



results shown in figure 27 in which all the data (in this 
case for the molecular orbital charge transfer excitation 
in La2Cu04), approximately collapse onto a single curve. 
They argued that the resulting function was a projected 
response function, that is the valence part of the dynamic 
structure factor, S'q(cj), weighted by the form factor of 
the core state. 

Doring et al. took a similar approach, deriving equa- 
tion 54 (again with Uc = 0) again from third order per- 
turbation theory for shakeup processes (Doring et al. 
[2004). In looking at CuO they observed that the dis- 
persion of the shakeup features with incident energy re- 
quired that two intermediate state resonances were in- 
volved, with distinct excitation energies, labeled A and 
B. Interestingly, they found Tk' for resonance A 

at 8981.4 eV, but Tk = Tk' for resonance B, at 8995.8. 
With these parameters, they could simulate the incident 
energy dependence of their features, both in terms of in- 



tensity and energy position ( Doring et al. 2004 ) 



Lu et al. { 2006 ) followed the same approach to model 
their data in La2Cu04. They found 4 energy loss features 
at 2.3, 3, 4 and 7.2 eV, respectively and showed that their 
experimental data could be modeled with the incident 
energy dependence of Eq. (54). They found Tk = Tk' 



and Uc = 0, for the 3 lower excitations and Tk — 1.7 eV, 
Tk' = 1.5 eV and Uc = 1 eV for the high energy, 7.2 eV 
feature. Similar results were obtained with the incident 
polarization in the CuO plane or normal to it, apart from 



an overall 4 eV shift in the incident energies (Lu et al. 



2006) 



Following this, and with a view to testing the predic- 
tion that the underlying response function is closely re- 



lated to Sq{u}) ( 


Ament et al. 


2007 


Van den Brink and 


van Veenendaal 2006 


), Kim et al. carried out a similar 



procedure, dividing the resonant prefactor out for a range 
of cuprates in the hopes of obtaining the material depen- 
dent response function. A range of cuprates were stud- 
ied; Bi2Cu04, CuGeOa, Sr2Cu304Cl2, La2Cu04 and 
Sr2Cu02Cl2 (iKim et al.\ |2009a|). Taking the peak in 



the XAS spectrum for each material, ojxas as being the 
incident energy resonance (i.e. ujxas = ^res + \Uc\) ■, 
and forcing F = Ffe = Tk' , there are three fitting param- 
eters for each compound. These were chosen to optimize 
the data collapse. Values for F ranged from 2.5 - 3.5 
eV for these materials and Uc from 2.6 - 5.7 eV (Kim 



et al. 



in figure 28 



2009a I. The results of this procedure are shown 
The data collapse is reasonable, though 
not perfect. The collapsed data are then compared with 
measurements of the loss function from EELS and ellip- 
sometry experiments on the same materials. Again the 
agreement is reasonable, for example, the leading edge is 
more or less reproduced for all materials, though, it is 
not perfect. Notably, it is better for the zero dimensional 
Bi2Cu04 and edge-sharing CuGeOs cuprates than for 




FIG. 28 (Color online) Comparison of S(a;) with the mea- 
sured /m(— l/e(a;)) for various cuprates. DifTerent symbols 
are used to denote RIXS spectra obtained with different inci- 
dent energies. EELS data from the literature are also shown 
for comparison. These incident photon energies (in eV) are 
(a) 8997.5 (■), 8998.5 (O), 8999.5 (A); (b) 8988 (■), 8990 
(O), 8992 (A), 8994 (□); (c) 8992 (■), 8993 (O), 8994 (□); 
(d) 8991 (■), 8992 (Q), 8993 (A), 8994 {£>)■ and (e) 8993 
(■), 8994 (O), 8995 (□). The EELS data for CuGeQ3, 
La2CuQ4, a n d Sr2Cu02C12 are taken from Refs. (|Atzkern 
et al.\ [2MT| IMoskvin eFoTI |2002l [Terauchi and Tanakaf 
19991, respectively. Although there is no EELS spectrum for 



Sr2Cu304Cl2 available in the literature, that of Ba2Cu304Cl2 



is shown ( Knauff et al. 1996 1 . Figure from Kim et al. ( 2009a I . 



corner-sharing geometries. This suggests that the ultra- 



short core-hole lifetime expansion (Ament et al. 2007 



Van den Brink and van Veenendaal 20061 works better 



in more localized systems. 
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C. Charge-density waves 



books ( Griffith 1961 Sugano et al. 1970). Let us in 



Recent work by Wakimoto et al. (Wakimoto et al. 



2009) has indicated the possibility of RIXS to pro- 
vide information on charge-density waves. Charge ex- 
citations were studied in stripe-ordered 214 compounds 
La5/3Sri/3Ni04 and 1/8-doped La2--K(Ba or Sr)^; CUO4. 
They studied the momentum-dependence of RIXS at the 
Ni and Cu K edges. At a momentum transfer corre- 
sponding to the charge stripe spatial period, an increase 
in inelastic spectral weight around transferred energies 
of 1 eV is observed for both the diagonal (nickelate) and 
parallel (cuprates) stripes. These features, which are be- 
low the regular charge-transfer gap, are interpreted as 
collective stripe excitations or anomalous softening of the 
charge excitonic modes of the in-gap states. These modes 
are related to fluctuations of the charged stripes, i.e. col- 
lective charge stripe motion in the form of meandering 
and compression modes. Another plausible explanation 
is that the observed anomaly is due to the in-gap state of 
the striped system, in particular, because the excitonic 
mode of the in-gap state may have anomalous softening. 



D. Crystal field and Orbital Excitations; Orbitons 

The RIXS spectra show the presence of local dd transi- 
tions at the K, L, and M-edges. In optical spectroscopy 
local dd transitions are forbidden because they are not 
dipole allowed. Due to an additional coupling to the 
lattice can merely become weakly allowed. However, the 
transition matrix element in L- and M-edge RIXS is very 
different since it is the product of the reduced matrix el- 
ements for the absorption and emission steps, which are 
both allowed. The main body of RIXS work on dd ex- 
citations is therefore done at the transition metal L and 
M-edges, to be surveyed in Section [V.D.1| 

At the transition-metal K edge weak dd excitations 
are observed in the pre-edge region governed by Is 3d 
quadrupolar transitions. The presence of dd excitations 
in the Is — Ap region in transition-metal compounds is 
less easily understood since it is not a priori clear how 
the intermediate state exciton exactly couples to dd ex- 
citations, see Section [V.D.2[ 

At the L and Af-edges the dd excitations are well un- 
derstood. In general they can be described in a relatively 
local approach consisting of a single ion in a crystal field. 



or an ion hybridizing with its surrounding ligands ( Grif 



fith 1961 Sugano et al. 1970). The observed spectral 



lines are related to crystal-field excitations, e.g., t2g — >■ Cg 
in transition-metal compounds in octahedral symmetry 
and multiplet terms due to higher-order effects in the 
Coulomb interaction. In order to describe the dd excita- 
tions, a many-body approach is essential. 

A detailed description of Coulomb multiplet effects is 
beyond the scope of this review and is given in several 



stead consider the prototypical example of a system with 
just two spins. The underlying reason that independent- 
particle models fail to describe Coulomb multiplet is that 
they are single-Slater determinant approaches. 

First consider the spin part. In an independent- 
particle picture, two spins can only interact via 
SizSjz, which only makes a distinction between 
parallel (|tt) and and antiparallel (|ti) and lit)) 

spins. The energies are respectively. In a many- 

body framework, however, spins can interact via a 
Heisenberg exchange —J'^ij Si • Sj which has as eigen- 
states a triplet with S = 1 (|tt), 75 (Itl) + lit)): III)) 
and a singlet with 5 = (^(Iti) + lit)))i with energies 
— jJ and I J, respectively. Note that the 5 = and 1 
states both have an Ms — component, in which the 
spins are antiparallel. These terms are described as a 
linear combination of Slater determinants. 

In addition to the dependence on the spin, the 
Coulomb interaction has a strong orbital component, 
leading to the well-known LS coupling schemes to de- 
scribe the eigenstates of the Coulomb interaction. For 
example, an atom/ion in spherical symmetry with two 3d 
electrons or two 3d holes has five distinct multiplets: ^F, 
^P, ^G, ^D, and ^S. Hund's rule indicates that the is 
the lowest term. Within an independent-particle frame- 
work without orbital polarization effects, one is only able 
to describe two different energies: spins parallel and spin 
antiparallel. When lowering the symmetry, the differ- 
ent Coulomb multiplet split. For example, the ^F term 
branches into ^Ti, ^T2, and ^A2, where for two electrons 
the ^Ti{tlg) is the lowest m energy and for 3d^, the 
with two Bg holes is the lowest. 

The dd excitations are often rather well described 
with numerical codes using a single ion, small cluster. 



or Anderson-impurity type calculations (De Groot and 



Kotani 2008 Kotani and Shin 2001). Despite an of 



ten lower resolution, the advantage of RIXS over optical 
spectroscopy is that the dd transition are allowed and 
the relative intensities can be described in a straightfor- 
ward fashion. This allows a better identification of the 
dd excitations. 



1. dd excitations with direct RIXS 

The creation of dd excitations with direct RIXS is best 
understood in an ionic framework. The absorption of an 
x-ray photon starting from Sd" ground state produces a 
c3c?"+^ intermediate state, where c denotes a hole in a 
core level. The excited electron is strongly bound to the 
core hole, forming a local exciton. It therefore does not 
easily delocalize to neighboring sites. 

The subsequent deexcitation can produce an excited 
3d" configuration. RIXS at transition-metal L and M 
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FIG. 29 Polarization dependent soft-x-ray Raman scattering 
spectra of Ti02 plotted relative to the elastic peaks located 
at zero. Five vertical dotted lines between 5 eV and 10 eV 
indicate Raman scattering to the multiple splitting 3d^L~^ 
non-bonding states, whereas a higher energy shift line around 
14 eV indicates Raman scattering to the antibonding state. 
Solid bars on each spectrum indicate peak positions of the 
spectrum labeled n fluorescence ( [Harada et al.\ |2000| ) . 



edges has been measured for 3(f systems (Agui et al. 



2009 Butorin et al. 1997 Harada et al. 2000), early- 



transition-metal compounds (Braicovich et al. 



and Callcott 



2006 



ganese compounds (Agui et al. 2005 Ghiringhelli et al. 



Kurmaev eral] [20031 JUlrich et aL[ |2009[ |2008l 



Zhang et al.\ 2002) via man- 



2007a 



ZhangI 



2008, 2006b) to late-transition-metal systems, such as 
CoO Ichiuzbaian et al. ,2008.; Magnuson et aLl|2002b|, 



NiO (Chiuzbaian et al. 



2005 



Ghiringhelh et al.\ 2005 



Ishii et al. 


2001| Magnuson et al. 


2002a| Matsubaral 


etal. 2005 


and copper compounds ( 


Butorin et al 


.1995 


Duda et al. 


2000a Ghiringhelh et al.\ 2004 2007 


Kuiper 


et al. 1998 


). 



The RIXS spectral hne shapes often show a strong de- 
pendence on the polarization of the incoming and out- 
going X-rays. The angular dependence can be under- 
stood by looking at the angular dependence of the effec- 
tive transition operators for direct RIXS as described in 
Section IV. B Let us first consider the simplest effective 
operator with x — 0, which e qua l the num ber operator 
for holes, w^'^'^ — n^, see Eqs. (33 ) and (34 ). The angular 
dependence is given by e' 



e, see Eq. (15). In an atomic 



viously, this scattering is strong for early transition-metal 
compounds that have many empty 3d states. The depen- 
dence on the polarization vectors can be most clearly seen 
is a system such as Ti02 (3d°) (Harada et al. 2000), al- 
though it is still present for late transition-metal systems. 



such as the cuprates ( Ghiringhelli et al. 2004 ) 



In RIXS experiments, one often employs a (close to) 
90° scattering condition. When the outgoing polarization 
is not measured, the outgoing beam can be described by 
two polarization vectors e'^ and e^^ in and perpendicular 
to the scattering plane, respectively. When the incoming 
light is linearly polarized, the direction of the polarization 
vector can be rotated around the incoming beam. When 
the incoming polarization vector is in the scattering plane 
(TT-polarized) , it is perpendicular to both outgoing polar- 
ization vectors at this 90° scattering angle and e^* • e = 
for i = 1,2. This means that there is no contribution of 
the scattering operator v/lf^. 

This is clearly seen in for instance Ti02, where for tt- 
polarization, the elastic line is close to zero, see Fig. 29 
When the incoming polarization vector is perpendicular 
to the scattering plane (tr-polarized) , it is still perpen- 
dicular to but parallel to ej. Hence, one should get 



a contribution from 



,,ddQ 



and indeed for Ti02 a large 



elastic line is observed, see Fig. [291 f Harada et a/.l|2000|). 



A similar effect is observed in the cuprates where a 
decrease in the elastic line is observed when changing 
from (T to TT polarized incoming x-rays, see N1,V— > NI,H 
in Fig.[30|and = 90° ^- 0° and 6* = 0° in Fig. ^ Note, 
the ct/tt polarization condition is sometimes referred to in 



the literature as vertical/horizontal (Ghiringhelli et al. 



2004) and polarized/depolarized (Harada et al. 2000). 
The small elastic line in RIXS for 7r-polarized incoming 
light combined with the reduction of the zero-loss peak 
due to explain the popularity of the 90° 7r-polarized 
scattering condition. 

When employing the commonly-used 90° 7r-polarized 
scattering condition, the scattering operator w'^'^^ {x = 
0) is not relevant, and we have to consider the w^'^^ 



{x 



9 



1) operator. The angu lar dependence is given by 
|e-* X e, see Eq. (15). In the absence of a core- 
hole spin-orbit coupling the effective transition operator 



IS w 



ddl 



= L^, where is the angular momentum with 



II = 1,0,-1 = x,z,y. Let us clarify this with the ex- 
ample of a divalent Cu^+ ion in I?4ft, symmetry, which is 
typical for the cuprates. 

The ground state of Cu^+ in has a hole in the 
x^ — orbital (/i = 2). Let us take a 90° 7r-polarized 
scattering condition with the incoming wavevector along 
the z direction. In this geometry, we have e = x, e[ — y, 
and e'2 — z. Taking the outer product gives Ti^ ^ z, — y 
corresponding to ^ = 0, —1. For ^ = 0, the selection rules 
simplify to ^' — — /i, see Eq. (41), where the minus sign 



model, this term only gives rise to elastic scattering. Ob- 



is due to the fact that the orbital moment is an odd oper- 
ator. Since the ground state is = 2, the final state that 
can be reached is /i' = —2 or the xy orbital. For 63 = z. 
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FIG. 30 Upper panel: RIXS spectra of Sr2Cu02Cl2 measured 
at grazing and normal incidence with V and H polarizations 
([GI,V], [NI,V], [NI,H]). Bottom panel: calculated spectra for 
the same scattering geometries. Upper panel: RIXS spectra 
of Sr2Cu02Cl2 measured (iGhiringhelli et aZ. 1 120041). 
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FIG. 31 RIXS intensity for scattering from a single Cu^+ 
ion with initial state 3d^2_ 2 and spin up along the 2- 



ajcis ( Van Veenendaal 2006 1 



we need to consider ^ = — 1. Applying the selection rule 
gives /i' = 1 or the zx orbital. From Fig. 31 we see in- 
deed that for (p = 0° and 9 = 0°, the strongest intensity is 
observed for x'^ — ti yz/zx f- Additional inten- 

sity conies from spin-dependent and higher-order {Q = 2) 
terms. 

At this point it is useful to make a comparison with 
nonresonant inelastic x-ray scattering (IXS), where dd 



son et al. 2007). Whereas in RIXS the dd excitations 
are a result of two consecutive p • A processes, in non- 
resonant IXS, the local excitations are due to higher- 
order terms in the expansion of e"^'^ of the term. 
Finite scattering amplitudes are only obtained for even 
opera tors x = 0,2,4 (Van Veenendaal and Haverkort 
2008). Therefore, the scattering from x-^ — y'"'{fi = 2) ^ 

The transition is due 
4) which are generally 
2) transitions in non- 
When the trans- 



xy{^ = —2) requires |A^| = 4. 
to hexadecapolar transitions {x - 
weaker than the quadrupolar (a; 



2007) 



resonant IXS ( |Haverkort et al. 
ferred momentum is in the xy plane, odd ^ are zero due 
to the angular distribution. Scattering can occur with 
^ = 2. The selection rule in Eq. (39), then gives, for the 
cuprates, quadrupolar transitions from x'^—y^ — >■ 3z^ — 
(/i = 2 -J> 0) ([Van Veenendaal and Haverkort[ |2008|. 



This example clearly shows that, although nonresonant 
and resonant IXS probe the same excitations, their in- 
tensities can be strongly different due to the different 
selection rules involved. 

For undoped high- Tc cuprates, RIXS has clearly shown 
the feasibility of studying the variations in crystal field 



among the various superconductors ( 


Butorin et al.\ 1995 


Duda et al. |2000a| IGhiringhelh et al.[ 2004 


2007 


Kuiper| 


et al. 


1998). The strong polarization dependence of 



the spectral features aids in the peak assignment. The 
RIXS intensities are generally well described within a 
CUO4 cluster model. For doped cuprates, larger clusters 



are necessary ( Ghiringhelli et al. 2007 ) to properly de- 



scribe the effects of nonlocal screening (Van Veenendaal 



etal. 1993 Van Veenendaal and Sawatzky 1993). RIXS 



has also been consistently interpreted using small cluster 



models for TiOs (Harada et al. 2000) and FeTiOg (Agui 


et al. 


2009) 


NaV205 (Van Veenendaal and Fedro 2004 1, 


MnO 


Ghiringhelli et al. 2008,,2006b|), manganites (Agui 


et al. 


2005 


, CoO (Chiuzbaian et al. 


2008 Magnuson 


et al. 


20021 


3), NiO (Chiuzbaian et al. 


2005 De Groot 


et al. 


19981 IGhiringhelh et al. 2005| Huotari et al. 


2008 


Magnuson et al.\ 2002a| Matsubara et al. 20051. 



This clearly indicates a good understanding of the RIXS 
process at the L and M edges in transition-metal com- 
pounds. 

A very interesting feature of RIXS at the L and M 
edges is that the large core-hole spin-orbit coupling en- 
ables spin flip processes, since spin is not a good quantum 
number in the intermediate state. This is reflected i n th e 
presence of spin-dependent operators vZ^^ in Eq. (35 1. 



This opened the way for the measurement of magnons 



and their dispersion ( Ament et al. 2009 Braicovich et al. 



2010), which will be discussed in Section V.E.I 

Before this a crystal-field excitation combined with a 

In 



spin-flip was observed by Chiuzbaian et al 

7^ 



(20051 



in octahe- 



transitions are also observed ( Haverkort et al. 2007 Lar- 



NiO, starting from the Sdx^-y^t—s^^-r^t^' 
dral symmetry) with 5*2 = 1, transitions were observed 
to Sdf^ ^3dg .j. (with Ti or T2 symmetry) with Sz = 
flnal states, demonstrating the applicability of RIXS to 
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study magnetic excitations. Following |De Groot et al. 
(1998), the authors employ a crystal field model for a 



single Ni^+ ion and extract a superexchange constant 
6 J = 125 ± 15 meV. However, it is not entirely clear 
whether the superexchange for excited states is equiva- 
lent to the superexchange causing the antiferromagnetic 
order. In addition, exciting a combined orbital and spin 
excitation makes the observation of a single-magnon dis- 



persion impossible (Ament et al. 2009). 



2. dd excitations with indirect RIXS 

The RIXS cross section at the transition-metal K-edge 
is dominated by charge-transfer excitations due to the 
strong core-hole potential. However, dd excitations are 
also observed. In the pre-edge region, the dd excitations. 



as seen for example in NiO (Huotari et al. 20081 and 



the 4p-3ci Coulomb interaction gives an electron opera- 
tor instead of a hole operator. The angular dependence 
is again given by — [e'*,e]^. 

The dd excitations created change the symmetry of 
the valence shell. This change in symmetry has to be re- 
flected in a change in symmetry of the excited 4p electron. 
Therefore, the 4p electron is no longer a spectator as is 
often assumed in K-eAge RIXS. However, care must be 
taken with the reverse conclusion. A change in symmetry 
in the Ap state (as occurs in a 90° 7r-polarized scattering 
condition) does not directly imply that the symmetry of 
the valence band has changed in the scattering process. 
The symmetry relation are local point group arguments 
that do not hold throughout the Brillouin zone, or, in 
other words, in the wide 4p band the different 4p states 
are mixed allowing excitation in one particular orbital 
and de-excitation from another. 



cuprates ( Kim et al. 2004b ) are a result of direct RIXS 



Ishii et al. ( Ishii et al. 2010 ) measured the polariza 



due to quadrupolar transitions. 

The dd excitations observed when the incoming photon 
energy is tuned to the main absorption line (Is — > 4p), 
for example, in cuprates and NiO (Van Vcenendaal et al.\ 
2010 ) are more difficult to understand since direct tran- 
sitions into the 3d shell are very weak. Although, in 
principle, these excitations can be due to hybridization 



tion of the outgoing photons and performed a symmetry- 
based analysis of the dd excitations at the Cu K-eAge of 
KCuFs. The Cu ions are in a Sd^ configuration, where 
the hole has eg character in the ground state. Although a 
quantitative evaluation of the intensity requires a micro- 
scopic model for the scattering process, the authors point 
out that certain RIXS selection rules are in qualitative 
agreement with their data: t2g-eg excitations appear in 



with the 4p band ( [Hancock et d] |2010[ ), that coupling ^11 scattering configurations that they used, while e 



is forbidden in most symmetry directions and very weak 
elsewhere in the Brillouin zone (unless the inversion sym- 
metry of the crystal is broken). Additionally, it does not 



reproduce the right intensities for NiO (Van Veenendaal 



et al. 2010). The most likely origin of the dd excita- 
tions in the Is — 4p excitation region is the interaction 
between the excited 4p electron and the open id shell. 



Using the methodologies developed in Section IV. C 



effective transition operator can be obtained (Van Vee- 



nendaal et al. 2010[ ) 



= E^K,^k')c. {f\wf^ \g) (55) 



where the operators related to the Is and 4p level can 
be removed leading to an effective id transition op- 
erator w'''*^. The resonance function P(a;k,i^k') is of 
the form of Eq. (46 1 or (53). The coefficients Ck — 

npO 2^^1 a^^TTTr'l M^^n^ A p2 Lr*! 3_^3 

^-^ 15^ 35^ ' 5^ ""^ 35*-^ ' 35-^ 15^ 245^ 

for k = 0, 1, 2 and F'^, G^, and are expressed in 



excitations only occur in some. No dispersion is seen 
within the experimental resolution (600 meV). 



3. Orbitons 

Traditional dd excitations, the subject of the previous 
sections, arise because of a strong crystal field splitting of 
the d levels due to a local lowering of the symmetry of the 
lattice. In cases where crystal field splitting is small and 
superexchange dominates the interaction between orbital 
degrees of freedom, the orbital excitations have a com- 
pletely different character: they can in principle move 
through the system and thus have an appreciable dis- 
persion, in contrast to dd excitations. Such dispersing 
orbital excitations are often referred to as orbitons or 
orbital waves. The experimental observation of such dis- 



persive orbital modes is a hot topic (Polli et al. 2007 



Saitoh et al. 2001 Tokura and Nagaosa 2000 ) - for in 



terms of radial integrals of the id-4p Coulomb interac- stance the observation of orbitons in LaMnOs with Ra- 



tion (Griffith 1961) 



man scattering remains controversial ( Griininger et al. 



There is a remarkable similarity between the effective 
spin-independent operator from two consecutive dipolar 
transitions in direct L- and M-edge RIXS, see Eq. (33 1, 



2002 Saitoh et al. 2001 2002). RIXS is a promising 



and dd excitations created by the id-Ap Coulomb inter- 
action as observed with indirect RIXS. The resonance 
function, obviously, is different for direct and indirect 
RIXS and an additional constant Cx enters. However, the 
same one-particle tensor operator is obtained, although 



technique for measuring orbitons because it is both di- 
rectly sensitive to dipole forbidden dd excitations and has 
enough momentum to probe their dispersion in a large 
part of the Brillouin zone. 

Orbitons can occur in highly symmetric lattices, where 
the 3c? orbital degrees of freedom can remain active. The 
interaction between orbital degrees of freedom gives rise 
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FIG. 32 Schematic illustration of the orbital wave. Top, the 
electron cloud in the ground state. Bottom, a snapshot of the 



electron cloud when an orbital wave is excited ( ,Saitoh et al. 



to an orbital ordered ground state. There are two types 
of interactions between orbital degrees of freedom on dif- 
ferent sites, one related to the lattice and the other to su- 



perexchange ( Kugel and Khomskii 1982 ) . A cooperative 



Jahn- Teller lattice distortion removes the local orbital de- 
generacy and stabilizes a certain type of orbital ordering, 
concomitant with a lattice distortion of the same symme- 
try. Superexchange is well-known to generate magnetic 
interactions between spins. It relies on the presence of 
virtual hopping processes that, on 180° bonds, promote 
antiferromagnetic alignment of electron spins. In the vir- 
tual hopping processes, not only the electron's spins can 
be effectively interchanged, but also the orbitals that the 
electrons occupy. Superexchange can thus also drive or- 
bital ordering. The type of orbital order promoted by su- 
perexchange is closely related to the magnetic order, be 
it ferro- or antiferromagnetic, because the hopping am- 
plitudes depend both on the spins and orbitals involved, 
as stated by the Goodenough-Kanamori rules. If the su- 
perexchange interaction between orbitals dominates, the 
couplings to the lattice can be viewed as perturbations 



to the orbiton dynamics, see e.g. ( 


Van den Brink 


Van den Brink et al. 


2004 


Khaliullin 


2005 





Orbiton Scattering in RIXS There are three different 
mechanisms by which photons can scatter from orbitons 
in RIXS. The first and simplest one is by direct RIXS: 
the incident X-ray photon promotes a core electron to a 
certain valence orbital, and an electron from a different 
orbital on the same site fills the core-hole. This process 
is discussed in Section fV.D. II 

In indirect RIXS, there are two coupling mecha- 
nisms: the on-site shakeup mechanism (Ishihara and 



Maekawa 


2000 1 and inter-site superexchange bond mod- 


ulation ( 


Forte et al. 




2008b 




Ishihara and Maekawa 



2000). In the shakeup mechansim, the potential of the 
core-hole and photo-excited electron locally change the 
orbital of one of the valence electrons. 

In the superexchange bond-modulation mechanism the 
excitation of orbitons is due to the core-hole potential 



modifying the strength of the superexchange bonds in the 
intermediate state. This modification generates orbital 
flips on the two sites that form the bond. This scattering 



process was analyzed within the UCL expansion (Forte 



et al. 2008b), which in this case is well-controlled be- 
cause the energy associated with the modulation of su- 
perexchange bonds is much smaller than the inverse 
core-hole life-time F at TM iiT-edges. This photon- 
orbiton scattering mechansim is similar to the process 
by which bi-magnons are excited in indirect RIXS, see 
Section rVETbl 

The symmetry of the intermediate state potential de- 
termines which on-site orbital shakeup transitions are 
possible. It is argued by lAment and Khaliullinl (|2010|) 



that multiplet effects at TM X-edges, like relativistic core 
spin-orbit coupling, drive the rapid evolution of the inter- 
mediate states and thus tend to wash out any particular 
symmetry of the potential when the core-hole lifetime 
broadening F is much smaller than the multiplet energy 
scales. Transitions of different symmetries are all ex- 
pected to become equal and weak in strength, except for 
the ^i(g) component, which becomes dominant. Con- 
sequently the shakeup channel gives low-energy orbital 
spectra that are very similar for different L edges. 



a. Titanates. Orbital excitations were investigated at 
the Ti La-edge RIXS in the titanates LaTiOg and YTiOg 



(Ulrich et al. 20091. In these materials, the orbitally ac- 
tive Ti ions form a cubic lattice. Every Ti ion has a single 
electron in the 3d t2g subshell. The La-compound is an- 
tiferromagnetic, while YTiOs is a ferromagnet. Two sce- 
narios have been proposed for the orbital ground state of 



LaTiOa: a traditional static orbital ordered one ( Pavarini 



et al. 2005 ) and a fluctuating orbital state ( Khaliullin 



and Okamoto 2003). YTiOa has an orbital ordered 
ground state, where the order is being driven by either a 
cooperative Jahn- Teller distortion or superexchange in- 
teractions. 

At incident energies of 452.0 and 454.6 eV, correspond- 
ing to excitations into the t2g and eg levels, experiments 
have found a loss feature at 250 meV in both compounds. 



33 


(Ulrich et al. 


2009 


2008 



orbital excitations. Upon lowering the temperature well 
below the magnetic transition temperature, the orbital 
excitations gain approximately 50% in spectral weight 



(see Figs. 33 (a) and (b)), similar to the behavior of 
magnons in Raman spectra of the cuprates. This tem- 
perature dependence is hard to explain in the scenario of 
a lattice driven orbital ordering. In addition the spectral 
weight increases with increasing transferred momentum 
in YTiOa, while it decreases in LaTiOa. The authors 
note that this behavior is opposite to that of the magnetic 
excitations, which is natural for superexchange driven or- 
bital ordering. The q-dependence at the t2g edge is shown 



in Figs. 33 (c) and (d). The orbital excitation shows lit- 
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FIG. 33 RIXS spectra at the Ti L3 edge of LaTiOs 



YTiOs ((Ulrich et al. 



and 

2009 1). For the La compound, q is 
directed along the [110]-direction, while for the Y compound, 
q is along the [100]-direction. The incident photons were tt- 
polarized. Figs, (a) and (b) show the temperature dependence 
at scattering angle 26 = 90°, and Figs, (c) and (d) display the 
dependence on scattering angle 20. 



tie dispersion, ruling out that the 250 meV feature is a 
single orbiton. Comparison of the Cg edge data of YTiOa 
to theory gives a best description of the data by single- 
site coupling to fluctuating orbitals. However, the theory 
also predicts a single orbiton peak, which is not observed 
in the experiment. 

A detailed theoretical analysis of RIXS on YTiOs was 



given by Ament and Khaliullin (2010), who considered 



separately single site shake-up scattering and two-site 
bond-modulation processes. A situation in which single- 
site RIXS processes are considered in a superexchange 
driven orbital ordered ground state gives the best fit to 
experiments. The 250 meV feature is described by two- 
orbiton modes, which disperse very little. 

When the isotropic symmetry-component of the core- 
hole potential dominates over all others, the single-site 
RIXS processes are suppressed and two-site bond mod- 
ulation processes can become visible. It is found that 
the two-site channel yields two-orbiton excitations, repro- 
ducing the peaks at 250 meV and capturing the trend of 
increasing intensity with increasing q, albeit the q depen- 
dence is too strong compared to the experimental data. 
From a fit of their theory to Raman and RIXS spec- 



tra, the authors obtain an orbital superexchange constant 
J„b = 60- 80 meV in YTiOg. 

Calculations for a lattice distortion dominated scenario 
generate a two-peak feature at 250 meV that does not 
capture the intensity trend correctly. Further, the au- 
thors note that this scenario cannot explain both RIXS 
and Raman scattering spectra with the same crystal field 
parameters. 



b. Manganites. LaMnOa is the mother compound of the 
colossal magnetoresistance manganites. The Mn ions in 
LaMnOa are in a Sd** configuration and form a cubic lat- 
tice. Three electrons occupy the t2g orbitals and form 
an S = 3/2 spin. The fourth electron occupies one of 
the two Cg levels and is therefore orbitally active. A 
strong cooperative Jahn- Teller distortion stabilizes an al- 
ternating 3x^-r^/ 3y^-r^ orbital ordering pattern in the 
ab plane, where spins order ferromagnetically. Superex- 
change therefore also contributes to the orbital ordering. 
Optical Raman scattering features at around 150 meV, 



ascribed to orbitons by Saitoh et al. (2001), were inter- 
preted as multi-phonon excitations by [Griininger et al.\ 
(|2002|), see alsolSaitoh et aL|(|2002|). Whether RIXS can 



observe orbitons in this manganite was analyzed theoret- 
ically - low-energy RIXS data on LaMnOa has not been 
reported yet. 

Orbital RIXS at the Mn L and K edges were calculated 



by Ishihara and Maekawa (2000) for a model Hamilto- 



nian where orbital and spin degrees of freedom interact 
via superexchange. The L edge is treated in the fast 
collision approximation and only single orbitons appear 
in this case. In the A-type magnetically ordered ground 
state phase of LaMnOs , the orbiton modes have approxi- 
mately q-independent intensity. For indirect RIXS at the 



Mn K edge, Ishihara and Maekawa (2000) consider both 



possible excitation mechanisms. The shakeup process, 
mediated by the 3d-4p interaction, gives the same result 
as for the L edge, except for the polarization dependence. 
The inter-site process yields one- and two-orbiton exci- 
tations. The calculated single-orbiton intensity shows a 
strong q-dependence. 

The indirect RIXS spectrum of an orbitally active Cg 
system such as LaMnOa was calculated by jForte et aL\ 
(2008b) by considering a model Hamiltonian for the Cg 



orbital degrees of freedom interacting by superexchange. 
The authors calculate the bond modification scattering 
channel microscopically by applying the UCL expan- 
sion (see Section IV.C.3). The orbital RIXS spectrum. 



shown in Fig. |34[ shows both one- and two-orbiton exci- 
tation features. The one-orbiton spectrum is sharp with 
a strongly q-dependent intensity peaking at the orbital 
ordering vector. The two-orbiton contributions cause 
broader features at somewhat higher energy, the inten- 
sity of which is also strongly momentum dependent. At 
q = (7r,7r) and at the F-point, the intensity vanishes. 
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FIG. 34 K-edge RIXS spectrum of orbiton excitations in 
the A-type antiferromagnetic phase of a 2D eg system ( Forte 
et a/.||2008b|). 



E. Magnetic Excitations 

Measuring magnetic correlations in materials is tradi- 
tionally the domain of inelastic neutron scattering. With 
great precision it can probe lower energy magnon dis- 
persions and, more generally, the dynamic magnetic sus- 
ceptibility x(<3: ^) of spin-solids and liquids. It is a very 
powerful and mature technique, but it also has systematic 
difRcultics. As a rule it requires for instance large sam- 
ple volumes and experiments with materials that contain 
neutron absorbing elements (e.g. Cd or Gd) are difficult. 

As neutrons carry a spin a projected angular mo- 
mentum of ASz = or 1 can be transferred to the ma- 
terial in inelastic neutron scattering. The spin of the 
photon is often called polarization. Although the photon 
spin couples only very weakly to the spin of the electron, 
it does affect the orbital angular momentum of the elec- 
tron. The spin can be changed due to the presence of 
additional interactions, such as the spin-orbit coupling, 
which mixes the electron's orbital angular momentum 
and its spin. It was already observed in the early days of 



RIXS (Tsutsui et al. 1999j ) that as a photon has an an- 
gular momentum L = 1, it can transfer AL^ = 0, 1 or 2 
to the solid, potentially creating not only single magnon 
excitations with AS'^ = 1, but also two-magnon excita- 
tions with ASz = 2. Conservation of angular momentum 
in principle allows for such excitations to be created in 
photon scattering. 

In the past few years the resolving power of RIXS has 
rapidly increased, reaching combined energy resolutions 
better than 0.1 eV. This has made it feasible experi- 
mentally to probe elementary magnetic excitations with 
both direct and indirect RIXS in a number of transition 
metal compounds. A large effort has been devoted to 
benchmarking magnetic RIXS by using the technique on 
paradigmatic magnetic materials such as NiO, La2Cu04 
and Sri4Cu24 04i. These materials are close realizations 
of the elementary spin Heisenberg model in 3, 2 and 1 
dimensions, respectively. The successes of these studies 
helped to establish RIXS as an inelastic magnetic scatter- 



ing technique. This made it possible to take the next step 
and address with RIXS open physics questions on mag- 
netism in more complex materials, for instance in doped 



high temperature superconducting cuprates (Braicovich 



et al. 


2009 




2010 


Ellis et al.\ 


Schlappa et al. 


2009 


)■ 



From a theoretical perspective, however, it is impor- 
tant to note that transfer of angular momentum requires 
a microscopic coupling between the spin moment of the 
electrons in a solid and the orbital angular momentum 
introduced in the solid by a photon. For valence elec- 
trons the spin-orbit coupling is normally weak, but for 
electrons within the atomic core it can be very large, up 
to tens of eV's. In RIXS this strong spin-orbit coupHng of 
the intermediate state core-hole can generate a large cou- 
pling between the orbital moment induced by a photon 
and the valence electron's spin moment. In what follows 
we will review this spin excitation mechanism in detail. 

Magnetic excitations can also be probed in the AL^ — 
0, channel i.e. without transfering photon angular mo- 
mentum L to valence electron spin S. This is a well es- 
tablished process in optical Raman scattering, where the 
electromagnetic radiation couples to two-magnon excita- 
tions (an Sz — +1 plus —1 magnon combined to a total 
Sz = excitation). This coupling is due to the elec- 
tromagnetic field that modulates the strength of bonds 



between spins (Devereaux and Hackl 


2007 


Elliott and 


Loudon 


1963 


Fleury and Loudon 


1968 


I . Two-magnon 



scattering is also weakly allowed in optical conductivity 
experiments, where this ASz = magnetic mode can be 
excited together with a phonon. Being low energy pho- 
ton techniques, optical probes have the advantage of an 
outstanding energy resolution, but the disadvantage of 
being restricted to essentially hq — 0. 

The ALz = channel is also present in RIXS. In fact 
it is the dominant magnetic scattering channel if the in- 
termediate state core-hole is not subject to a spin-orbit 
interaction. This is the case at a transition metal or oxy- 
gen K-edge, where the core hole is in a Is (L 0) shell. 
In this type of RIXS experiment, as in optical Raman 
scattering, the microscopic coupling between the electro- 
magnetic radiation and the Sz ^ two-magnon excita- 
tion is due to a modulation of the strength of magnetic 
bonds. As we will review in the sections to come, the 
difference is that in RIXS the modulation is caused by 
the intermediate state core-hole modifying the magnetic 



exchange processes ( 


Van den Brink 


2007 Forte et al. 


2008a 


Hill et al. 


2008 


) . This gives rise to transition ma- 



trix elements that are very different from optical Raman 
scattering. Compared to optical probes the resolution 
of K-edge RIXS is far poorer, but it has the advantage 
of being able to probe the dispersion of the two-magnon 
excitations throughout the Brillouin zone. 

This section is organized as follows. In Sec. |V.ETT1 the 
various ways in which the RIXS process can excite mag- 
netic degrees of freedom are discussed. The recent theo- 
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retical and experimental progress is discussed in parallel 



in Sec. V.E.2 through V.E.4 where we focus on magnetic 
materials with interactions in 1, 2 and 3 dimensions. In 
this section we focus on purely magnetic excitations, i.e., 
magnetic excitations without, for example, accompany- 



ing orbital excitations as described in Sec. V.D 



1. Coupling to magnetic excitations witii RIXS 



tering, the RIXS cross section therefore consists of a local 
structure factor (depending on polarization, experimen- 
tal geometry and the excitation mechanism) multiplied 



by the appropriate spin susceptibility ( Ament et al. 2009 



Haverkort, 2009). For non-interacting spins, the suscep- 
tibility is uniform and featureless, but for systems with 
interatomic spin-spin interactions the susceptibility ac- 
quires a strong q-dependence, which may be measured 
in momentum-resolved RIXS. 



As indicated above, both direct and indirect RIXS can 
probe magnetic excitations, be it of different kind. In di- 
rect RIXS, single spin-flip excitations can be made at the 
2p, 3p — >■ 3d edges of Ni and Cu because of the large spin- 



orbit coupling of the 2p core-hole ( 


Ament et al. 


2009 


De Groot et al. 


1998 




Van Veenendaal| 2006|). For cop- 



per these edges split into intermediate states with total 
angular momentum j = 1/2 (L2,M2 edges) and j — 3/2 
(L3, M3 edges). The edges are well- resolved in XAS spec- 
tra (particularly at the L edge), and therefore in RIXS 
one selects either the intermediate states with j = 1/2 or 
j — 3/2. This implies that the spin and orbital angular 
momentum separately are no longer good quantum num- 
bers in the intermediate state. Therefore orbital and spin 
orbital angular momentum can be exchanged and direct 
spin- flip processes can in principle be allowed in RIXS, 
unlike in optical spectroscopy. 

For indirect RIXS at transition metal K-edges (Is — > 
4p), the core hole couples to the spin degree of free- 
dom by locally modifying the superexchange interac- 



tions (Van den Brink 2007 Forte et al. 2008a). In such 



a process, the total spin of the valence electrons is con- 
served, and only excitations where at least two spins are 
flipped (with total AS^ = 0) are allowed. A similar in- 
direct process can also occur at the oxygen K edge. At 
the copper and nickel L and M edges, both direct and 
indirect RIXS processes are allowed. 

It should be noted that the RIXS process can also lo- 
cally change the size of the magnetic moment. Transi- 
tions with AS* — 1 and 2 are allowed for certain mag- 
netic ions, for instance a transition of high spin Ni^''" 
into its low spin state. The energy scale of excitations 
with AS* 7^ is set by Hund's rule coupling Jh, which is 
typically of the order of 1 eV. These excitations therefore 
have a much higher energy than the A^* = ones, whose 
energy scale is usually set by the superexchange constant 
J - 100 meV. 

When single spin- flip scattering is allowed in RIXS, 
one can also measure the dispersion of magnons. This 
might not be obvious because magnons are after all not 
local spin-flips but highly delocalized magnetic excita- 
tions. But one should keep in mind that the incident 
photon can be scattered at any equivalent site, leading to 
a final state that is a superposition of spin-flips at equiva- 
lent sites. Such a final state carries a non-local magnetic 
excitation with momentum hq. Just as in neutron scat- 



a. Direct RIXS. Here we give an overview of the cases in 
which direct RIXS can cause spin-flip excitations without 
involving any additional change in the orbital component 
of the valence electrons. Combined excitations of that 
kind were discussed in Sec. IV. E. Pure spin-flip scattering 
creates a single magnon quasiparticle and causes a change 
of angular momentum of the spin system ASz = 1. Pro- 
cesses of this type were studied theoretically for Ni^"*" 



and Cu2+ by ( 


Ament et al. 


2009| |De Groot et al. 


1998 


Haverkort 


2009 


Van Veenendaal 2006 


)■ 


De Groot et al. 


(1998 


) showed numerically that there are no pure spin- 



flips present for a Cu^+ ion when the spin is aligned along 
the 2:-axis and the hole occupies the 3d^2 
can also be seen in Fig, 
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y2 orbital. This 
where the 3^2,2 _y 2 initial state 
has spin-up along the z-axis, and the ion is positioned in 
a crystal field. The prefactor for the effective op- 
erator that creates spin flips S± is — 4, which means 
that spin flips are not allowed for ^ = ±2 = — ,xy 



( Van Veenendaal 2006 ) , as is confirmed numerically. 

In contrast, for Ni^+ in an octahedral crystal field and 
spin parallel to the z-axis, pure spin-flip transitions with 
ASz = 1,2 are possible w ithout altering the orbit al part 
of the ground state (iDe Groot et al.\ |1998|). Note 



that the crystal field is not crucial in these observations: 
it only singles out the particular ground state used in the 
calculations. 



It was subsequently shown theoretically ( Ament et al. 



2009) and experimentally (Braicovich et al. 2010) that 



a pure spin-flip transition for a Cu^"*" ion with a hole in 
the 3d^2_y2 orbital does occur if the spin is not parallel 
to the z-axis. In this case, single spin-flip excitations are 
a result of the L^S^ operator in the core hole spin-orbit 
coupling. When looking at the effective 3c? transition 
operators, spin flips are a result of S^. This operator has 
a prefactor /x^ — 1 and spin flips are therefore allowed for 
— orbitals if the spin is in the xy plane (on the other 
hand, they are now forbidden for fi ~ ±1 zx, yz). For 
the cuprates these results put L edge RIXS as a technique 
on the same footing as neutron scattering. 



Using linear spin wave theory, Haverkort (2009) cal- 
culated the direct magnetic RIXS ^2,3 spectra for Cu^"*" 
{S = 1/2) and {S = 1) ions emlaedded in 1, 2 and 

3D systems. Although pure Heisenberg spin chains do 
not intrinsically order, Neel order could be stabilized by 
external forces like inter-chain coupling, and linear spin 
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wave theory would then be a vahd approach. The RIXS 
spectrum is factorized into a local transition operator and 



Initial 



Intermediate 



Final 



a spin susceptibility factor, following Anient et al. (20091. 



Haverkort found that the local operator can be described 
in terms of fundamental X-ray absorption spectra, like 
the magnetic linear and circular dichroic spectra. 

In principle the same holds at the Cu M edge (~75 
eV). However there the photon momentum is small, so 
that only magnons in a very small portion of the Brillouin 
zone can be probed. In general, for 3d transition metal 
ions, the 2p core-hole has both a spin-orbit coupling and a 
lifetime about an order of magnitude larger than 3p core- 
hole. Therefore, the spin-orbit split L2 and L3 edges are 
better resolved than the M2,3 peaks. Consequently the 
L edges offer a better opportunity to measure magnons 
with direct RIXS. 



b. Indirect RIXS. For indirect RIXS, double spin- flip 
excitations can be made because the core hole locally 



modifies the superexchange constant J (Van den Brink 



2007 Forte et al. 2008a I. This perturbation does not 



commute with the Hamiltonian of the magnetic system, 
and consequently magnetic excitations are created. 

This is best illustrated by considering in detail a su- 
perexchange process involving two ions that both have 
a singly occupied orbital. If an electron hops from one 
ion to the other a doubly occupied site is created, at the 
cost of a Coulomb energy U. In an exchange process the 
second electron subsequently hops in the opposite direc- 
tion. The two electrons are exchanged, and the effective 
amplitude for this process is, in second order perturba- 
tion theory, J oc /U , where t is the hopping amplitude 
of the participating electrons. The superexchange J cou- 



pling is modified by the presence of a core-hole ( Van den 



Brink 2007). 



At the transition metal K edge the photo-excited elec- 
tron ends up far (some eVs) above the Fermi level. The 
valence electrons tend to screen the core hole, see Fig. 35 
Screening affects the superexchange coupling: in the in- 
termediate state of the superexchange process, there is a 
doubly occupied site, and the electrons can be either both 
screening the core hole, or both occupying a neighboring 
ion. In the first case, the double-occupancy energy U 
is decreased by the core hole potential [4, while in the 
latter case it is increased by Uc- Because in this way U 
is effectively modified, J is also different at the core hole 
site. 

RIXS experiments at the Cu K edge of La2Cu04 and 
Nd2Cu04 showed a feature at 500 meV, which was at- 



tributed to two-magnon modes (Ellis et al. 2010 Hill 



et al. 2008). This discovery triggered theoretical work: 



detailed calculations of the magnetic response functions 
within the Ultrashort Core-hole Lifetime expansion (see 



Sec. IV. C. 3) were initiated in (Van den Brink 2007) and 



expanded upon in ( Forte et al. 2008a ) by including ef 
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FIG. 35 Indirect, magnetic RIXS mechanism for the Cu^^ K 
edge ( [Forte et aZ.[|2008a| . In the intermediate state, the 3d 
holes either screen or not screen the core hole. The Coulomb 
repulsion U is effectively replaced by U' , incorporating the 
effect of the core hole potential Uc. Note that the open circles 
represent holes, and the arrows indicate the spin direction. 



fects of longer range magnetic exchange interactions and 
finite temperature, which give rise to a generally weak 
singie-magnon contribution to the RIXS signal. Nagao] 



and Igarashi (20071 computed the magnetic RIXS spec- 
tra within this theoretical framework, also taking the 
magnon-magnon interactions into account. These inter- 
actions modify the detailed line-shapes of two-magnon 
RIXS spectra and cause a transfer of spectral weight to 
lower energies. 

At the transition metal L edges, a similar mechanism 
occurs (in addition to the direct single spin-flip scattering 
discussed above) because the photo-excited electron in 
the 3d subshell frustrates the local superexchange bonds. 
For Cu^"*" ions, the intermediate state configuration is 
3c?^° and superexchange is blocked completely: locally, 
J = (iBraicovich et 1 120091). 



Also at the oxygen K edge of transition metal oxides a 
similar modification of J is expected because the photo- 
excited electron is again disturbing the superexchange 
process. Transitions at the oxygen K edge are possible 
because of hybridization of the O 2p orbitals with the 
transition metal ions' unoccupied orbitals. 



2. ID systems - Chains and ladders 

When magnetic ions are Heisenberg-coupled together 
in a (ID) chain, collective magnetic excitations (spinons) 
emerge. Two parallel spin chains can be coupled to- 
gether to form ladders. If the couplings between the 
chains (i.e. the rungs of the ladder) are dominant, sin- 
glets will form on the rungs in the ground state. The sin- 
glet to triplet excitations on top of such a ground state 
are S = 1 triplons, which disperse because of coupling 
along the chains. The term triplon is used to distin- 
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FIG. 36 RIXS data after subtraction of the elastic signal 
( Schlappa et al. 2009 1 . The open black diamonds indicate 



the center of mass, obtained by fitting a Gaussian to the spec- 
tra. Theoretical one-triplon (red dashed line) and two-triplon 
(black full lines) dispersions are also shown. The right axis is 
scaled in units of J along the rungs, extracted from a single 
band Hubbard model. 



guish these three-fold degenerate elementary triplet ex- 
citations from magnons which are the elementary exci- 



tations of long-range ordered antiferromagnets ( Schmidt 



and Uhrig 2003 ) . There are many variants of these sys- 
tems, for instance, 'ladders' with three or more chains, 
but they are not yet investigated with RIXS. 

Sri4Cu2404i was the first one- dimensional system that 
has been studied with RIXS down to energies low enough 



to resolve the elementary magnetic excitations ( Schlappa 
et al. 2009). This compound has alternating layers of 



chains and ladders of spins. It superconducts under pres- 
sure when doped, and for that reason the system has at- 
tracted considerable interest, especially because it does 
not have the typical two-dimensional perovskite Cu02 
structure common to the high-T^ superconductors. The 
ladder has a singlet, spin-liquid ground state that can 
carry spin 1 triplon excitations. Superexchange interac- 
tions in the edge-sharing chains are weak because they 
are mediated by a 90° bond, and the resulting excitations 
cannot, at present, be resolved from the elastic peak in 
RIXS experiments. 



Schlappa et al. (2009) investigated the Cu L3 edge of 
Sri4Cu2404i. The authors see two-triplon excitations 
(on the spin-liquid ladder) which are gapped by 100 ± 30 
meV at the T-point, see Fig. [36] This region of the Bril- 
louin zone is inaccessible to inelastic neutron scattering 
due to lack of intensity, although it is possible to see the 
single triplon gap with neutrons, which is located away 
from the F-point. 

To investigate spin chains, Sr2Cu03 is an ideal com- 
pound. It is a corner-sharing spin chain system with 
short-range AFM order. The superexchange constant J 
along the chain is extremely large (~ 200 — 250 meV), 
so that excitations are easy to resolve in RIXS. Further, 
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FIG. 37 Single magnon RIXS data at the copper Z/3 edge of 
La2Cu04 with an energy resolution of 140 meV. The trans- 
ferred momentum (projected on the copper-oxide plane) is 
along the (vr, 0) direction. The decomposition of the low- 
energy spectrum is shown, and it is found that the maxi- 



mum follows the single magnon dispersion. ( Braicovich et al. 
20T0I. 



the inter-chain coupling is about an order of magnitude 
smaller, ensuring ID behavior. 

Seo et al. (2006) measured the RIXS spectrum of 
Sr2Cu03 at the Cu Is ^ id quadrupole transition, with 
an energy resolution of ^ 450 mcV. The authors observe 
a feature attributed to dd-cxcitations at an energy loss of 
2 eV. They claim that this feature disperses by ~ 0.2 eV 
with a periodicity of tt (half the lattice periodicity), and 
attribute this dispersion to two-spinon excitations on top 
of the dd-excitation. 



3. 2D systems - Magnons and Bimagnons 

Of the wide variety of (quasi) two-dimensional mag- 
netic materials, so far only the two-dimensional antifer- 
romagnetic cuprates have been studied experimentally 
with RIXS. 2D irdates have also been considered from a 
theoretical viewpoint (Ament et al. 2010). That the at- 



tention focusses on cuprates is not surprising, given their 
central role in high-Tc superconductivity. The undoped 
cuprates are described well in terms of 2D Heisenberg 
antiferromagnets. When doped with either holes or elec- 
trons, the magnetic order melts and superconductivity 
appears. In undoped cuprates the elementary low-energy 
magnetic excitations are well separated from the gapped 
charge excitations, but upon doping these excitations mix 
and thus so will their RIXS response. 

RIXS can create single and double spin-flip excitations, 
which are equivalent to single and bimagnon excitations 
in Heisenberg antiferromagnets. Both have been ob- 



served (Braicovich et al. 2010 Hill et al. 2008), and are 



reviewed below. In addition also three-magnon scattering 
is predicted to contribute substantially to the magnetic 



spectral weight (Ament and van den Brink 20101 
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FIG. 38 Comparing single magnon neutron data and RIXS 
measurements at the copper L3 edge of La2Cu04. Details as 
in Fig. [37| ([Braicovidi et ar||2010[). 



There are many more 2D materials of interest that so 
far have not been studied with RIXS, like frustrated mag- 
nets, doped cuprates, spin-orbit coupling dominated ma- 
terials, etc. This is not only due to scarcity of beam-time, 
but also because for RIXS materials with large magnetic 
exchange integrals causing relatively high-energy mag- 
netic excitations are preferable in view of the restricted 
energy resolution. 



a. Single magnon. In many two-dimensional cuprates, 
single spin-flip excitations are allowed in direct RIXS, i.e., 
at the Cu L and M edges (Ament et al. 2009). For in- 
stance, La2Cu04, CaCu02, Sr2Cu02Cl2, and Nd2Cu04 
all have the copper spins (residing on the 'idr^-2_y2 orbital) 



oriented in the xy-plane. The above-mentioned materials 
are two-dimensional Heisenberg antiferromagnets, and in 
such systems a spin-flip is described in terms of the single 
magnon final states. Also in strong spin-orbit coupled 
iridates, e.g., Sr2lr04, RIXS can probe the elementary 
magnetic excitations ( Ament et al. 20To[ ) . 



The direct RIXS cross section consists of the atomic 
spin-flip factor multiplied by the spin-flip susceptibility, 
which accounts for the collective response. The atomic 
factor depends on the polarization vectors and the spin 
orientation. Haverkort 2009 relates these atomic factors 
to fundamental x-ray absorption spectra. The suscepti- 
bility factor depends only on momentum transfer. In the 
antiferromagnetic cuprates, the latter factor peaks at the 
magnetic ordering vector q — (tt, tt), while going linearly 
to zero around the F-point. 



It is remarkable that RIXS allows the measurement of 
a magnon dispersion such a thin film sample with high 
accuracy. 



Guarise et al. (2010) measured the magnon disper- 
sion in 2D antiferromagnetic insulator Sr2CuCl202. The 
spin-waves display a large (70 meV) dispersion between 
the zone-boundary points (7r,0) and (7r/2,7r/2). From an 
analysis of the dispersion in terms of t-t'-t"-U Hubbard 



model, Guarise et al. ( 2010 ) conclude there is a significant 



electronic hopping beyond nearest-neighbor Cu ions. 



Braicovich et al. (2010) also studied the RIXS response 
of underdoped La2-2:Sra;Cu04 with a hole concentration 
of a: = 0.08. Due to the observed asymmetry in the cross 
section when changing q to — q, the RIXS spectrum could 
be assigned predominantly to magnetic excitations, fol- 



lowing theory (Ament et al. 2009 ). At low temperatures. 



two dispersing branches are observed: the high-energy 
one (not previously observed) follows the magnon dis- 
persion of the undoped compound within 10%, while the 
low-energy one coincides with neutron measurements of 
doped systems around (7r,7r). The authors caution that 
q = (7r,7r) no longer has to be equivalent to (0,0) since 
there is no long range magnetic order at this doping. 
The two branches merge at high temperature. |BraicovicE| 
et al. (2010) suggest that the two branches signal phase 
separation on a microscopic scale. These branches are 
visible down to a transfered momentum of 0.44 A~^, and 
the authors argue that sustaining coherent excitations of 
such momentum, requires magnetic patches much larger 
than the corresponding length scale of 15 A. The patch 
size has also been estimated with other techniques that 
have a much longer measuring time scale. These mea- 
surements find a smaller patch size, leading the authors 
to conclude that the underdoped compound is dynami- 
cally phase separated. 

2. Iridates The relativistic spin-orbit coupling in 
Ir^+ dominates the spin-orbital dynamics, and the ionic 
ground state is a Kramers doublet which can be described 
by a pseudo-spin S = 1/2. In Sr2lr04, the t2g analog 
of La2Cu04, these spin-1/2 moments couple by superex- 
change and yield an antiferromagnetic Heisenberg model. 



Ament et al. (2010) have calculated the Ir L-edge RIXS 
cross section for Ir'^"'" ions in an octahedral field. They 
have theoretically shown for Sr2lr04 that, next to transi- 
tions between spin-orbit multiplets, RIXS can also probe 
the magnon spectrum, as shown in Fig. |40) At the L2- 
edge, no intensity is expected, in contrast to the L^-edge. 



1. Cuprates Single-magnon excitations were ob- 
served in RIXS experiments at the Cu L3 edge on a thin 
film of La2Cu04 by Braicovich et al. (2010). The au- 



thors found that the maximum of the magnetic RIXS 



spectrum (Fig. 37 1 follows the magnon dispersion mea- 



sured with neutron scattering very closely, see Fig. 38 



b. Bi-magnon. For direct RIXS the projected spin an- 
gular momentum of the material may be changed by 
1 or 2 units by a transfer of angular momentum of the 
photon to the electron spins in the system. For indirect 
RIXS such a direct transfer of angular momentum is not 
possible and the spin angular momentum of the mate- 
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FIG. 39 The RIXS spectra of the underdoped cuprate 
Lai,92Sro.o8Cu04 shows two dispersing branches at low tem- 
perature (a), which merge at room temperature (b). The 
green line in Fig. (c) displays the corresponding neutron data. 
The data were taken at the copper L3 edge with an en- 
ergy resolution of 140 meV. The transfered momentum (pro- 
jected on the copper-oxide plane) is along the (tt, 0) direction 
|Braicovich et a/.,|2010( ). 
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FIG. 40 Theoretical RIXS spectrum of Sr2lr04 at the Ir L- 
edge for tt and a polarization in 90° scattering geometry, and 
where q is in the first Brillouin zone. The broad peak at 
lOJ is due to transitions between the spin-orbital multiplets. 
From this peak one can determine the spin-orbit coupling A. 
At low energy, strongly dispersive single- and double-magnon 
features are seen (Ament et al. 20101. 



rial cannot be changed in the scattering process. There- 
fore, the leading order allowed magnetic excitations that 
result from indirect processes are bimagnon excitations. 
These have been observed with RIXS in numerous two- 
dimensional cuprates and are theoretically well under- 
stood. 

The first experimental observation of the momentum 
dependence of bimagnon excitations in the cuprates, at 



the copper K edge, was performed by Hill et al. (20081. 
The measurements were done at T = 20 K, and with 
an energy resolution of 120 meV. The authors measured 
indirect RIXS spectra at the Cu Is Ap^c edge on 



La2Cu04 and Nd2Cu04, and both materials showed a 
peak at momentum transfer q = (tt, 0) around an en- 
ergy loss of 500 meV, which they interpreted as a bi- 
magnon excitation. This resonant feature disappeared 
quickly with hole doping, similar to bimagnon Raman 
scattering. Further, when the polarization was changed 
from parallel to the z-axis to perpendicular to it, the fea- 
ture vanished. The authors measured at five points in 
the Brillouin zone, and found zero intensity at q = (0, 0) 
and (7r,7r), while the bimagnon DOS is highest at these 
points. The RIXS intensity peaked at q = (tt, 0), see 
Fig.|4l] 

In a follow-up study, the energy and intensity of the 
q = (tt, 0) peak in La2Cu04 at the Cu if -edge are 
tracked as a function of temperature and doping, pro- 
viding evidence that the peak is a two-magnon excita- 
Data was taken at three beam- 



2010) 



tion (Ellis et al 
lines. At the APS 30ID beamline, the incident photons 
were 7r-polarized and the scattering geometry ^ 90°. A 
combined energy resolution of 120 meV (FWHM) was 
achieved. At the APS 9ID beamline, the incident photons 
were u-polarized and the combined energy resolution 130 
meV (FWHM). Measurements at SPring-8's BLllXU 
beamline were done in the same geometry as at the 30ID 
beamline; the resolution was 440 meV. The q-dependence 



of the peak as found by |Hill et al. ( 2008 1 is confirmed, see 
Fig. [41] The authors find that the intensity of the peak 
decreases with temperature, vanishing only at T = 500 
K, far above the Neel temperature = 320 K. This be- 
havior is similar to Raman data on various cuprates, and 
also consistent with neutron measurements. The width 
of the peak was estimated to be ~ 160 meV, although 
it could not be established firmly because of the energy 
resolution of 120 — 130 meV. The peak shift as a func- 
tion of temperature is identical to Raman data on other 
cuprates when the temperature is renormalized to the 
Neel temperature. The authors argue that this is evi- 
dence that magnon-magnon interactions are the source 
of this shift. The decreasing intensity on the other hand 
is attributed to phonons. Doping with Sr also causes the 
500 meV peak to shift to lower energy and lose inten- 
sity rapidly following a phenomenological formula. One 
problem in attributing the peak to two-magnon excita- 
tions is its energy: 500 meV where interacting magnon 



theory predicts 400 meV. Ellis et al. (2010) suggest a 



strongly momentum-dependent magnon-magnon interac- 
tion to explain this fact: it should be strong in the Raman 
case of q = (0, 0) and weaker at (tt, 0). 

From the theoretical side, it was noticed already very 
early on that excitations related to spin degrees of free- 
dom, such as two-magnon Raman scattering, are present 
in RIXS spectra computed for small Hubbard Hamilto- 



nian clusters for RIXS on the K-edge of copper ( Tsutsui 



et al. 19991 and oxygen (Harada et al. 2002). Van den 



Brink ( 2005 2007 ) proposed the microscopic mechanism 



by which the bimagnons couple to the intermediate state 
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FIG. 41 RIXS spectra at the copper Is — > 4p edge for a 
few points in the two-dimensional Brillouin zone of La2Cu04 



(Ellis et al. 2010 Hill et al. 20081. Fig. (b) is taken with 



TT-polarization at room temperature, (c) with cr-polarization 
at 45 K. The excitation energy is 8994 eV. The solid lines 
are Lorentzians fitted to the data. The elastic peak has been 
subtracted from these data. 



core-hole in RIXS 



by locally modifying the superex- 
This approach 



change constant J (see Sec. V.E.I I. 
yielded a momentum dependence of the cross section 
in agreement with experiment and in particular repro- 
duced the lack of intensity at q = (0,0) and (tTjTt). The 
magnetic RIXS response within this framework includ- 
ing magnon-magnon interactions was computed by |Na-| 
gao and Igarashi ( 2007 1 and is shown in Fig. 42 Magnon- 



magnon interactions are observed to broaden the loss fea- 
tures and and at the same time enhance lower-energy 
spectral weight with respect to the non-interacting the- 



ory of Van den Brink (2005 20071 



Later theoretical work took the approach of extending 



1963 Fleury and Loudon 1968) to finite momentum 



transfer (Donkov and Chubukov 2007 Vernay et al. 



2007 ) . This approach is in essence non-resonant because 



it omits the core-hole (and its effect on the valence elec- 
trons) from the analysis. 



Vernay et al. (20071 consider a single-band, half-filled 



Hubbard model, and assume that the incoming x-ray 
photon excites a valence electron to a neighboring site, 
making it a doubly occupied one in the intermediate 
state. In the final state, one of the two electrons on 
the doubly occupied site has hopped back, emitting the 
outgoing photon. They classify spectra according to po- 



bimagnon Raman scattering theory (Elliott and Loudon [eFoI] 120091) 



larization geometries and find that the fully symmetric 
component transforming as Aig vanishes at q = (0, 0) 
and both Aig and Big geometries vanish at (7r,7r). For 
Big geometries they identified how the bi-magnon inten- 
sity disperses from the zone center (Raman) at ^ 2.8J 
to 4 J at (7r,0) due to the reduction of magnon-magnon 
interactions with increasing q, while the Aig bi-magnon 
intensity is strongest at (tt, 0) appearing at 4 J. 



Donkov and Chubukov ( 2007 ) take a similar approach 



Starting from a half-filled single-band Hubbard model, 
they make a mean-field approximation and introduce an 
antiferromagnetic order parameter, obtaining an upper 
and lower Hubbard band. They also assume that, instead 
of making core transitions, the X-rays excite electrons 
from the lower Hubbard band to the upper one. These 
electron-hole excitations couple to magnons via the full 
Hubbard U interaction, which was initially treated on a 
mean-field level. The authors find that in the Big exper- 
imental geometry the intensity peaks either at q = (0, 0) 
for small magnon-magnon interactions or at (7r,0) for 
large interactions. The bimagnon peak disperses to lower 
energies away from q = (0,0). For the Aig geome- 
try, both the intensity and dispersion go up away from 
q = (0,0). The intensity peaks at (tt, 0) and is zero at 
(0,0). For both geometries, the authors predict a finite 
intensity at q = (7r,7r). 

calculated the RIXS bimagnon 



Forte et al. (2008a 



scattering via the core-hole at finite temperatures, con- 
sidering in addition longer-range superexchange interac- 
tions, comparing the theory in detail with the K-edge 



experiments of Hill et al. (2008). The authors employed 



the Ultra-short Core-hole Life-time (UCL) expansion in 
the framework in which the bimagnons are created by 



the core-hole (Van den Brink 2007) and also calculated 



the sub-leading order of the UCL expansion, which leads 
to finite bimagnon spectral weight at q = (0, 0). 

At the Cu L and M edges, the coupling mecha- 
nism to bimagnon excitations is different, as clarified in 
Sec. |V.E.1[ but the resulting cross section is the same to 
first order in the UCL expansion, see Fig. 42 This comes 



about because, instead of modifying the superexchange 
constant to some new value J', superexchange is locally 
blocked at the Cu L and M edges: J' = (Braicovich 



Measuring 
along the q = 



at many different transferred momenta 
(0,0) — (tt, 0) direction, (Braicovich et al. 



2009 ) report bimagnon excitations at the Cu L3 edge of 
La2Cu04 and CaCu02. At T = 30 K and at room tem- 
perature (which are both below the Neel temperature), 
they see a feature dispersing from meV at q = (0, 0) to 
approximately 400 meV at q = (27r/3,0) in both a and 
TT polarization. The authors assign it to a bimagnon ex- 
citation, which disagrees with theoretical work that pre- 
dicts that a strong single magnon peak should be present 
as well (Ament et al. 2009). In following experimental 



work, where the energy resolution was improved from 400 
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FIG. 42 Bimagnon RIXS spectra obtained with interacting 
spin wave theory for a nearest-neighbor S = 1/2 Heisenberg 
model ( Nagao and Igarashi 2007[ ) . The vertical axis indicates 
transferred momentum projected to the Cu02 plane and uj is 
the energy loss. For the 2D square lattice, z = 4. 



meV to 140 meV, the authors reinterpret the dispersing 



feature as a single magnon (Braicovich et al. 2010). The 



remaining spectral weight at slightly higher energies is 
now interpreted as bimagnon excitations. 

The effect of chemical pressure on 2D antiferromag- 
netic cuprates is studied by comparing the Cu L3 RIXS 
spectra of CaCu02 and BaCu02 (Bisogni et al. 2009). 



By substituting Ba for Ca, the in-plane lattice parameter 
increases from 3.84±0.01 A to 3.96 ±0.01 A. This causes 
the hopping integrals to decrease, and the authors ob- 
serve that, correspondingly, the superexchange constant 
in BaCu02 is 0.67 ±0.08 times the CaCu02 value. They 
extract this ratio by comparing the position of the mag- 
netic peak in the two compounds, which scales linearly 
with J . Putting J at 120 meV in CaCu02, they obtain a 
superexchange constant of 85±10 meV for BaCu02. The 
spectra are normalized so as to keep the integrated weight 
of the dd-cxcitations the same in both compounds. With 
such a normalization, the spectral weight of the magnetic 
peak in the Ba compound is found to be to less than 
half the Ca value. This is in agreement with what one 
would expect for bimagnons, since the bimagnon spectral 
weight scales as . The magnetic peak was therefore in- 
terpreted by the authors as a bimagnon peak, although 
they later changed this interpretation to single magnon 



(Braicovich et al. 2010). 



At the Cu M edge, experiments are more difficult be- 
cause of the large elastic peak, even though the experi- 
mental resolution is usually higher than at other edges. 
RIXS spectra taken at the Cu M2 edge of CaCu02 with 
a resolution of 140 meV show two low-energy features 
(Freelon et al. 2008). In the data dd-excitations are sup- 



pressed and it is uncertain what causes this effect. The 
authors observe a loss feature at 390 meV when the in- 
coming polarization is parallel to the c-axis. This being 
a single magnon excitation can be excluded as at the M 
edge q « and for these momenta the single magnon 
energy vanishes. The authors interpret it as a bimagnon 
excitation. 

The bimagnon excitation can in principle also be ob- 
served at the oxygen K edge (Is — ?> 2p transition). Al- 
though the cross section for the oxygen K edge is rel- 
atively small, this case is interesting because the single 
magnon excitation is forbidden here: there is no spin- 
orbit coupling for Is core orbitals. In this sense the situa- 



tion is similar to the copper K edge. Harada et al. ( 2002 ) 



investigated Sr2Cu02Cl2 at the oxygen K edge with an 
energy resolution of 0.5 eV. The transferred momentum 
is fixed by the 90° scattering geometry and near-normal 
incidence of the x-rays. For scattering with a polariza- 
tion a feature is visible around 0.5 eV which the authors 
interpret as a as a two-magnon excitation. A calculation 
for a CU4O13 cluster reproduces the feature at 0.5 eV 
and its polarization dependence, but underestimates the 
intensity, possibly due to finite size effects. 



4. 3D systems - Magnons and Bimagnons 

A number of 3D magnetic Mott insulators have been 
studied with RIXS: in particular NiO, NiCl2 and CoO. 
The prototypical strongly correlated antiferromagnet 
NiO has been investigated most extensively, and we 
therefore start by briefly summarizing its magnetic prop- 
erties. 

NiO has the same crystal structure as NaCl, with Ni^+ 
ions (3^^^ configuration) centered in an octahedron of oxy- 
gen ions, generating a crystal field of that same symme- 
try. In the ground state the t2g levels of Ni^"*" are fully 
occupied, with the remaining two electrons occupy the 
eg orbitals, with parallel spins to that the Ni^+ ion has 
an S* = 1 ground state. The spin points along the (11-2) 
direction. Superexchange interactions mediated by the 
oxygen ions give rise to four (weakly coupled) AFM sub- 
systems, ordering at Tm — 523 K and described at low 
temperature by an S' = 1 Heisenberg model on a cu- 
bic lattice. Magnons are the elementary excitations of 
this Heisenberg antiferromagnet. The magnon density 
of states is strongly peaked around 6 J, where 6J = 114 



meV in NiO (Hutchings and Samuelsen 1972) 



Atomic calculations for a Ni^"'" ion in an octahedral 



crystal field by De Groot et al. ( 1998 ) show that at the 
Ni L-edge AS'z = 1,2 transitions are possible. These 
RIXS final states correspond to single (AS'z = 1) and 
bimagnon (AS'^ = 2) excitations. 

Final states with AS'^ = 2 are interesting because 
they cannot be reached with neutron scattering, optical 
Raman scattering or indirect RIXS. In addition to the 
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ASz = 2 bimagnons, also ASz = bimagnons can be ex- 
cited in NiO via the indirect RIXS process. In this case 
the bimagnon scattering comes about by the core-hole 
in the intermediate state locally disrupting the magnetic 
bonds, as discussed in Section [V.E.l.b| for the cuprates. 

Besides the mentioned low-energy magnon excitations 
- in which only the z-projections of the local spin mo- 
ments changes on the Ni sites, RIXS can also cause local 
spin transitions, i.e. a transition from a high-spin triplet 
{S = 1) state of Ni^+ to a low spin singlet S* = state. 
Such spin state transitions are at much higher energy 
than the magnons as they are at an energy scale set by 
the atomic Hund's rule coupling Jh ~ 1.3 eV. 



a. (Bi)-Magnons. The first measurement of magnetic ex- 
citations in NiO with RIXS was performed on the nickel 



M edge by Chiuzbaian et al. ( 2005 ) . Because of the over- 



whelming elastic peak, the observation of pure spin-flip 
transitions was impossible with their energy resolution of 
130 meV. 

In follow-up work on NiO, it proved to be possible to 
observe two purely magnetic excitations at ~ 95 and 190 



meV ( Ghiringhelli et al. 2009 ) . These experiments had 



nearly the same energy resolution (120 meV), but prof- 
ited from the fact that at the Ni L3 edge the elastic peak 
is much smaller than at the M edge. These excitations 
were not observed to disperse. 

To establish that the peaks are of magnetic nature, the 
authors also measured NiCl2. Like NiO, NiCl2 also has 
Ni^"*" ions in a (slightly weaker) octahedral crystal field 
of Cl~ ions. The superexchange constant is an order of 
magnitude smaller than in NiO, and therefore any mag- 
netic peak will be buried in the elastic line, while crystal 
field excitations should remain almost unchanged. The 
two low-energy peaks seen in NiO are absent in the NiCl2 
spectra, indicating that they are indeed magnetic excita- 
tions. 

The authors model the RIXS spectrum using a single 
Ni^"*" ion in an effective exchange field. The resulting 
spin-fiip transitions take the 5 = 1 ground state into an 
S — 1 final state, where the projection on the magneti- 
zation axis a can change with AS'q = 1,2. This local 
model describes the data reasonably well, although the 
single spin-flip peak deviates from the calculations. 

Measuring the dispersion of magnon and bimagnon ex- 
citations in a 3D system is challenging, but a dispersion 
should be present just as it is in the magnetic 1 and 
2D cuprates. The single magnon dispersion is expected 
to be similar to the one observed with neutron scatter- 
ing. The bimagnon signal (with AS'z = 2) is found to 
hardly disperse in a model calculation on a 3D spin one 
nearest-neighbor Heisenberg model on a cubic lattice. 
Rather it peaks around 12 J, following the bimagnon DOS 



(Haverkort, 2009) 



been published of the 3D antiferromagnet LaTiOa and 
ferromagne t YTiOg with Tjy = 146 K and Tc = 27 K re- 
spectively (Ulrich et al. 2009). In both compounds, the 



Ti"*" ions are in the 3(i^ configuration: the spin S = 1/2. 
Spectra with an energy resolution of 55 meV were taken 
at the Ti L3 edge at room temperature and at 13 K, 
i.e., above and below the magnetic ordering tempera- 
tures. The authors point out that in the magnetically 
ordered state, the low-energy region is dominated by a 
single magnon peak at approximately 40 meV for LaTiOa 
and 16 meV for YTiOa. The energies agree with neutron 
scattering measurements. This magnetic peak could only 
be distinguished from the elastic line in 90° scattering ge- 
ometry with TT-polarization. At room temperature, where 
both compounds are paramagnetic, the peaks shift to 22 
meV for LaTiOa and zero for YTiOa. 



b. Local Spin Transitions. NiO has also been studied at 
the O K edge, where direct RIXS spin-flips are forbidden 
because of the absence of spin-orbit coupling in the 1s 
core orbital. 



Duda et al. (2006) measure an excitation 



around 1.9 eV, in cr-polarization only. The authors inter- 
pret this feature as a 'double singlet' excitation on two 
neighboring Ni ions: Ittiii) ^ Itiiit), leaving both Ni 
ions in an 5 = state. The energy of this excitation 
- twice Hund's rule coupling Jh ~ 1-3 eV - should be 
offset by 0.7 eV to account for the hybridization energy. 
Note that this energy scale is much higher than the ~ 190 
meV for the bimagnon excitation seen at the Ni L3 edge 



(Ghiringhelli et al. 2009). The polarization dependence 



of the double singlet excitation seen in experiment agrees 
with calculations on a NieOig cluster in which the cen- 
tral O ion is excited by an X-ray photon. This model 
yields a double singlet only in u-polarization, but not in 
TT-polarization. 



Going one step back in the periodic table, Magnu- 



son et al. ( 2002b ) and Ghiuzbaian et al.\ ( 2008 ) studied 
CoO. The Co^"*" ions are in a 3d' configuration and in 
the ground state assume spin S = 3/2. A cubic crys- 
tal field {lODq « 0.5 eV) combined with Hund's rule 
selects the configuration t^g e^, and the orbital ground 
state has Tig symmetry with substantial charge-transfer 
character. The spectrum consist of several spin-quartets 
and -doublets within ~ 3 eV ( Magnuson et al. 2002b ) . 



Chiuzbaian et al. (2008) measured Co M2,3 and Co L3 
RIXS, while Magnuson et al. (2002b) measured the Co L3 



edge only, see Fig. |43) The lowest energy excited states 
(below ^ 3 eV) are dominated by crystal field states de- 
rived from atomic spin-quartet states, while above ~ 3 
eV, states derived from spin-doublets dominate the spec- 



trum (Chiuzbaian et al. 2008). On the theory side, both 



authors simulate the RIXS and XAS spectra within a sin- 
gle impurity Anderson model including multiplet effects. 



Chiuzbaian et al. (2008) performed a calculation on a 



Of the early transition metal oxides, RIXS spectra have CoOg cluster including 3d spin-orbit coupling but with- 
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FIG. 43 Experimental (dots) and theoretical (solid line) RIXS 



spectra of CoO at the Co L edges ( Magnuson et al. 2002b I 



The theoretical curve is a broadening of the vertical bars. The 
excitation energy runs from 778.1 eV (A) to 785.2 eV (E) at 
th e Z/3 edge, and 795 . eV (F) is at the L2 edge. The vsrork 
of Chiuzbaian et al. (20081 includes similar spectra with a 



slightly better combined energy resolution. 



out exchange field. Including exchange (9J — 0.3 eV) 
improves agreement between theoretical and experimen- 
tal RIXS spectra, as shown by iMagnuson et al.\ ( 2002b I , 



but reduces the agreement for the Co L3 XAS spectra. 
All experiments were performed at room temperature, 
which is far above the Neel temperature Tjv ~ 288 K. 



VI. SUMMARY AND OUTLOOK 

In this review, we have surveyed Resonant Inelastic 
X-ray Scattering with a focus on the creation of ele- 
mentary excitations. Recent years have seen a tremen- 
dous improvement in instrumentation for the sources, for 
the beamlines and for the spectrometers. As a result, 
there has been a tremendous advance in the technique. 
Whereas early work primarily dealt with excitations on 
the scale of electron volts, such as dd and charge-transfer 
excitations, more recent experiments have resolved low- 
energy features such as magnons. 

It would be interesting to see if this can be pushed 
down to electronic pairing gaps, such as the supercon- 
ducting and charge-density wave gap, or more exotic 
ones, such as pseudogaps in cuprates. RIXS provides 
an alternative and complementary look at the charge dy- 
namics of a material. In the coming years, the experimen- 
tal RIXS work in the L and M edge regions likely will 



be focused on establishing itself as a viable technique to 
measure elementary excitations at low energies, and in- 
vestigate in detail the momentum dependent elementary 
charge response, viz magnons, orbitons and phonons, of 
correlated and uncorrelated materials. At the K edge, 
the majority of the work has been performed on cuprates 
and this research is likely to continue. However, more- 
detailed extensions to other transition-metal compounds, 
such as manganites, pnictides, and ferroelectric materi- 
als, can undoubtedly be expected. Furthermore, the work 
on other strongly correlated systems, such as the edges 
of the rare-earths, is currently extremely limited and is 
expected to grown. In addition, with the advent of new 
synchrotron sources, there will be a push to improving 
the experimental resolution towards the 10 meV level. 
As has been demonstrated many times to date, each im- 
provement in resolution brings with it the discovery of 
new details in the excitation spectrum of these materi- 
als. In addition to searching for electronic gaps, men- 
tioned above, it would also be of particular interest to 
challenge theory to reproduce the finest of details seen in 
these response functions with this improved resolution. 

Apart from bulk samples, it would also be interesting 
to make a comparison of lower-dimensional systems, such 
as thin films, interfaces, and nanostructured materials, 
which have not been covered in this review. Each of 
these are exciting areas in which new excitations and new 
physics is to be anticipated. RIXS is likely the only tool 
that can probe these systems over the relevant energy and 
momentum scales. Another fruitful area of growth for the 
technique will be in the study of materials where other 
spectroscopies may not be able to provide information on 
elementary valence excitations. Examples here include 
complex materials where the chemical selectivity aids in 
the analysis, or materials in extreme environments, such 
as high pressure or strong electric/magnetic fields, that 
are difficult to access with techniques such as neutron 
scattering and photoemission. 

Concurrent with the experimental advances, theory 
has made great progress in analyzing the RIXS cross 
section. Determining the charge dynamics in strongly 
correlated systems is a deep and complex problem. In 
elucidating this problem, it is very helpful to have effec- 
tive methods that establish the correlation function that 
RIXS measures under certain conditions, and its relation 
to other spectroscopic techniques. Additional theoreti- 
cal methods can aid in providing a deeper understand- 
ing of the charge dynamics which is relevant in under- 
standing the optical and conductive properties of mate- 
rials. Here advancements in the theory would include 
techniques which do not suffer from finite size effects, are 
able to address many intertwined degrees of freedom, and 
have as few approximations as possible. 

The core-hole present in RIXS represents in many cases 
a strong local perturbation to the electronic system. As 
a result, dynamical mean-field theory (DMFT) may be 
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shorter than what can be produced by any other x-ray 



source available now or in the near future, see Fig. 44 



FIG. 44 Left: Peak brightness and Right: Average bright- 
ness of present and future x-ray sources. Present third gen- 
eration synchrotron light sources average brightnesses are in 
the I0^° phot/ s/Q.l7oBW/mm^ /mradP range. A brilliance 
enhancement by two orders of magnitude is envisaged for ul- 
timate x-ray storage ring-based light sources and energy re- 
covery linacs (ERLs). Free electron lasers (FELs) based on 
self-amplification of spontaneous emission (SASE) are pulsed 
and can enhance the peak brightness by ~ 8 orders of magni- 
tude (see|Saldin et al. (20001 and references therein). Feasi- 



bility studies show that X-Ray FEL Oscillator can greatly en- 
hance the average x-ray brightness (Kim et al. 20081. Cour- 
tesy of K.-J. Kim and R. Lindberg 



an important technique to apply to this scattering pro- 
cess as it can treat local dynamical many-body effects ex- 
actly. DMFT (and its cluster extensions) is well suited to 
include both local and non-local electronic correlations, 
and for certain model Hamiltonians this can be done ex- 
actly. As a result the multi-particle correlation functions, 
including multi-particle vertex renormalizations, can be 
formulated to calculate various spectroscopies, such as 
RIXS, as well as x-ray absorption and emission, and core 
level photoelectron spectroscopies. Progress has been 
made for non-resonant inelastic x-ray scattering across 

2003a|b 



a metal-insulator transition (Devereaux et 



Freericks et al. 2003) and in a CDW phase (Matveev 



et al. 2009). Moreover recent extensions to treat reso- 



nant Raman scattering with the Hubbard bands ( Shvaika 



et al. 2004 2005 1 point towards a future applicability for 



RIXS. One of the key ingredients will be the development 
of a formalism within DMFT to calculate the core-hole 
propagator. It is to be hoped that such development lies 
in the near future. 

One of the newest developments is the extension of 
RIXS into the time-domain. The very recent availabil- 
ity of ultra-short and ultra-bright x-ray pulses generated 
by free electron lasers such as FLASH at DESY and the 
Linac Coherent Light Source (LCLS) at SLAG National 
Accelerator Laboratory offers unique opportunities for 
investigating the temporal dynamics of electronic states 
using time-resolved x-ray scattering and spectroscopies. 
Pulses of x-ray laser light from LCLS are many orders of 
magnitude brighter and/or several orders of magnitude 



These characteristics will enable frontier new science in 
areas that include discovering and probing new states of 
matter, understanding and following chemical reactions 
and energy conversion processes in real time, imaging 
non-crystalline biological materials at atomic resolution, 
and imaging chemical and structural properties of mate- 
rials on the nanoscale. In addition, the strong coherence 
of these sources open the door to performing coherent 



stimulated RIXS ( jHarbola and Mukamel[ |2009 1 . 

From a theoretical point, concepts developed for RIXS 
are close to pump-probe experiments. In fact, a num- 
ber of groups have advocated the concept of viewing 
RIXS as an effective pump-probe experiment, where the 
absorption is the pump and the emission is the probe. 
The core-hole lifetime is the effective clock. These time- 
dependent measurements require the inclusion of non- 
equilbrium dynamics in spectroscopy. Although current 
theoretical frameworks describe charge and spin dynam- 
ics in the discussion of RIXS, after excitation, the total 
energy of the system is always conserved. When study- 
ing the time domain, the system will dissipate energy by 
coupling to the surroundings. The inclusion of nonequil- 
brium effects in x-ray spectroscopy is still largely an un- 
explored field. 

Other source developments may provide additional op- 
portunities for inelastic x-ray scattering techniques. Re- 



cent feasibility studies ( Kim et al. 2008 ) for X-ray free- 



electron oscillators, first proposed by Colella and Luccio] 
(1984), show the possibility of producing x-ray beams 
with a bandwidth of 2 meV in the 10 keV region. Re- 
lated possiblities include so-called seeded-FELs, which 
also produce transform limited photon pulses at high rep- 
etition rates. The presence of such sources would come a 
long way towards solving the major experimental RIXS 
issue of a limited photon flux, holding out the promise 
of many orders of magnitude (up to 10^) gain in incident 
flux. This is a tantalizing prospect allowing one to imag- 
ine studying single magnons in small exchange systems, 
full dispersion curves, superconducting gaps, etc. 

Summarizing, with the current steady progress in in- 
strumentation, combined with the increasing theoretical 
expertise, the use of resonant inelastic x-ray scattering 
techniques is becoming increasingly prevalent. The fu- 
ture for RIXS as a characterization tool able to probe 
elementary excitations in real materials under real con- 
ditions looks bright. 
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Appendix A: Spherical Tensors 

The use of spherical tensors is prevalent in x-ray 
spectroscopy resulting from expansions around the sites 
where the absorption takes place. The angular depen- 
dencies are typically expressed in terms of spherical har- 
monics Yimi^t, 'fit)- Since the angles 0} and (pi- deter- 
mine the direction of the unit vector f = r/r, we can 
also use the shorthand Y/m (f ) for the spherical harmonic 
in the direction of f. The spherical harmonics for a 
particular I form a tensor C'(f) with normalized com- 
ponents C^(f) = y^4:TT /[l]Yi„i{r) using the shorthand 
[l^.. ■ ;„] = (2/i -t- 1) • • • {2ln + 1). The rank of the tensor 
is I and its components are m = 1,1 — 1, ■ ■ ■ , —I. Note 
that the components of a vector form the tensor rC^(f). 

It is often useful to couple tensors together to from a 
new tensor. This can be done via the tensor coupling 

[T'',v']- = ^<'-,,(r^,)^KL, (Ai) 

mm' 

where the tensors V' and T' of rank I and V, respectively, 
are coupled to a new tensor of rank x. The coefficients 
^mm'q defined as 

-L% = (-l)''-™'(j;, ^^)n,TL (A2) 



the term between parentheses is a 3j symbol, which is di- 



rectly related to the Glebsch-Gordan coefficients (Brink 



and Satchlerl [19621 [Valshalovich et al| [19881 [Yltsis et al. 



1988). We use normalized tensor products that have the 



advantage that more nicely defined products are obtained 



when dealing with physical properties ( Thole and van der 



Laan 1994 Thole et al. 1994 1. This is achieved by defin- 



ing the normalization constants 
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{l + iy. l{L-2l)\{L-2l')\{L-2x)\ 

with L = l+V +x. Important examples of tensor products 
are the inner product [T', V']q = T' • V' and the outer 
product [T^,V^] ~ x V^. Without normalization 
factors, these tensor products would contain square roots. 

An identity used in this paper is the recoupling of the 
product of two inner products to the inner product of 
two tensors: 

l+V 

(a'.b')(A''.B'')^ ^ [x]nl,A^\A''Y -[h^B^'r 

x=\l~V\ 

(A3) 

This recoupling allows us to combine similar quantities 
to a total tensor of rank x. 
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Appendix B: Table of symbols 



Symbol 


Significance 


Iff) 


Ground state of the system 




\n) 


Intermediate state of RIXS process 




\f) 


Final state of the system 




Eo 


Ground state energy of the system 


s 


En 


Intermediate state energy of the system 


s 


Ef 


Final state energy of the system 


s 


|g) = l5;ke) 

|f) = l/;kV) 


ground state of system plus incident photon 




final state of system plus scattered photon 




Es 


= Eg + ftOJk 


s 


Ei 


= Ef + ftcjk/ 


s 


r„,r 


Inverse lifetime of intermediate state 


s 


Kk / hk' 


momentum of incident/scattered photon 
transferred momentum hk' — hk 


V 




V 


K 


reduced wave vector m = q — nG 


V 


G 


Reciprocal lattice vector 


V 




Incident/scattered photon energy, k = fck 


s 


hw 


Transferred energy hw^ — huj^^i 


s 


e / e' 


Polarization of incident /scattered photon 


V 


G / TT 


Polarization vector is perpendicular (a) 
or parallel (tt) to the scattering plane 




29 


Scattering angle (between k and k') 


s 


R 


Radius of curvature of grating 


s 


I 


Double-differential RIXS cross section 


s 


w 


Transition rate 


s 


dfl 


Solid angle of scattered photons 


s 


^fg 


RIXS scattermg amplitude 


s 


r 


Fundamental RIXS scattering amplitude 


t 




KlXh scattering amplitude 
in components of Vq 


t 


J. 


Angular and polarization dependence of 


t 




Resonance function 


s 


vv 


Effective transition operator 





X 


rank of tensor of fundamental 
spectra and operators 


s 


- clL 


Core to valence level orbital 





E 


Average energy 


s 


Mm 


H-H nfbitfil i"Vfln<^if inn nnprfifnr 


Q 


—9 ' "9 


d-d spin-orbital transition operators 







Is — > 4p dipole transition amplitude 


s 




Normalization constant 


s 




= [lL]nli^k' /[l]\\ 


s 


[X] 


=2x+l 


s 


[xiX2 . . . X„] 


= (22:1 -f- 1){2X2 + 1)... {2xn + 1) 


s 


[,r 


Tensor product, see Appendix IaI 


t 


Ylm 


Spherical harmonics 


s 


Rni 


Radial wave functions 


s 


ji 


Spherical Bessel functions 


s 




= r'f(",f(^) = V47r/HY,„(f) 


t 


/A B C\ 


3j symbol 




\a b c J 


s 



r n 


Electron position, momentum 


V 


t 


Time 


s 


R 


Ion position 


V 


A 


Electromagnetic vector potential 


V 


<f> 


Electric Potential 


s 


E, B 


Electric, magnetic field 


V 


a 


Vector of Pauli matrices: cr = ((Ta^, ct„, Uz) 


V 


p 


Charge density 


s 


H 


Hamiltonian 





G(z) 


Green's function 







^ Eg + hui, + ir 


s 


/ a 


Photon creation/annihilation operator 





V 


Electromagnetic transition operator 





D 


Electromagnetic dipole operator 


V 


Vf 


Fermi velocity 


s 


A 


Charge transfer energy 


s 


t,t',t" 


l"", 2"^^ and 3'''^ neighbor 






hopping amplitudes 


s 




On-site d-d Coulomb energy 


s 




Core hole-valence electron Coulomb energy 


s 


J 


Superexchange constant 


s 




Dynamic structure factor 


s 


3d" 


3d levels contain n electrons 




L 


Hole at the ligand site 




m 


Mass of electron 


s 


-e 


Charge of electron 


s 


c 


Speed of light in vacuum 


s 


h 


Planck's constant 


s 


ao 


Bohr radius 


s 


a 


Fine structure constant /A-jreohc 
Classical electron radius e^/(47reomc^) 


s 


Te 


s 


Z 


Atomic number 


s 


eo 


Vacuum permittivity 


s 


V 


Volume 


s 


iV 


Number of electrons 


s 



TABLE IV Table of symbols used in this article, the physical 
quantity they denote and if they are used as a scalar (s), 
vector (v), tensor (t) or operator (o). 
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